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I.  Introduction 

Theoretical  techniques  for  describing  laser-induced  molecular  dynamics 
combined  with  multiphoton  processes  are  reviewed,  with  emphasis  on  our  own 
work  carried  out  during  the  past  five  years.  The  presentation  is  divided 
according  to  Sections  II  and  III  as  Gas-Phase  Processes  and  Surface  Processes, 
respectively.  The  former  section  is  further  divided  into  two  parts,  where 
Part  A  is  Unimolecular  Dynamics  and  Part  B  is  Molecular  Collisions.  Within 
Part  A  we  address  five  topics:  1.  heat-bath  models  (both  classical  and 
quantum  mechanical);  2.  the  rate  equation  approach  using  a  quasi-continuum 
model;  3.  the  excitation  of  intramolecular  modes;  4.  selectivity  in 
multi  photon  processes;  and  5.  multiple-frequency  laser  excitation.  Within 
Part  B  we  address  two  topics:  1.  resonance  formation,  chemical  reactions 
and  transition-state  spectroscopy;  and  2.  isotopic  selectivity. 

Section  III  is  also  divided  into  two  parts,  where  Part  A  is  Energy 
Flow  in  Adspecies-Surface  Systems,  and  Part  B  is  Laser  Applications  to  Sur¬ 
face  Chemistry.  Part  A  constitutes  most  of  Section  III  and  is  further  di¬ 
vided  into  five  topics:  1.  microscopic  treatment  of  single-phonon  and  multi - 
phonon  processes;  2.  selective  excitation  and  thermal  effects,  including 
(a)  the  use  of  the  memory  function  to  describe  feedback  mechanisms  from  the 
substrate  to  the  adspecies,  (b)  the  competition  between  selective  and  nonse- 
lective  processes  in  a  multilevel  system,  and  (c)  adbond  excitation  in  the 
hydrogen/tungsten  system;  3.  photon  energy  population;  4.  isotope  effects; 
and  5.  dynamic  effects  associated  with  migration  and  desorption. 


II.  Gas-Phase  Processes 


A.  Unimolecular  Dynamics 

Since  the  first  report  of  infrared  laser  dissociation  of  polyatomic 

molecules  in  1971,*  the  phenomena  of  multiphoton  excitation  (MPE)  and 

dissociation  (MPD)  have  been  extensively  studied  and  have  been  shown  to  be 
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a  novel  method  for  vibrational -mode-controlled  infrared  photochemistry. 

Most  of  the  major  aspects  of  MPE  and  MPD  are  now  qualitatively  understood 
by  experimental  measurements  and  related  theoretical  studies.  However, 
the  quantitative  understanding  of  multiphoton  processes  and  their  appli¬ 
cations  to  selective  photochemistry  is  still  in  its  infancy.  The  complexity 
of  the  aforementioned  processes  results  not  only  from  the  coherent  properties 
of  the  laser  radiation  but  also  the  states  of  the  excited  molecular  system. 
Some  of  the  important  questions  which  have  been  (or  are  expected  to  be) 
addressed  are  the  following: 

(i)  For  a  large  molecular  system,  what  is  the  nature  of  the  excited 
vibrational  state  and  how  does  it  change  with  laser  parameters,  such  as 
intensity,  fluence  (intensity  x  time)  and  frequency? 

(ii)  What  is  the  magnitude  of  the  absorption  cross  section? 

(iii)  What  is  the  interplay  between  coherent  and  incoherent  processes 
when  the  molecular  system  is  excited  and  makes  the  transition  from  discrete 
states  to  the  so-called  quasi-continuum? 

(iv)  What  is  the  role  of  intramolecular  vibrational  relaxation  (IVR) 
in  MPE,  and  what  determines  the  rate  of  IVR?  Is  it  always  fast  enough 

to  randomize  the  energy  absorbed  by  the  intramolecular  modes,  so  that  a 
a  statistical  theory  may  be  used  to  describe  MPD? 
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(v)  How  is  the  degree  of  isotopic  selectivity  determined  by  the  IVR 
rate,  the  laser  pumping  rate  and  the  dissociation  rate? 

(vi)  What  are  the  dynamics  of  the  dissociation  event? 

(vii)  How  fast  is  the  randomization  of  intramolecular  energy 
compared  to  a  chemical  process  such  as  fragmentation? 

Quantitative  answers  to  the  above  questions  will  provide  us  with  an 
understanding  of  the  photophysical  (excitation  and  relaxation)  as 
well  as  the  photochemical  (decomposition  and  reaction  of  fragments) 
processes.  However,  the  theory  developed  so  far  and  the  experimental 
information  are  not  sufficient  for  us  to  construct  these  answers.  We 
therefore  note  that  the  discussions  in  this  Part  A  of  Section  II 
only  provide  either  qualitative  or  semi -quantitative  descriptions  of  the 
above  features  of  MPE  and  MPO. 

Two  distinctly  different  approaches  will  be  discussed:  (i)  heat- 
bath  model  (HBM)  and  (ii)  quasi-continuum  model  (QCM).  In  the  first 
approach,  HBM,  the  active  mode  of  the  molecular  system  is  singled  out, 
and  the  coupling  between  this  mode  and  the  bath  modes  provides  the  relaxa¬ 
tion  mechanism.  With  a  consequent  energy  leakage  from  the  active  to  the 
bath  modes,  the  molecule  is  heated  up  and  the  increase  of  the  mode-mode 
coupling  results  ina  red  shift  and  smearing  out  of  the  absorption  spectrum 
of  the  pumped  mode.  The  second  approach,  QCM,  is  based  on  a  picture  in 
which  the  polyatomic  molecule  forms  a  quasi-continuum  (QC)  at  high  energies 
due  to  its  high  density  of  states.  In  this  approach,  all  vibrational  modes 
of  the  system  are  treated  on  equal  footing,  and  each  individual  state  of 
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the  QC  is  generally  a  superposition  of  all  the  modes.  To  describe  the 
dynamics  of  MPE  based  on  HBM,  we  shall  start  with  a  simple  classical 
model  and  a  generalized  Langevin  theory.  Then  the  results  of  a  quantum 
treatment  will  be  discussed  and  compared  with  those  of  the  classical  approach. 
The  second  approach,  QCM,  which  involves  the  rate  (master)  equation  for 
the  energy  populations,  will  be  used  to  describe  the  phenomena  of  MPO, 

To  demonstrate  the  dynamics  of  MPE,  a  quantum  treatment  of  IVR  will  be 
discussed  and  the  rate  of  IVR  for  SFg  molecule  will  be  estimated. 

Different  types  of  selectivity  and  the  possibility  of  bond-selective 
photochemistry  will  be  investigated.  Finally,  the  cooperative  effects 
in  MPE  via  multiple-frequency  lasers  will  be  analyzed  based  on  both  a 
quantum  treatment  and  a  classical  treatment. 
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1.  Heat-Bath  Models 

Both  classical^"^  and  quantum^"^  treatments  have  been  employed 
to  calculate  the  photon  energy  deposited  in  a  damped  anharmonic  classical/ 
quantum  oscillator.  The  total  energy  of  a  laser-excited  molecular  system 
may  be  described  either  by  a  classical  Lagrangian  or  a  quantum  Hamiltonian, 
while  the  energy  transfer  dyanamics  is  governed  by  Newton's  or  Heisenberg's 
equations  of  motion.  The  absorption  cross  section  may  be  obtained  from  the 
classical  power  absorption  (force  x  velocity)  or  from  the  quantum  rate  of 
excitation  (expectation  value  of  the  energy  operator).  Immediately  below, 
classical  and  quantum  models  are  shown,  and  the  results  will  be  compared 
later  in  Part  A. 5  of  this  section,  where  the  absorption  cross  section 
of  two-laser  excitation  will  be  investigated. 

a.  Classical  Models5"8 

The  classical  Lagrangian  describing  a  polyatomic  molecule  subject  to 
infrared  radiation  may  be  written  in  terms  of  the  normal  coordinates 
as  follows: 

L(Q1,Q2,...Qi)  =  i-0  -  V  +  f(t)QA,  (II. 1) 

where  Lq  is  the  unperturbed  Lagrangian,  and  V  is  the  interaction  potential 
of  the  molecular  system  in  which  one  of  the  normal  modes,  QA,  is  infrared 
active  and  coupled  to  the  laser  driving  force  f(t).  The  corresponding 
Newton's  equations  of  motion  are 
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i  —  =  0.  (II. 2) 

dt  3Q1  3Qi 

For  a  given  potential  energy  V,  the  equation  of  motion  for  the  pumped  mode 

(normal  coordinate)  may  be  derived  by  Eq.  (II. 2).  However,  the  solution 

of  that  is  in  general  not  available  due  to  the  set  of  coupled  equations 

whose  dimensional ity  is  proportional  to  the  numbers  of  normal  modes  in  the 

system.  It  is  possible  that  for  a  simple  functional  form  of  V,  e.g., 

linear  coupling  case  involving  terms  like  z  Q.Q. ,  the  many-body  problem 

i  1 

of  the  system  may  be  reduced  to  a  single-body  problem  in  which  the 
ensemble-averaged  equation  of  motion  of  the  active  mode  coordinate  is 
given  by,8  letting  QA  =  Q, 

••  o 

<Q>  +  Y<0>  +  Wg-f-f  <Q>  =  6E  cos(o)t)/m.  (II. 3) 

<...>  denotes  the  ensemble  average  over  the  bath-mode  coordinates;  y  is 
a  damping  factor  resulting  from  the  active-  and  bath-mode  interactions; 

E  is  the  electric  field  of  the  IR  radiation  linearly  polarized  in  the 
direction  of  QA  with  a  circular  frequency  to;  m  and  e  are  the  reduced  mass 
and  the  effective  charge  of  the  active  mode  which  is  treated  as  an  anharmonic 
oscillator  with  an  effective  frequency  given  by8 

Ugff  s  “0  ”  *  ( II  *4  ) 

where  <og  is  the  harmonic  frequency  of  the  active  oscillator  and  is  red 
shifted  by  the  anharmonic  term  proportional  to  the  anharmonicity  K*  and 
the  steady-state  amplitude  A. 
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By  ignoring  the  transient  damping  part  of  the  solution  of  Eq.  (II. 3), 
the  steady-state  solution  gives  us  the  instantaneous  power  absorption  of 
the  oscillator  defined  as  force  times  velocity.  For  near-resonance 
excitation,  the  time-  and  ensemble-averaged  power  absorption  is  found  to 


(a-K*A2)2+y2 


(II. 5) 


which  is  a  Lorentzian  with  FWHM  =  2y.  The  above  expression  provides  us 
with  a  classical  absorption  cross  section  a  =  <dE/dt>/I,  where  I  is  the  laser 

O 

intensity  related  to  the  electric  field  by  I  =  E  c/8n,  c  being  the  speed 
of  light.  The  steady-state  average  excitation  (numbers  of  photon  absorbed) 
of  the  active  mode  defined  by  X  =  mu)A2/2ti  =  <dE/dt>/(vM  corresponds  to 
the  roots  of  the  cubic  equation 


aX3  +  bX2  +  cX  +  d  =  0, 

(II. 6a) 

a  *  (K*/mw0)2  =  (K*)2, 

(II. 6b) 

b  =  -2aK* 

(II. 6c) 

c  *  (A2+r2)/(K*)2 

(II. 6d) 

d  =  -  e2E2/(4mfkj). 

(II. 6e) 

[The  coefficient  "c"  in  Eq.  (II. 6d)  is  obviously  not  the  speed  of  light.] 
We  note  that  the  above  expression  shows  an  asymmetric  Lorentzian  when 
K*  f  0  and  there  is  an  optimal  detuning  a*  =  K*X*,  where  X*  is  the  maximum 
value  of  X  given  by  (dX/dA)  =  0.  Furthermore,  for  a  fixed  laser  intensity. 
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the  excitation  profiles,  X  vs.  a,  show  the  "bistability"  feature  given 
by  the  criterion  K*>K*W  =  y^/(2d)^^. 

We  have  shown  that  the  absorption  cross  section  is  governed  by  a 
Lorentzian,  Eq.  (II. 5),  where  a  constant  damping  factor  (or  the  energy 
relaxation  rate),  y»  is  assumed.  For  a  realistic  system,  the  damping 
factor  could  be  time  dependent,  and  the  phase  of  the  oscillator  cuuid 
also  relax  due  to  the  incoherent  behavior  of  the  excited  molecular  system. 

For  a  more  rigorous  description  of  the  dynamics  of  MPE,  instead  of  the 
simple  equation  like  Eq.  (II. 3),  we  consider  the  generalized  Langervin 
equation  (GLE)7 

t  - 

Q  +  /  dt ' [K, (t-t * )  +  «fffK2(t-t,)]G(t')  =  [f (t)  +  R(t)]/m.  (II. 7) 

0 


Here  we  have  introduced  the  energy  damping  kernel  and  the  dephasing 
kernel  IC,.  f ( t)  and  R(t)  are  the  laser  driving  force  and  the  bath- 
induced  random  force.  Considering  a  Markoff  process  with  K^(t)  =  2y16(t) 
and  an  exponential  dephasing  kernel  Kg(t)  =  exp(-2Y2t),  we  find,  from  the 
Laplace-Fourier  transform  of  the  velocity  autocorrelation  function,  the 
average  energy  absorption  rate  to  be  (assuming  that  <R(t)>  =  0) 


<I>  -  <<*)2p<t0) 


YgC+uiD 

— 2 — 1 

c  +u 


(II. 8a) 


P(T0)  =  (2ha))"1(kT0/m)Cl-exp(-WkT0)], 
2  2 

C  =  -  U)  +  Y1Y2* 

D  =  w(y1  +  y2)> 


(II. 8b) 
(II. 8c) 


(II. 8d) 


where  Tg  is  the  initial  temperature  of  the  molecule.  For  the  situation  of 
2  2  ? 

near- resonance,  -  w  s  2<dA  and  YjY2<<w  »  Eq.  (II. 8)  becomes 


dE 

*  ' 

eE 

2 

p<v 

2 

(a-K*A^)^+(y2+Y2^/^ 

Since  at  steady  state  the  rate  of  energy  deposited  in  the  active  mode 
equals  the  rate  of  energy  relaxation,  the  above  expression  also  repre¬ 
sents  the  rate  of  energy  absorbed  by  the  whole  system  (active  plus  bath 
modes).  By  conservation  of  energy,  we  then  readily  obtain  the  rate 
equation  for  the  active  mode  energy  as 


d<E.> 

~at“  =  aI  ‘  yi(<ea>-<eb>)» 


(II. 10) 


where  ol  a  <dE/dt>  given  by  Eq.  (II. 9)  is  the  energy  absorption  rate  of 
the  whole  molecule  with  absorption  cross  section  a  and  laser  intensity 
I,  and  <E^>  and  <Eg>  are  the  ensemble-averaged  energies  of  the  active  mode 
and  the  bath  modes.  We  note  that  in  this  classical  rate  equation  the 
net  energy  transfer  rate  of  the  active  mode  is  characterized  by  the 
energy  relaxation  rate  (y^)»  while  the  total  absorption  cross  section 
of  the  molecule,  ol,  is  governed  by  the  total  relaxation  rate,  y1  +  y2* 

We  shall  show  later  that  the  above  derived  classical  rate  equation  can 
also  be  derived  from  a  microscopic  quantum  Hamiltonian. 

In  addition  to  the  above  phenomenological  models,  a  classical 

g 

trajectory  calculation  of  MPE  was  recently  made  for  SFg.  The 
results  show  that  the  energy  transfer  and  dissociation  rate  depend 


on  laser  fluence  (intensity  x  time),  and  the  intramolecular  energy 

relaxation  rate  is  estimated  to  be  on  the  order  of  a  picosecond. 

Another  classical  treatment  of  MPE  and  MPD  for  a  system  of  two 

nonlinearly  coupled  oscillators  was  studied  by  means  of  Krylov- 
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Bogol iubov-Mitropolsky  theory.  The  results  show  two  regions  of 
behavior  for  the  exchange  of  energy  between  the  system  and  the  laser 
field:  (i)  the  regular  region  with  well-defined  frequencies  and  (ii) 
the  erratic  region.  Motion  in  the  latter  region  leads  to  dissociation 
of  the  pumped  molecule.  For  a  single  oscillator,  the  motion  is  always 
periodic,  while  for  a  system  with  many  oscillators,  e.g.,  a  polyatomic 
molecule,  the  motion  of  the  system  becomes  ergodic  when  it  is  highly 
excited  due  to  the  strong  anharmonic  coupling  among  the  intramolecular 
modes. 

b.  Quantum-Mechanical  Models 

In  the  quantum  heat-bath  treatment,  the  MPE  of  a  polyatomic  molecule 
may  be  described  by  the  total  Hamiltonian 

H  =  HA  +  HB  +  HAF  +  HAA  +  HBB  +  HBF*  (u* 

H^  g  represent  the  unperturbed  Hamiltonians  of  the  active  (A)  mode 
and  of  the  bath  (B)  modes  with  the  interaction  Hamiltonian  HAg; 
represents  the  A  mode  and  laser  field  interaction;  and  the  last  three 
terms,  which  are  often  neglected  in  theoretical  studies,  represent  the 
mode-mode  couplings  among  the  active  modes  ( HAA ) ,  the  bath  modes  (Hgg) 
and  the  direct  laser  excitation  of  bath  modes  (Hgp).  For  an  isolated 


molecule  or  collisionless  process,  is  negligible.  The  mode-mode 
coupling  among  the  bath  modes  usually  is  very  strong  such  that  the  time 
scale  for  reaching  a  thermal  equilibrium  among  these  modes  is  much 


faster  than  that  of  the  active-mode  excitation  and  relaxation.  There¬ 


fore,  the  rapid  dynamics  within  the  bath  modes  governed  by  Hgg  may  be 
absorbed  into  Hg  in  which  the  occupation  number  of  these  modes  is 
treated  as  a  constant.  Finally,  the  direct  excitation  of  the  bath 
modes  through  Hgp  is  not  significant  when  their  frequency 
spectrum  is  far  off-resonant  with  the  laser  frequeny.  Within 


the  spirit  of  the  heat-bath  treatment,  we  assume  that  the  active  mode 
can  always  be  singled  out  from  the  other  modes,  which  serve  as  an  infinite 
energy  sink  defined  by  a  constant  temperature.  However,  we  may  have  the 
situation  that  the  IVR  rate  is  very  fast,  particularly  when  the  molecule 
is  highly  excited,  such  that  the  photon  energy  is  rapidly  randomized  to 
all  the  bath  modes.  For  a  finite  heat  bath,  we  expect  to  have  energy 
feedback  from  these  thermal  bath  modes  to  the  active  mode.  Furthermore, 


when  the  active  mode  is  highly  excited,  its  vibrational  frequency  may  be 
comparable  to  that  of  the  bath  modes,  and  the  total  absorption  cross 
section  of  the  molecule  reflects  both  H^p  and  Hgp.  In  the  above 
described  situation,  the  heat-bath  mode  model  may  not  be  an  appropriate 
treatment  of  MPE  processes.  A  more  relevant  treatment  based  on  an 
"equal  footing"  of  all  the  modes  will  be  discussed  later.  In  the 
following  discussion,  we  shall  neglect  the  last  three  terms  in  Eq.  (II. 11) 
and  focus  on  the  excitation  of  the  active  mode  while  treating  the  bath 


modes  as  an  infinite  sink. 


The  second  quantization  expressions  of  the  components  of  the 
Hamiltonian  of  Eq.  (II. 11)  are  given  by12,14 

Ha  =  fl(co0-e*a+a)a+a, 

HB  ■  ?*“jbjbr 

J 

H.n  =Tiy(G  B  a+  +  G*B+a), 

AB  L  v  v  v  v 

v 

HAF  =  TiV  cos(uit)  (a+  +  a) 

V  =  (2fima,0)‘1/2w'(0)E 

4* 

a  ,a  and  b.,  b.  are  the  usual  harmonic  vibrational  ladder  operators  for 
J  J 

the  active  mode  and  the  j-th  bath  mode  with  frequencies  u)Q  and  , 
respectively;  e*  is  the  anharmonicity  of  the  nonlinear  quantum  oscillator 
(active  mode)  with  the  derivative  of  the  diple  moment  u'(0)  evaluated  at 
the  equilibrium  point;  and  B^  =  nbj  represents  a  multimode  operator  with 
an  order  p  defined  by  p  s  oig/uj. 

In  the  Heisenberg  picture,  the  equations  of  motion  for  the  operator 
0  are  given  by 

-  W 

for  0  =  a+,  a,  B^,  Bv  etc.  Here  b{ t)  =  exp(iHt/h)0  exp(-iHt/fi)  is  a  time 
dependent  Heisenberg  operator  of  the  Schrodinger  operator  0  which  is  time 
independent.  In  general,  tractable  solutions  of  the  above  coupled 
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equations  are  not  available  due  to  the  many-body  effects  and  the  nonlinear 
couplings  of  the  bath  modes.  However,  by  employing  the  so-called  Weiskoff- 
Wigner  single-pole  approximation  or  assuming  a  Markoff  process  for  the  bath 
modes,  the  many-body  effects  resulting  from  H^g  may  be  reduced  to  a  damping 
factor  and  a  frequency  shift  for  the  active  mode.  The  above  treatment  is 
similar  to  the  usual  phenomenological  approach,  where  a  complex  frequency 
of  the  active  mode  is  assumed,  except  that  the  damping  factor  may  be 
rigorously  derived  based  on  the  microscopic  Hamiltonian.  Beyond  this 
heat-bath  model,  the  actual  mode-mode  coupling  should  be  rigorously 
treated  without  assuming  a  continuum  spectrum  for  the  bath  modes.  This 
will  be  investigated  later  in  the  intramolecular  mode-mode  coupling  within 
the  SFg  system. 

Within  the  heat-bath  model,  we  obtain  a  set  of  equations  of  motion 
for  the  active  mode  operator  and  the  average  excitation  n  =<a+a>  in  the 
rotating-wave  approximation, 

(II. 14a) 

(II. 14b) 

<<...>>  represents  the  ensemble  average  over  the  active-  and  bath¬ 
mode  coordinates;  Yj  2  are  the  bath-mode- induced  energy  and  phase 
relaxation  factors;  "N  is  the  multimode  occupation  number  of  the  bath 
modes;  and  all  other  parameters  are  as  previously  defined.  In  the  adia¬ 
batic  approximation,  [ a  1  (>'i+'>'2 ) >:> I ®  1  *  we  °btain  the  rate  quation 
for  the  active-mode  excitation. 


d<a> 

dt 


i(A-2e*<a+a>)  +  (y1+y2)/2  <a>  -  iV/2, 


d<n> 

dt 


=  -  (iV/2)<<a+-a>>  -  Y^(<f1>  “  N). 
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d<n>  _  al 
dt  fill 


Y1(<n>-N). 


(11.15) 


This  is  seen  to  be  the  quantum  version  of  the  classical  rate  equation 
(II. 10)  by  noting  that  <E^>  = 'hoxn>  and  <Eg>  =  fiiJT.  The  quantum 
absorption  cross  section  of  the  molecule  is  given  by 


Using  Eq.  (11.16),  a  cubic  equation  for  the  steady-state  excitation 
X  =  <n(t)>  is  readily  found  from  Eq.  (11.15), 

5  •  5  • 


AKI+Tj/Yj) 

X - 5 - +  N, 

(a-2e*X)  +(y1+y2)  /4 


(11.17) 


which  is  identical  to  the  classical  results  in  Eq.  (II. 6)  [or  Eq.  (II. 9) 
with  y2  included]  if  we  ignore  the  initial  heat-bath  occupation  number 
IT  and  make  the  corresponding  relations:  2e*  =  K*  and  AI  =  d.  The  above 
results  show  the  power  law  for  the  steady-state  excitation  vs.  laser 
intensity  X  «  Ip,  with  p  *  1  in  the  low-excitation  harmonic  region, 
2e*X<<(y1+Y2)/2»  and  1/3  <  p  <  1  in  the  high-excitation  anharmonic 
region.  Furthermore,  from  the  time  integral  of  the  total  cross  section 
of  the  molecule,  ol/fia),  we  obtain  the  power  law  for  the  mean  numbers  of 
photons  (<N>)  absorbed  by  the  molecule  (active  mode  plus  bath  modes)  vs. 
laser  fluence  (41)  as 

<N>  =  Jdt(al/fiu))  «  4>^» 

4>  =  /dt  I  , 


(II. 18a) 
(II. 18b) 
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where  the  nonlinear  power  q<l  results  from  the  anharmonicity  2e*<n>  of 
the  absorption  cross  section.  The  calculated  power  law  may  be  used  to 
fit  the  experimental  results  in  which  the  average  excitations  are  deduced 
from  the  direct  attenuation  of  the  laser  beam  propagating  through  the 
SFg  gas  cell.*^®)  In  Fig.  1,  we  see  that  in  the  low-excitation  region 
the  data  are  well  fitted  by  Eq.  (11.18)  with  q  =  1,  whereas  in  the  high- 
excitation  region  the  curve  is  bent  to  a  power  law  with  q  =  0.65.  An 
alternative  measurement  of  <N>  vs.  laser  intensity,  rather  than  the 
fluence,  is  shown  in  Fig.  2,  which  indicates  the  power  law  of  <N>  <*  Iq 
with  q<l.15(b) 

From  the  above  calculated  results  based  on  the  heat-bath  model, 
the  following  features  of  MPE  may  be  summarized: 

(i)  The  steady-state  excitation  of  the  active  mode  is  governed 

by  the  laser  intensity,  rather  than  the  laser  fluence,  by  a  power 
law  <n>  «  Ip  with  1/3  <p<l  . 

(ii)  When  the  excitation  of  the  active  mode  is  saturated,  i.e., 
the  pumping  rate  is  balanced  by  the  relaxation  rate,  the  number  of 
photons  absorbed  by  the  molecule  (active  plus  bath  modes)  is  characterized 
by  the  laser  fluence  with  a  power  law  <N>  «  $q  with  q«l. 

(iii)  Mode-selective  excitation  is  possible  by  a  high-power  short- 
pulse  laser  in  which  the  excitation  rate  of  the  selectively-pumped  mode 

is  comparable  or  faster  than  the  relaxation  rate,  i.e.,  the  energy  leakage 
rate  from  the  pumped  mode  to  the  bath  modes. 

(iv)  For  long-pulse  laser  excitation  of  a  system  with  very  strong  mode¬ 
mode  coupling,  the  photon  energy  is  rapidly  randomized  to  all  the  modes, 
and  an  effective  vibrational  temperature  of  the  excited  molecule  is  defined 
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by  the  conservation  of  energy  as  follows: 

Sflm  =  <N>fiu, 

where  S  is  the  total  number  of  vibrational  modes  (e.g.,  S  *  15  for 
the  SFg  molecule),  and  7^  is  the  overall  occupation  number  of  the 
molecule  given  by  the  Bose-Einstein  function 

Nm  =  [exp(H«/kTeff)  -  l]’1. 

Here  all  the  intramolecular  modes  are  assumed  to  be  in  thermal 
equilibrium  with  an  effective  temperature  Tg^  and  a  mean  vibrational 
frequence  uj.  We  note  that  the  average  molecular  excitation  <N(t)> 
given  by  Eq.  (11.18)  in  general  is  time  dependent  through  the 
temporal  profile  of  the  laser  pulse,  I(t).  Therefore  Tfi^(t)  is  also 
time  dependent  and  reaches  a  peak  value  at  the  end  of  the  pulse,  e.g., 
a  square  pulse,  and  gradually  decays  to  the  initial  temperature  after 
the  pulse  is  off.  Further  discussion  of  will  be  given  later 
when  a  thermal  population  is  introduced  by  the  solution  of  the  rate 
equatiun. 

The  heat-bath  model  may  be  extended  to  a  system  consisting  of 
several  subgroups  in  which  the  mode-mode  coupling  is  more  significant 
within  each  subgroup  than  between  the  subgroups.  In  this  situation 
the  specifically-excited  subgroup  may  be  characterized  by  a  higher 
T  ff  than  the  others  and  in  turn  opens  the  selective  channels  of 
laser-induced  rate  processes  which,  e.g.,  are  governed  by  an 
Arrhenius  form  exp(-EA/kT  ff) ,  where  EA  is  an  activation  energy  of 
the  specific  channel.  We  note  that  the  rate  processes  are  determined 


(II. 19a) 


(II. 19b) 
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by  the  ratio  of  E/^AT  and  not  simply  by  EA,  so  that  the  "easier  channel" 
is  not  necessarily  the  one  with  a  lower  E^. 

2.  Quasi-Continuum  Model  —  Rate  Equation  Approach 

We  now  discuss  the  second  quantum  approach, with  treats  all  vibrational 
modes  on  an  equal  footing,  with  the  eigenstates  of  the  system  being  mixtures 
of  all  the  normal  modes.  Due  to  the  high  density  of  vibrational  states  at 
high-excitation  energy,  the  molecule  can  easily  absorb  more  photons 
through  resonant  incoherent  transitions  between  so-called  quasi -continuum 

1C  po 

levels,  where  rate  equations  have  been  widely  used  for  collisionless  and 
collisional24-26  MPD. 

A  well -accepted  qualitative  model  for  MPD  of  SFg  has  provided  a 
general  picture  where  the  molecular  energy  levels  are  divided  into  three 
regions  (see  Fig.  2):  (I)  the  discrete  region  characterized  by  coherent 
excitation,  where  the  v3  (active)  mode  absorbs  3  to  6  photons  and  the 
anharmonicity  of  the  vibrational  potential  is  nearly  compensated  for  by 
allowed  rotational  transitions; (II)  the  quasi-continuum  region  where  the 
level  density  is  very  high  and  incoherent  excitation  is  essential;  and 
(III)  the  true  continuum  region.  Region  (I)  processes  are  responsible 
for  isotopic  selectivity,  coherent  effects  (multi photon  resonance, 
photon  echoes,  coherent  wave  propagation,  etc.)  and  the  intensity 
dependence  of  MPE  with  high  selectivity.  For  the  excitation  processes 
in  regions  (II)  and  (III)  it  has  been  shown  experimentally  for  SFg  that 
high  laser  fluence  (energy),  not  high  laser  power  (intensity),  is 
necessary  for  driving  the  molecule  through  the  quasi-continuum  and 
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is  the  important  parameter  for  determining  the  dissociation  yield. ^ 

However,  later  experiments  on  other  molecules  have  shown  that  both 

2(c) 

laser  fluence  and  intensity  determine  the  dyanmics  of  MPD.  '  ' 

For  multiphoton  excitation  in  the  quasi-continuum  region,  the 
Fermi  Golden  Rule  is  valid  and  the  full  Schrodinger  equation  reduces 
to  a  set  of  incoherent  rate  (master)  equations.  For  the  collision- 
free  condition  (low  pressure),  the  rate  quation  for  the  energy  popu¬ 
lation  is 


Ur 

-r£-  -  Wa  ,Pn  .  +  .  -  (W*  +  W®  ,)Pn  -  kP, 

dt  n-1  n-1  n  n+1  n  n-1  n  n  n 


(11.20) 


where  Pp  is  the  population  in  the  n-th  level  (n  photons  absorbed). 

Wa  (We)  is  the  transition  rate  constant  for  absorption  (emission) 

from  level  n  to  n+1  (n+1  to  n)  and  is  related  to  the  absorption 

cross  section  (an)  and  frequency  of  the  field  (with  intensity  I) 

by  Wa  =  onI/fia)  and  W®/Wa  =  9n/9n+1>  where  gn  is  the  molecular 

density  of  states  at  energy  rrfiu>.  kn,  the  unimolecular  decomposi- 

27 

tion  rate  constant,  can  be  calculated  by  RRKM  theory  or  quantum 
RRK  theory.28  From  Eq.  (11.20)  we  realize  that  collisionless  MPE 
and  MPD  are  characterized  by  the  laser  intensity  and  frequency,  the 
absorption  cross  section,  the  density  (or  degeneracy)  of  states  and 
the  unimolecular  reaction  constant.  During  the  past  several  years, 
the  rate  equations  describing  MPD  have  been  studied  by  different 
approaches,  such  as  a  thermal  model  for  Boltzmann-type  energy 
populations,2^  diffusion  model  for  continuum  populations,22  exact 

25 

stochastic  model,  the  model  of  restricted  intramolecular  relaxation 

29 


and  the  random  coupling  model. 
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Depending  on  the  form  of  the  absorption  cross  section  and  the 
density  of  states  in  the  quasi -continuum,  the  solution  of  the  above 
rate  equation  gives  us  different  populations  and  corresponding  dis¬ 
sociation  probabilities.  For  example,  various  types  of  n-dependence 
of  the  absorption  cross  section  have  been  considered:  an  =  oQ/n, 
an  =  exp(-an+e)  and  on  =  constant  in  Refs.  19,  20  and  16,  respectively. 
A  "thermal  model"  with  suitable  functional  forms  for  o  and  q  (and 

neglecting  the  dissociation  rate  kn)  has  been  studied  in  which  the 

21 

population  is  governed  by  the  rate  equation 


dPn  AI 

nr  *  f  + 


T  '"+1)Pn 


(2n+S), 


(11.21) 


where  the  constant  A  is  proportional  to  the  absorption  cross  section  and 
S  is  the  number  of  vibrational  modes  involved  in  the  process.  The  above 
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rate  equation  can  be  solved  in  closed  form  to  give  the  thermal  population 


pn  =  N9n  exp[-rrfWkTeff], 


(11.22) 


where  N  is  a  normalization  constant  and  gn  is  the  quantum  degeneracy  given  by 


N  -  IP 


n=0  n 


[1  -  exp(fiu/kTeff]' 


(II. 23a) 


9n  -  (n+S-l)!/[n!(S-l)!], 


(II. 23b) 


and  the  effective  temperature  T  ^ 
conservation  of  energy 


as  a  function  of  time  is  given  by  the 


n  =  [exp(-fWkTeff )  -  1]  ,  ( 

where  o  is  the  absorption  cross  section  and  $  is  the  fluence  of  the 
laser  radiation.  For  multiphoton  absorption  processes,  for  kTe^»fitu 
and  for  n»S,  we  obtain  the  high-temperature  or  classical  limit  of  the 
population, 


s-i  n 

r 

n 

exp 

nS 

_(S-1) ! (<n>/S)S  _ 

<n> 

where  <n>  is  the  mean  numbers  of  photon  absorbed  per  molecule,  i.e., 

<n>  =  I  nP  .  We  note  that  in  this  classical  limit  kT  «  =  fiw<n(t)>/S, 
n=0  n  eTT 

which  follows  the  time  dependence  of  <n(t)>,  and  we  may  obtain  a  power  law 

according  to  Eq.  (11.18).  Furthermore,  the  corresponding 

dissociation  probability  is  found  to  be 


Pd  3  1  pn«  exp(-Sn*/<n>),  ( 

n=n* 

which  is  the  usual  Arrhenius  form. 

A  "continuum  model"  was  also  developed  for  more  general  functional 

forms  of  a„  and  g  : 

n  Jn 

an  *  °o(n+C)a'  ( 

9n  =  90(n+C)S  1 ,  ( 


where  Oq,  gg,  C  and  a  are  suitable  constants.  For  high  excitations  n>>l, 

the  difference  equation  (11.20)  reduces  to  a  classical  diffusion  equation 
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where  n  becomes  a  continuous  variable, 


fa 

3<j> 


-1  3  3P- 

^n  3n  ^nan  3n 


?)■ 


Using  the  expression  of  Eq.  (11.27),  the  population  can  be  solved  in  a 

closed  form,  and  the  corresponding  dissociation  probability  is  also 
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given  in  closed  form  as 


pd  *  «P(-T-)C1  +  I 

d  *  j-i  + 

♦0  s  (n*)e/(02ao), 


J  =  (S/B)-l,  B  -  2-o, 


for  integer  values  of  S/e.  This  model  results  in  a  set  of  "universal 
curves"  which  may  be  used  for  data  analysis  in  various  experimental 
situations.  We  note  that  the  thermal  results  shown  in  Eq.  (11.25)  repre¬ 
sent  a  special  case  of  Eq.  (11.29)  with  a  =  1,  i.e.,  a  harmonic 
model.  From  the  solution  of  Eq.  (11.28),  we  may  calculate  the  average 
excitation  (photon  quanta  per  molecule).  The  result  shows  a  power  law 

<n>«=$q,  where  q  =  (2-a)"*.  It  is  seen  that  <n>«4>  for  a  harmonic  model 

1/3  22 

(o=l)  and  <n>«$  '  for  an  anharmonic  model  (a=  -1).  These  results  are 

in  agreement  with  that  of  the  heat-bath  model,  Eq.  (11.17). 

The  dynamics  of  MPE  and  MPD  has  been  investigated  by  a  molecular 


beam  method,  and  the  experimental  results  have  been  analyzed  by  an 
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improved  rate  equation  where  the  intensity  dependence  of  excitation  over 

the  discrete  levels  is  included.  These  dynamical  features  are  summarized 

as  follows:  (i)  Both  the  laser  intensity  and  the  fluence  are  important  in 

determining  the  MPD  yield,  (ii)  Excitation  over  the  discrete  levels  into 

the  quasi-continuum  is  mainly  governed  by  the  laser  intensity,  (iii) 

Excitation  through  the  quasi-continuum  to  the  dissociation  level  depends 

only  on  the  laser  fluence.  (iv)  Above  the  dissociation  level,  the  excess 

energy  of  the  dissociated  molecule  is  determined  by  the  balance  between 

the  up-excitation  rate  and  the  dissociation  rate.  Therefore,  in  the  case 

of  long  laser  pulses,  the  excess  energy  is  determined  by  the  laser  intensity, 

while  for  short  pulses  case  it  is  determined  by  the  laser  fluence.  For  SFg, 

the  excess  energy  was  shown  to  be  laser  fluence  dependent  when  the  laser 

2 

pulse  has  energy  higher  than  5  J/cm  .  (v)  The  photon  energy  deposited  in 

the  active  mode  is  rapidly  randomized  among  all  the  vibrational  modes  on 
-8 

a  time  scale  of  10  s,  which  is  much  faster  than  the  dissociation  life¬ 
time.  The  observed  dissociation  rates  and  the  overall  dynamics  of  MPD  are 
in  good  agreement  with  the  results  predicted  by  RRKM  theory,  (vi)  The 
population  distribution  Pn  in  the  quasi-continuum  predicted  by  a  model 

with  a  *  exp(-an)  is  narrower  than  a  thermal  distribution.  Another  model 
n 

based  on  a  dephasing  broadening  in  the  quasi-continuum  also  shows  this 
narrower  population  which  is  governed  by  a  Poisson  distribution  (see 
Fig.  3).  More  recently,  a  model  of  quasi-continuum  absorption  has 
been  developed  to  provide  a  more  precise  meaning  for  the  term  "quasi - 
continuum".  This  theory  accounts  for  unimolecular  dipole  dephasing 
and  partial  dipole  recorrelation,  and  it  predicts  the  collapse  and 
revival  phenomena  based  on  the  concept  that  the  level  spacing  of  the 
quasi-continuum  is  close  enough  that  some  continuum-like  behavior  is 
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evident  in  the  transition  dynamics,  but  also  far  enough  apart  so  that 
the  actual  discreteness  of  the  upper  “band"  of  the  final  levels  is  not 
completely  smeared. 

12 

3.  Excitation  of  Intramolecular  Modes 

In  the  previously  discussed  heat-bath  model  of  MPE,  the  active  mode 
is  singled  out  as  the  pumped  system  and  all  other  modes  are  treated  as 
a  heat  bath.  The  coupling  between  the  system  and  the  bath  modes  provides 
the  relaxation  mechanism  where  the  many-body  effects  due  to  the  anharmonic 
mode-mode  coupling  are  phenomenologically  treated  by  the  energy  (y^)  and 
phase  (yg)  relaxation  factors.  This  damping  model,  while  providing  a  simple 
picture  of  the  energy  transfer  processes  and  focusing  on  the  active-mode 
excitation,  is  not  able  to  provide  the  detailed  individual  excitations 
within  the  intramolecular  modes.  To  describe  the  dynamical  excitations 
and  evaluate  the  intramolecular  vibrational  relaxation  (IVR)  rate,  we 
shall  now  present  a  microscopic  model  which  includes  the  anharmonic 
mode-mode  coupling,  in  contrast  to  the  damping  model  where  only  the 
active-mode  excitation  was  included. 

Before  investigating  IVR  for  the  SFg  molecule,  we  shall  develop 
general  formulas  for  describing  the  dynamics  of  intramolecular  mode¬ 
mode  coupling  for  any  polyatomic  system.  The  total  Hamiltonian  for 
an  anharmonic  quantum  oscillator  subject  to  infrared  radiation  may 
be  written  as 

V  *  H0ff  *  XHanh  ■ 


H(t,  Q^,  Q2 »  ••• 


(III. 31) 
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eff 

where  Hq  is  the  effective  unperturbed  Hamiltonian  of  the  system 

(with  N  normal  modes),  xHan^  is  the  intramolecular  anharmonic  coupling, 

and  H'(t)  is  the  laser  field  interaction  Hamiltonian.  In  a  second- 

quantization  representation,  the  individual  terms  of  the  total  Hamiltonian 

12 

may  be  expressed  by 


,eff 

0 


at 

0) .  a  •  a  . 

j  j  j 


iJ*r./2, 
J  J 


ann  £,m,n=l  1  1  J  J  11 


In  Eq.  (11.32)  the  imaginary  part  of  the  complex  effective  Hamiltonian 
ri  takes  account  of  the  level  broadening  of  the  open  system  due  to 
factors  such  as  collisional  (pressure)  broadening  for  a  homogeneous 
system  or  surface-phonon-induced  damping  for  a  heterogeneous  system. 

In  Eq.  (11.33),  two  typical  terms  of  the  anharmonic  coupling  are  shown: 
the  first  term  representing  phase  relaxation  processes,  where  no 
vibrational  energy  among  the  intramolecular  modes  is  exchanged,  and 
the  second  term  representing  the  energy  relaxation  processes  where 
the  photon  energy  is  deposited  in  the  active  mode(s)  and  populated  in 
all  the  other  inactive  mode  via  anharmonic  coupling.  In  Eq.  (11.34), 


(11.32) 


(11.33) 


(III. 34) 


V.(tMj0cos(cot)  is  proportional  to  the  derivative  of  the  dipole  moment 
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and  the  effective  local  electric  field  of  the  j-th  mode  with  laser 

frequency  w.  We  note  that  the  anharmonic  coupling  factors  and 
ilk 

G&nn  are  related  to  the  p-th  derivatives  of  the  anharmonic  potential 
energy  (with  respect  to  the  normal  coordinates)  with  p  =  l  +  m  +  n 
being  the  order  of  the  multiquantum  coupling. 

The  solution  of  the  time-dependent  Schordinger  equation 


ifi[3*/at]  =  Hf, 


(11.35) 


with  the  total  Hamiltonian  given  by  Eq.  (11.31),  may  be  written  as 


14 


y(t)  =  exp 


i  Heff. 
*H0  t  + 


A(t)  +  £(Bj(t)at  -  B|(t)3j) 


*(0), 


(11.36) 


where  A(t)  and  B(t)  are  time-dependent  coefficients  to  be  determined 
by  the  equations  of  motion.  Since  the  energy  population  and  the 
excitations  of  the  system  are  only  characterized  by  the  amplitude 
functions  B - ( t)  and  are  independent  of  A(t),  we  shall  focus  on  the 

J 

equation  of  motion  for  B.(t),  which  can  be  obtained  by  the  substitu- 

J 

tion  of  Eq.  (11.31)  into  Eq.  (11.35)  together  with  the  expression  of 
Eq.  (11.36)  for  the  wave  function: 


(11.37) 
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where  the  external  detuning  (a^)  and  the  internal  detuning  (Aj)  are 

defined  by  a^  =  2ir(vj+^  -  \>j)  -  u  and  Aj  =  2w(£vi  +  mv^  -  nv^).  In 

deriving  Eq.  (11.37),  we  have  used  the  usual  rotating-wave  approximation 

14 

and  operator  algebra  such  as 


[exp(B*a)]F(a+)[exp(-B*a)]  =  F(a++BJ), 


(11.38) 


where  F(a+)  is  any  operator  function  and  B.  are  c-numbers. 

J 

The  important  features  of  the  quantum  equations  of  motion 
[Eq.  (11.37)]  are:  (i)  The  j-th  normal  mode  of  the  system  is  strongly 
coupled  to  the  laser  field  when  Vjq  ?  0  and  the  external  detuning  is 
not  far-off  the  resonance  condition  Aj|  =  0  (due  to  the  anharmonicity, 
we  shall  find  later  that  the  optimal  detuning  is  red  shifted,  i.e., 
A^>0).  (ii)  The  intramolecular  coupling  is  governed  by  the  coupling 
strengths  XF^,  xG^  and  the  internal  detuning  Aj.  For  the  intra¬ 
molecular  coupling  to  be  significant,  we  shall  expect  a  near  resonance 
Aj  :  0,  which  was  also  shown  by  a  classical  treatment-  (iii)  By 
Eq.  (11.37)  and  its  complex  conjugate,  we  find  that  the  equation  of 
motion  for  B*B.  contains  no  T2  dephasing  term,  i.e.,  no  vibrational 
energy  is  exchanged  among  the  intramolecular  modes  due  to  the  T2 
dephasing  processes,  (iv)  The  amplitude  function  is  related  to  the 
average  excitation  in  the  Poisson  distribution 


Pnj(t>  = 


(11.39) 


where  Pnj  is  the  energy  population  of  a  harmonic  oscillator  (for  x  =  0), 

and  n-  and  n",  are  the  vibrational  quantum  number  and  the  average 
J  J 
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excitation  of  the  j-th  mode,  respectively,  where 

nj(t)  =  |Bj(t)|2  =  £  njPnj(t)‘  (II. 40) 

n  j 

_  2 

Therefore  the  average  excitation  n.  =  | B - |  is  a  relevant  quantity  in 

J  J 

describing  the  energy  distribution  among  the  intramolecular  modes,  as 
we  initially  expected. 

We  have  shown  that  IVR  will  be  significant  only  when  the  internal 

resonance  condition  is  met,  i.e.,  the  photon  energy  deposited  in  the 

active  mode  will  be  rapidly  shared  by  those  modes  which  have  frequencies 

close  to  that  of  the  active  mode  with  the  frequency  combination  resonant 

or  near  resonant  with  that  of  the  active  mode.  Furthermore,  from  the 
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concept  of  the  energy-gap  law,  we  know  that  IVR  is  dominated  by  those 
low-order  combinations,  i.e.,  small  values  of  p,  since  it  is  a  strongly 
decreasing  function  of  p  when  p  increases.  We  therefore  may  truncate 
the  order  of  IVR  in  the  sum  of  Eq.  (11.37),  which  otherwise  cannot  be 
numerically  solved  with  a  reasonable  amount  of  computer  time,  particularly 
for  a  large  molecule  system  with  a  large  number  of  vibrational  modes  like 

SF6* 

To  show  the  role  of  the  external  and  internal  detunings  in  IVR,  we 
now  consider  the  SFg  system  in  the  collisionless  low-pressure  regime 
(i.e.,  r-  =  0)  and  focus  on  the  case  of  the  fourth-order  anharmonic 

J 

2 

coupling  ^Q^Qg*  P  3  4,  where  e  is  the  anharmonicity  and  Q3, 

Q3  and  Qg  are  the  normal  coordinates  of  the  (active),  and  vg 
modes.  It  is  known  that  SFg  has  an  0^  symmetric  group,  and  the  frequency 
spectrum  of  the  15  normal  modes  (some  are  degenerate)  is  (in  cm  *): 


=  773.6,  V2  =  642.1,  Vg  =  947.97,  v^  =  615.03,  vg  =  522.9  and 
vg  =  364.  The  corresponding  internal  detuning  then  is  Aj  = 

2Tr(vg-2v2-Vg)  and  the  external  detuning  is  Ag  =  2^  -  u>,  where  w 
is  the  frequency  of  the  tunable  COg  TEA  laser  which  is  near  resonant 
with  the  mode.  The  quantum  coupled  equations  of  motion  [Eq.  (11.37)] 
for  this  case  become 

ifiBj  =  A'cJbjb2  +  X'C^CgBgB*  exp(iAjt)  +  (V1Q/2)2C1  exp(iA^t),  (II. 41a) 


ifiB2  ^A'C^BjJB2  +  2A’C1c|c3B1B*B3exp(-i a* t) 


ifiBg  =  X'CgB^I  +  x'CjC^CgBJBg  exp(-iAjt), 


(II. 41b) 


(II.41c) 


Bj(j*l»2,3)  are  the  related  amplitude  functions  of  the  normal  modes 

J 

vg,  v2  and  Vg  respectively;  A'  =  12A,  aJ  =  2ir(vg  +  Vg  -  2v£)  and  Ag  = 

2irvg  -  oj,  where  vj(j=2,3,6)  are  the  anharmonic-corrected  frequencies, 

i.e.,  vl  =  v.  -  A 1 C i /fi ;  and  C-  are  the  quantization  constants,  C.  = 

J  J  J  J  J 

(fi/2m.a).)1//2.  The  above  system,  for  the  case  of  exact  resonance, 

J  J 
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Ag  =  0,  has  been  discussed  and  the  related  IVR  has  been  estimated. 

Here  we  shall  show  the  dynamical  features  of  IVR  and  find  the  optimal 
external  and  internal  detunings  by  solving  the  coupled  equations  of 
motion  numerically.  We  investigate  the  time  evolution  of  the  average 
excitations,  =  1 6 x  1 2 ,  n"2  =  | Bg | 2  and  rTg  =  |Bg|2,  governed  by  Eq.  (11.41) 
with  the  initial  Boltzmann  populations  at  room  temperature  (T  =  300K) 

B^ (0)  =  0.10,  B2(0)  =  0.214  and  Bg (0)  =  0.417,  and  the  initial  phases 
tan~*[I .(0)/R.(0)]  =  0,  where  I.  and  R.  are  the  imaginary  and  real 

J  J  J  J 

parts  of  B^(j  =  1,2,3). 


Fig.  4  shows  the  peak  values  of  the  average  excitations  of  the  active 
mode  (n£)  and  of  the  v2  mode  (n|)  vs.  the  external  (a£)  and  internal  (a{) 
detunings,  respectively.  It  is  seen  that  the  optimal  detuning  a|  =  17  cm-1, 
and  the  value  of  n^  decreases  for  a|  larger  than  a|.  Similar  behavior  can 
also  be  seen  for  n£.  This  "red-shift"  behavior  of  the  optimal  detuning 
is  a  general  feature  of  any  laser-excited  anharmonic  system  where  the 
energy  spacings  are  closer  for  higher  levels.  Furthermore,  these  non¬ 
linear  effects  due  to  the  anharmonicity  show  the  "bistability"  transi¬ 
tion  in  the  diagram  of  steady-state  average  excitation  vs.  laser 
intensity.  (See  Section  III  for  more  detailed  discussion.) 

Fig.  5  shows  that  the  time  evolution  of  the  excitations  of 
various  modes  for  exact  resonance,  far-off  resonance  and  the  optimal 
detuning.  It  is  seen  that,  by  comparison  of  A  with  B,  the  active-mode 
excitation  (n^  is  coupled  to  the  other  modes  more  significantly  for  the 
resonance  situation  than  for  the  far-off  resonance  situation,  as  we 
predicted  from  the  internal  resonance  criterion.  Note  that  the  peak 

value  of  the  v2  mode  (n|)  decreases  from  0.4  (for  Aj  =  0)  to  0.1  (for 
Aj  =  63  cm-*).  From  these  numerical  results  we  estimate  that  the  rising 

times  of  the  average  excitations  t|  [defined  by  the  first  peak  values  of 
the  average  excitation  "n|(t*)]  range  from  tf  =  27  picoseconds  (ps)  to 
t£  =  70  ps  depending  on  the  internal  detuning  Aj.  However,  the  rising 
times  of  the  v2  mode,  i.e.,  the  IVR  times,  are  on  the  order  of  t£  =  10  ps 
(or  3.3  cm"*).  A  typical  excitation  diagram  is  shown  in  Fig.  5(C)  for 
the  case  of  (a£/2tt,aJ/2tt)  =  (2.71,  1.00). 

4.  Selectivity  in  Multiphoton  Processes 

Selectivity  is  characterized  not  only  by  the  coherent  properties  of 
the  laser  field  but  also  by  the  molecular  properties  of  the  excited  system. 
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Therefore,  the  types  of  selectivity  in  infrared  MPE  can  be  classified 

according  to  the  relation  between  the  various  relaxation  times  of  the 

excited  system  and  the  energy  pumping  rate  of  the  laser  field.  Let  us 

first  define  R.  and  R.  as  the  intra-  and  intermolecular  vibra- 
mtra  inter 

tional  energy  transfer  rates,  respectively,  RyT  as  the  relaxation  rate 
for  molecular  vibrational-translation  coupling,  i.e.,  Ryy'1  is  the 
time  for  complete  thermal  equilibrium  to  be  reached  in  the  molecular 
mixture,  and  Wgxc  as  the  rate  of  vibrational  multiphoton  excitation  of 
the  molecule.  We  can  distinguish  four  different  types  of  selectivity 
depending  on  the  relative  magnitudes  of  the  relaxation  rates  and  the 
laser  excitation  rate  (note  that  the  relaxation  rates  Rjnter  and 
RVj  are  pressure  dependent,  although  for  low  pressures  we  shall 
expect  that  R„  «  «,„ter  «  IWa1' 

(I)  Mode  (bond) -selective  excitation  (Wexc  »  Rintra)-  A 
certain  mode  or  a  functional  group  of  a  polyatomic  mol  eel ue  is  in  a 
nonequilibrium  state  which  has  a  higher  vibrational  temperature  Teff 
[defined  by  Eq.  (11.24)]  as  compared  with  the  remaining  modes  or  func¬ 
tional  groups.  This  is  the  situation  of  a  long  lifetime  or  high  pumping 
rate. 

(II)  Molecular-selective  excitations  (Rintra  »  Wgxc  »  Rinter)- 
In  this  case, the  absorbed  photon  energy  is  rapidly  randomized  within 
the  excited  molecule  in  which  the  local  vibrational  temperature  is 
higher  than  the  overall  translational  temperature  of  the  mixture  of 
different  molecules. 

(III)  Vibrational -selective  excitation  (Rinter  >>  Wgxc  >>  Rvj). 

In  this  more  moderate  condition,  vibrational  equilibrium  among  all  the 
mixed  molecules  is  reached,  but  there  is  still  no  overall  thermal 
relaxation.  This  situation  prevails  for  low  pressures,  where 


vibration-translation  relaxation  rates  are  lower  than  vibration- 
vibration  relaxation  rates. 

(IV)  Nonselective  thermal  excitation  {R^  »  Wgxc).  This  is  the  situation 
of  thermal  excitation  of  all  the  molecules  in  a  mixture  by  a  low-power  CW 
laser.  This  area  of  infrared  thermal  chemistry  is  of  interest  for  a  hetero¬ 
geneous  system,  e.g.,  species  adsorbed  on  a  solid  surface,  where  the  laser 
radiation  is  used  to  excite  the  adsorbed  gas  molecules  without  significantly 
heating  up  the  solid  surface  (assuming  the  phonon  coupling  to  be  small). 

We  shall  discuss  this  type  of  excitation  in  Section  III. 

The  attraction  of  infrared  laser  chemistry  is  that  if  the  photon  energy 
can  be  deposited  and  maintained  in  a  specific  vibrational  mode  (or  functional 
group),  a  selective  reaction  channel  may  be  induced  by  that  mode  (functional 
group).  Therefore,  one  of  the  critical  questions  concerning  MPD  processes 
is  whether  the  photon  energy  remains  localized  in  the  pumped  mode  (or 
molecule)  long  enough  to  result  in  a  mode-selective  (or  molecule-selective) 
reaction.  Most  theoretical  models  of  MPD  assume  that  IVR  is  very  fast  (on 
the  order  of  a  picosecond),  and  statistical  approaches  such  as  RRKM  theory 
have  been  applied  successfully  to  a  number  of  experimental  results.*^’2® 
Recently,  however,  several  experiments  examining  product  branching  ratios 

in  relatively  complex  molecules,  e.g.,  cyclopropane,  have  shown  that  the 
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laser  selectivity  cannot  be  explained  by  the  statistical  theory.  The 

evidence  suggests  that  the  laser-induced  reactions  could  result  from  a 

nonergodic  or  partially  mode-selective  excitation.  Several  research 

groups  have  also  examined  the  bond  localization  character  of  large 

39-41  42 

molecules,  e.g.,  benzene,  both  experimentally  and  theoretically. 

This  nonstatistical  behavior  is  explained  by  a  local-mode  model  in  which 
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the  total  system  is  divided  into  two  groups  of  vibrational  modes.  Within 
each  group,  the  vibrational  modes  are  strongly  coupled  and  the  photon 
energy  quickly  randomized.  Between  groups,  however,  the  coupling  is 
considerably  weaker,  and  intergroup  randomization  rates  are  therefore 
appreciably  slower  by  the  concept  of  the  energy-gap  law. 

Another  theory  called  the  "restricted  IVR  model"  has  been  recently 
proposed  for  the  possibility  of  laser  selective  photochemistry.  This 
model  provides  a  basis  for  a  quantitative  description  of  MPD  when  IVR  pro¬ 
cesses  occur  on  a  time  scale  competitive  with  the  unimolecular  reaction 


and  radiative  pumping  rate.  Note  that  the  statistical  theories  of  uni¬ 
molecular  decay,  like  the  RRKM  theory,  preclude  any  possibility  of 
selectivity  because  of  their  assumption  of  very  fast  IVR  rates.  For 
MPD  of  the  SFg  molecule,  the  IVR  rate  is  estimated  on  the  order  of 
10  sec  ,  so  that  selective  bond  breaking  (or  mode  selection)  is 
possible  only  when  one  uses  laser  radiation  with  very  short  pulse  dura¬ 
tion  and  about  two  orders  of  magnitude  more  intense  than  currently 
20 

available  lasers.  However,  for  laser-induced  decomposition  of  cyclo¬ 
propane,  the  experimental  results  strongly  indicates  a  selective 
mechanism  and  may  be  analyzed  by  the  restricted  IVR  model  in  which 
a  slower  IVR  rate  (10*°  sec”*)  is  estimated.  According  to  this  model, 
the  vibrational  modes  of  the  system  are  divided  into  two  groups  with 


an  intragroup  relaxation  rate  R^ntra  and  an  intergroup  rate  R^nter  = 

8  R.  *  .  The  theory  then  allows  for  a  smooth  transition  from  the  RRKM 

mtra 

theory  (with  g  -*•  1)  to  the  Slater  theory  (with  g  >>  1),  and  group- 


selective  processes  are  possible  when  g  <<  1. 
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5.  Multiple-Frequency  Laser  Excitation 

Recent  experiments  using  two  infrared  lasers  simultaneously  have 

demonstrated  that  MPD  can  be  produced  at  significantly  reduced  intensities 

compared  to  single-laser  excitation,  and  consequently  a  sharper  isotopic 

44 

selectivity  can  be  achieved.  In  these  experiments,  a  lower  inten¬ 
sity  near-resonant  radiation  and  a  higher-intensity  off-resonant  radiation, 
neither  of  which  produces  dissociation  by  itself,  together  give  dissociation. 

Both  classical  and  quantum  models  have  been  proposed  to  analyze  MPD  enhanced 

45-48 

by  two  or  many  infrared  lasers.  The  effect  of  laser  polarization  was 

also  investigated,  and  it  was  found  that  the  excitation  is  more  effective 

when  the  polarizations  of  the  two  lasers  are  parallel  than  when  they  are 
49 

perpendicular.  More  recently,  the  interaction  of  a  diatomic  molecule 

50 

with  two  near-resonant  lasers  was  studied  by  a  classical  anharmonic  model. 

Within  a  heat-bath  treatment,  the  total  Hamiltonian  of  an  anharmonic 
quantum  oscillator  subject  to  two  infrared  lasers  may  be  expressed  in  the 
same  form  as  that  of  a  single-laser  excitation,  Eq.  (11.12),  except  that 
the  interaction  term  H^p  is  now  extended  to  the  form 


HAF  +  Vt)]  (a++a)’ 


(11.42) 


where  £(*)  are  the  pumping  rates  of  the  laser  radiations  with  fre¬ 
quencies  w1  By  the  same  procedures  as  previously  introduced  for 
Eq.  (11.15),  we  obtain  the  rate  equation  for  the  average  excitation 
via  two-laser  processes  as 


d<n>  =  qlIl 
dt  fiw  ^ 


q2*2 


F<tlI2> 


1/2 


Y1(<n>  -  N). 


+ 


(11.43) 
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1/2 

F( 1 1 12 )  /  1nter^erence  term,  with  an  appropriate  function  F 

related  to  the  oscillatory  functions  sin(A^t)  and  cos(A^t),  where 
A.j  =  (i  =  1,2),  and  the  absorption  cross  sections  are 

given  by,  similarly  to  Eq.  (11.16),  as 

Ai  (VY2)f,a)i 

1  (Ai-2e*<n>)2+(y1+Y2)2/4 

For  lack  of  temporal  mutual  coherence  of  the  two  laser  sources,  we  may 
neglect  the  interference  term  and  focus  on  the  cooperative  effects  of 
these  two  different  lasers  on  the  average  excitation  of  the  system. 

As  shown  in  Fig.  6,  we  see  that  at  the  optimal  detuning  with  (a^,  a2)  = 
(2,  8)  in  the  unit  of  ( y1+y2 ) »  the  system  has  a  maximum  excitation 
(curve  A).47  It  is  also  interesting  to  note  that  a  single-laser 
excitation,  where  the  pumping  rate  equals  the  sum  of  the  individual 
rates,  is  less  effective  than  the  two-laser  simultaneous  excita¬ 
tions.  These  cooperative  effects  of  MPE  by  two  different  lasers  are 
characterized  by  the  nonlinearity  of  the  absorption  cross  section 
through  the  anharmonic  potential  of  the  system. 

The  above  results  based  on  a  quantum  heat-bath  model  may  also  be 
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obtained  by  a  classical  approach.  The  equation  of  motion  is  given  by 
an  extension  of  Eq.  (II. 3), 

<Q>  +  y<Q>  +  co2ff<Q>  =  £  ^[eEi  cos(u)it)]/m. 


(11.44) 


(11.45) 


and  the  corresponding  power  absorption  is  given  by  the  extended 
expression  of  Eq.  (II. 5) 
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Y 

4m 

(a.-K*A2)2+y2 

+  Pj(t), 


(11.46) 


where  the  interference  term  is  given  by 


Pj(t) 


(•!, 


E1E2 


_  o  _  _ 

y  cos[Urw2)t]  -  (A.j-K*A^)  sinC^-w^tj 


(A, 


•K*A2)2+y2 


(11.47) 


To  compare  the  above  classical  results  with  those  of  the  quantum  models, 
we  note  that  the  interference  term  F(I1I2)1//2  of  Eq.  (11.43)  is  pro¬ 
portional  to  EjE2  in  Eq.  (11.47),  since  the  intensity  is  related  to 

2 

the  electric  field  by  I..  =  E.j/(8ir/c) .  Furthermore,  the 

classical  and  quantum  correspondence  may  be  easily  obtained  by  the 
following  quasi -quantum  relations: 

(II. 48a) 
(II. 48b) 

where  e  is  a  classical  charge,  u'(0)  is  the  matrix  eloment  of  the 
derivative  of  the  dipole  operator,  eval  /ated  at  the  equilibrium  point,  and 
<n>tiu)Q  is  the  photon  energy  absorbed  by  the  quantum  osci'lator  and  hence 
equals  the  energy  of  the  classical  osicllator  with  amplitude  A. 

The  enhancement  of  MPE  via  multiple-frequency  lasers  has  been 
investigated. it  was  suggested  that  if  a  single  laser  of  intensity 
I  produces  a  maximum  excitation  of  frequency  w,  then  the  same  maximum 


e  =  u'(0). 


1  2  2 

•^mugA  *  <n>fiojg, 
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excitation  may  be  produced  by  N  lasers  of  intensity  I/N^  and  frequencies 

=  u>Q  +  k((jj-u)g)/N,  with  k  =  1,  2,  ...N.  The  enhancement  of  MPE  or  the 

reduced  power  law  of  the  laser  intensity,  I/N  ,  compared  with  that  of 

the  separate  N-laser  excitation,  I/N,  in  principle  can  be  achieved  if 

one  can  modulate  the  phases  of  the  laser  pulses  such  that  the  Rabi 

frequencies  a.e  synchronized  with  the  excitations  between  pairs  of 

the  anharmonic  energy  levels,  i.e.,  tt  pulses  are  required  in  addi- 

45 

tion  to  the  resonance  condition  . 
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the  tantalising  prospect  of  selectively  altering  the  dynamics  of  chemical 
reactions,  it  remains  true  that  the  literature  (both  experimental  and 
theoretical)  is  generally  characterised  by  the  assumption  of  very  intense 
laser  fields,  and  these  are  somewhat  impractical  for  two  reasons.  First, 
intense  fields  have  a  habit  of  producing  stray  side  effects  that  may 
interfere  with  the  process  of  interest  (gas  breakdown  due  to  ionisation 
is  a  typical  example).  Second,  and  perhaps  more  important,  is  that 
intense  fields  are  consistent  only  with  pulsed  lasers,  and  these  are 
unfortunately  "switched  off"  most  of  the  time  they  are  operating:  the 
pulse  length  is  typically  extremely  short  compared  to  the  time  between 
pulses.  Thus,  any  cross  section  which  one  calculates  is  effectively 
"reduced"  by  many  orders  of  magnitude  in  a  realistic  experiment.  In 
our  opinion,  this  is  one  of  the  most  serious  limitations  posed  to 
conventional  laser-induced  chemistry. 

Of  course,  one  could  use  the  radiation  in  a  theoretically  trivial 
way,  to  excite  the  reactants  into  quantum  states  that  make  the  desired 
reaction  more  feasible.  However,  in  many  cases  this  may  not  be  desirable. 
For  example,  it  may  not  be  possible  to  excite  the  reactants  because  the 
available  lasers  are  not  resonant  with  the  transition.  Moreover,  even 


after  the  reactants  are  excited,  one  still  has  only  a  relatively  crude 
control  over  the  outcome  of  the  reaction:  undesirable  side  reactions 


may  also  take  place,  or  the  quantum  efficiency  for  the  reaction  of 
interest  may  turn  out  to  be  low. 
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One  can  to  a  large  extent  get  around  these  problems  by  taking 
advantage  of  the  huge  numbers  of  molecular  states  that  are  formed 
during  the  collision.  This  means  that  a  far  greater  range  of  exper¬ 
imental  conditions,  such  as  collision  energy  and  laser  frequency,  can 
in  principle  give  rise  to  a  reaction.  However,  one  again  runs  into 
the  problem  that,  since  many  collisions  are  short-lived,  any  photons 
that  are  to  be  absorbed  during  the  collision  have  to  be  absorbed  in  a 
relatively  short  period  of  time.  This  of  course  means  that  in  order  to 
observe  any  reaction  at  all,  one  must  provide  more  photons,  and  this 
raises  the  intensity  requirements  of  the  laser  to  levels  that  may  be 
impractical  for  the  reasons  we  have  given.  [The  kind  of  intensity  that 

has  been  discussed  in  the  literature  is  typically  greater  than,  and 

o 

usually  much  greater  than,  1  MW/ cm' .] 

There  are  two  related  factors  by  which  radiative  collision  cross 
sections  can  be  increased.  The  first  is  resonance  formation,  and  the 
second  is  low-energy  collisions.  They  are  related  since  the  effect 
in  both  cases  is  to  increase  the  time  of  collision,  thus  allowing  the 
colliding  partners  more  chance  of  absorbing  a  photon.  It  is  perhaps 
not  surprising  that  resonance  formation  leads  to  an  enhancement  of  the 
cross  section,  but  it  may  seem  rather  paradoxical  that  a  decrease  in 
collision  energy  should  have  a  similar  effect.  After  all,  the  Arrhenius 
Law  implies  that  chemical  reactions  proceed  at  a  slower  rate  as  the 
temperature  is  lowered.  However,  it  must  be  remembered  that  in  a 
chemical  reaction  there  is  always  an  activation  energy  barrier,  and 
the  reason  that  reactions  become  faster  at  higher  energies  is  simply 
that  more  reactants  have  sufficient  energy  to  overcome  this  barrier. 
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But  if  the  activation  energy  is  supplied  by  a  photon,  there  is  then  no 
need  for  a  high  collision  energy--indeed  low  energies  are  favored  since 
they  entail  longer  collision  times.  At  exceptionally  low  collision 
temperatures  (less  than  a  few  degrees  Kelvin),  the  collision  must  be 
treated  quantum  mechanically  since  the  scattering  wavelengths  become 
comparable  to  the  molecular  dimensions,  except  for  very  heavy  atomic 
and  molecular  reactants.  We  shall  first  consider  an  atom-atom  collision 
in  the  presence  of  laser  radiation  and  investigate  the  ramifications 
of  radiative  resonance  formation  when  more  than  one  laser  is  employed. 

Later,  we  shall  sketch  some  work  which  is  underway  to  investigate  the 
feasibility  of  driving  certain  chemical  reactions  with  lasers  of  low 
intensity. 

We  consider  an  atom-atom  collision  in  the  presence  of  two  lasers, 

LI  and  L2.  The  total  Hamiltonian  is  given  in  atomic  units  as 

H  =  -  -^  V2  +  He1  +  HS0  +  Hrad  +  Hjnt  +  Hjnt  ,  (11.49) 

pi 

where  y  is  the  reduced  mass  of  the  system,  and  H  is  the  electronic 
Hamiltonian  in  the  absence  of  spin-orbit  coupling,  with  eigenvalues 
(in  the  Born-Oppenheimer  approximation)  which  are  the  molecular 
potentials.  Hs0  is  the  spin-orbit  Hamiltonian.  Hrad  is  the  Hamiltonian 
of  the  radiation  fields. 
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in  spherical  coordinates:  o  =  0,  +1  and  -1  denotes,  respectively, 
light  plane  polarized  in  the  z-direction,  and  right  and  left 

circularly  polarized  and  propagating  along  z.  The  eigenvectors 

_ 

of  H  are  the  Fock  states  In.  >,  such  that  a.  |n.  >  =  /n .  In.  -1> 

1  la  la1  la  la1  la 

and  a.  In.  >  =  /n.  In.  +1>.  The  radiative  interaction  Hamiltonians 
la1  la  la1  la 

within  the  dipole  approximation  are 


H*nt 


r  f  *-i 

e.l[a.  y  +  a.  y  ] 
l  la  a  la  aJ 
a 


(11.51) 


where  yQ  is  the  space-fixed  component  of  the  dipole  operator  and 
2™.  1/2 

ei  =  Hj/'T')  *  with  Vi  as  the  quantization  volume  of  the  field  i. 

Note  that  yQ  can  be  related  to  the  molecule -fixed  components  of  the 
dipole  operator,  y  ,  by  the  transformation 


(II. 52) 


l  ^ 

where  D^R)  is  the  Wigner  rotation  matrix. 

We  are  particularly  interested  in  the  two  curve-switching  processes 

which  are  shown  in  Fig.  7.  Process  I  has  been  called  the  "photon-catalytic 
51 

effect".  This  is  because  there  is  no  net  loss  of  photons  from  the 
field,  since  for  each  absorption  there  is  an  emission  of  the  equivalent 
photon  (same  frequency,  polarisation  and  direction).  Such  a  process 
would  therefore  appear  to  hold  some  promise  for  laser-induced  chemistry. 
However,  from  a  practical  point  of  view,  laser  cavity  losses  are  far 
more  important  than  losses  to  and  from  a  molecular  system.  Therefore, 
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the  efficiency  of  the  process  is  hardly  affected  by  the  few  photons  that 

are  returned  to  the  radiation  field  by  the  molecular  systems.  Moreover, 

the  overall  cross  sections  are  not  in  general  very  high.  Process  II  on 

the  other  hand,  while  apparently  a  more  contrived  means  of  achieving  the 

same  result,  turns  out  to  have  two  advantages.  First,  it  allows  the 

colliding  species  to  approach  at  exceptionally  low  energies,  since  the 

energy  defect  for  curve  switching  is  made  up  by  absorption  of  the  second 

photon  from  L2.  Second,  it  provides  greater  scope  for  optimization  of 

the  collision  conditions:  by  changing  the  bound-to-bound  transition 

it  is  possible  to  optimize  the  bound-continuum  interactions  so  that 
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the  overall  cross  section  can  be  maximized. 

In  order  to  treat  resonance  scattering,  we  consider  a  Breit-Wigner 
separation  of  the  T-matrix.  The  radiative  interaction  gives  rise  to 
resonances  in  the  scattering  spectrum,  while  the  field-free  interaction 
may  be  considered  to  be  a  slowly-varying  function  of  the  scattering 
energy.  It  is  therefore  convenient  to  separate  the  T-matrix  into  slowly 
and  rapidly  varying  parts: 

T  =  TP  +  TR  , 

P  R 

where  T  is  the  slowly  varying,  potential  scattering  T-matrix  and  T 

represents  the  resonance  scattering.  Following  the  formulation  of 
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Feshbach,  we  define  the  following  projection  operators: 

P  =  I  /  |  ip  ,><4,  E ,  | dE '  (11.53) 

a 

Q  =  I|4>n><4>nl  • 
n 


(11.54) 
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Here,  a  is  a  complete  set  of  angular  momentum  quantum  numbers,  and  E'  is 
the  energy  of  the  continuum  state  specified  by  a.  The  n  index  is  a  set 
of  vibrational  and  angular  momentum  quantum  numbers  which  specify  a  bound 
state  of  the  system.  The  space  of  Q  can  be  further  subdivided  into  R  and 
S,  such  that  R  +  S  =  Q: 

R  =  Eiv<4J  (II. 55) 

m 


S  =  E  K'^n'  I  »  (11.56) 

n,  n  n 
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It  is  now  possible  to  write  the  T-matrix  for  resonance  scattering  as 

Tfi  =  ^f  lHPRS7HRpl^i>  (IL57) 

where  fl  =  (E-HRR-HRp  gJ  HpR)~l  (11.58) 


Hxy  =  XHY  . 


(11.59) 


Gp  is  the  Green's  function  for  potential  scattering,  i.e.,  with  the 

radiative  interaction  switched  off  and  the  scattering  determined  by 
e  1 

the  term  H  in  the  Hamiltonian.  This  term  has  off-diagonal  elements 
in  a  diabatic  basis  which  are  indirectly  due  to  spin-orbit  coupling, 
reflecting  a  curve-switching  even  in  the  absence  of  radiation.  The 
scattering  eigenfunctions  and  ify  are,  respectively,  the  out-wave 
for  this  potential  scattering  problem,  starting  in  the  initial  state  i, 
and  the  in-wave  starting  in  the  final  state  f.  These  are  total  potential 
scattering  wavefunctions  consisting  of  a  plane  wave  in  the  specified 
channel  and  spherical  waves  in  all  channels.  These  may  be  written 
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in  terms  of  the  standing-wave  solutions  of  the  potential  scattering, 

l*aE5  by  using  the  orthogonal  matrix  W  which  diagonal ises  the  potential 
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scattering  S-matrix.  The  result  is 


jR  =  y 

fi 


2  W,  H  Q,  H_  1 0  W 


mm'  a8 


fa  am  “mm1  nm‘8  "6i  ’ 


(II. 60a) 


i(DWBA) 


Hfm  nmrn'  Hm'i  * 


(II. 60b) 


where  m  and  m'  are  bound  intermediate  states  of  the  collision  complex,  and 


(11.61) 


where  is  a  standing-wave  eigenfunction  of  H.  In  the  distorted-wave 
Born  approximation  ( DWBA ) ,  W  is  the  unit  matrix,  and  furthermore, 

D 

Tfi-  =  H^.  .  Therefore,  we  have  for  the  total  T-matrix 

Tfi  <»*>  =  Hf1  *  E  .  "fmW  Vi  (n'62) 

mm 


The  quantity  of  central  interest  is  Q  since  it  contains  the  detailed 
information  about  radiative  couplings.  Possible  kinds  of  laser-induced 
resonance  scattering  illustrated  in  Fig.  8.  By  defining  we  also 
define  a  particular  pathway:  for  example,  in  Eq.  (II. 60b)  ^gs^stHta 
the  amplitude  for  Process  2  in  Fig.  8.  Notice,  however,  that  all  the 
ancillary  absorptions  and  emissions  induced  by  L2  are  implicitly  taken 
into  account  in  the  inversion  of  that  is,  depends  not  only  on 
|s>  and  |t>  but  also  on  |u>.  We  first  evaluate  Q-*  by  expanding  G*  in 
terms  of  its  open  channel  eigenfunctions: 
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where  V  is  the  nonradiative  interaction  operator  (e.g.,  spin-orbit 
coupling),  and  Tse  is  the  width  due  to  spontaneous  emission.  The 
3x3  matrix  ff*  is  now  inverted  to  give  the  following  result: 
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where  i.  *  E  •  E|  ■  F^,  and  E ^  is  the  unperturbed  energy  of  the 
bound  state  i . 

Using  Eqs.  (11.60)  and  (11.67),  we  can  now  write  down  the  T-matrix 
for  Processes  I  and  II  in  Fig.  7.  We  assume  no  predissociation  channel 
and  neglect  spontaneous  emission: 


(11.68) 


(11.69 


It  is  noteworthy  that  as  the  energy  defects  A  of  the  dressed  states 

approach  zero,  the  processes  have  the  same  order  of  magnitude.  Indeed, 

the  actual  magnitude  of  Process  II  might  even  greatly  exceed  that  of 

Process  I,  even  though  Process  II  is  a  three-photon  process.  We  have 

shown  by  means  of  a  detailed  calculation  on  the  XA  +  F  system  that  this 
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is  indeed  the  case,  provided  the  initial  kinetic  energy  is  low. 

Moreover,  the  intensity  requirements  on  the  radiation  are  also  very  low 

2 

LI  operating  at  3  kw/cm  can  give  rise  to  gas  kinetic  cross  sections  for 
Process  II.  On  the  other  hand,  calculations  have  shown  Process  I  to 
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give  rise  to  a  relatively  small  cross  section  (around  0.05  A  )  even 

2 

though  the  intensity  of  LI  is  large  (0.3  MW/cm  ).  We  emphasize  that 
large  effects  can  be  obtained  from  low-intensity  lasers,  provided 
resonances  are  involved  and  the  initial  collision  energy  is  extremely 
low. 

Recently,  there  have  been  a  number  of  interesting  experiments 

which  appear  to  have  demonstrated  chemical  reactions  that  are  induced 

by  lasers  resonant  with  any  asymptotic  states.8^"88  Although  this  is 

an  encouraging  start,  many  of  these  reactions  seem  to  have  a  very  small 

total  cross  section.  It  would  be  useful  to  have  a  theory  of  laser-driven 

chemical  reactions  which  would  enable  us  to  decide  what  kind  of  reactions 

may  be  feasible,  and  to  perform  some  quantitative  estimates  of  cross 

sections.  To  this  end  we  may  consider  a  model  of  chemical  reaction  which 

is  based  on  the  theory  of  the  atom-atom  collision  discussed  above  and 

which  conforms  to  the  criteria  we  have  established.  An  important  condition 

is  the  ability  to  form  resonances  during  the  collision.  We  have  so  far 

considered  radiative  resonance  formation,  but  it  may  be  useful  to  begin 

by  considering  systems  which  form  resonances  in  the  absence  of  radiation. 

An  important  class  of  such  systems,  which  are  readily  obtained  experimentally, 

are  van  der  Waals  molecules  (VDWM)  formed  in  their  ground  electronic  state. 

We  may,  for  example,  be  interested  in  the  reaction 

He  +  Ig  +  fiw  -  ■>  Hel*  +  I*  ,  (II  7C 

where  the  asterisk  signifies  an  excited  electronic  state. 

Consider  a  coll  inear  reaction,  A  +  BC  +  flw  — *  BC*  +  C*  , 
where  the  bonds  BC  and  AB*  are  strong  and  the  bonds  AB  and  BC*  are 
relatively  weak  (see  Fig.  9).  This  kind  of  condition  would  be  expected 
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in  a  typical  VDWM.  In  the  absence  of  radiation,  resonances  occur  in 

each  individual  electronic  state.  For  example,  the  BC  bond  may  gain 

a  quantum  of  vibrational  energy  through  coupling  to  the  AB  bond, 

which  then  enters  a  bound  state  of  motion.  The  atoms  A,  B  and  C 

then  co-exist  in  a  quasibound  state  which  is,  however,  unstable 

with  respect  to  the  reverse  process,  predissociation  of  atom  A. 
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This  effect  is  known  as  "vibrational  predissociation". 

Let  us  outline  the  construction  of  a  scattering  theory  in  which 

the  predissociation  takes  the  form  of  a  Feshbach  resonance.  The 

radiative  interaction  between  electronic  states  could  be  handled 

in  much  the  same  way  as  for  the  atom-atom  problem.  The  two  electronic 

states  have  the  coordinate  representation  Xj([p];  R,r)  and  x2([p] ;  R,r), 

respectively,  where  [p]  are  the  electronic  coordinates  and  r  and  R  are, 

respectively,  the  length  of  the  strong  bond  [I-I  in  the  above  reaction] 

and  the  distance  from  the  more  weakly-bound  atom  (He)  to  the  center 

of  mass  of  the  two  strongly-bound  atoms  (see  Fig.  9).  Averaging  the 

total  Hamiltonian  H  over  the  electronic  coordinates,  we  obtain  the 

effective  (Born-Oppenheimer)  Hamiltonian  for  the  nuclear  degrees  of 

freedom: 

<x1|H|x2>  =  V12(R,r) 

<X1  lHlxi>  =  Tab  +  TCjBA  +  mx(R)  +  p2(r)  +  V^R.r) 

<x2 1 H I x2>  ■  Ta>bc  +  TgC+  m2(R)  +  p2(r)  +  V2(R,r) 


(II. 71a) 
(II. 71b) 

(II.71c) 


Here,  T  stands  for  the  kinetic  energy  operator  of  the  subscripted  species, 
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and  m  and  p  are  Morse  functions  representing  the  potential  energy  associated 
with  the  two  "bonds"  (strictly  speaking,  the  weak  bond  is  not  a  conventional 
bond  since  its  length  is  measured  to  the  center  of  mass  of  the  other  bond). 

and  V ^  represent  the  three-body  portions  of  the  ground  and  excited  potential 
surfaces,  which  are  coupled  only  through  the  radiative  coupling  V^. 

We  wish  to  write  the  T-matrix  for  the  chemical  reaction  in  terms  of  the 
operators  P  and  Q  which  project,  respectively,  onto  the  open  channels  and 
the  eigenvectors  of  a  zeroth-order  Hamiltonian  defined  to  be  the  effective 
Hamiltonian  of  Eq. (11.71)  minus  V-j  and  Vg.  These  are  then  given  by 

P  =  Z  /dE1  In^Jx^Ejl  +  I  JdE2|n2E2><n2E"|  (II. 72a) 

nl  v2 

Q  =  £  | ni v, xn^Vp |  +  £  ln2V2><n2V2 I  ’  (II. 72b) 

nlvl  n2V2 

where  ) n jEj>  «  lnj>]Ej>»  etc.,  j  being  the  electronic  state  index.  |E^>  is 
the  state  vector  for  inelastic  scattering  in  the  ground  electronic  state,  con¬ 
sisting  of  a  plane  wave  in  the  channel  labeled  by  the  asymptotic  energy  E,  plus 
outgoing  spherical  waves  in  all  other  channels.  Similarly,  [E^  is  an  incoming 
wave  for  the  rearrangement  channel.  The  quantum  numbers  n-  and  v •  label, 

J  J 

respectively,  the  quantum  states  of  the  Morse  oscillators  m.  and  p..  We  consider 

J  J 

in  addition  the  projectors  R  and  S  such  that  R  +  S  =  Q.  The  space  of  R  is 

defined  by  those  bound-state  vectors  which  give  rise  to  strong  resonances  (these 

66 

can  be  called  "doorway  states"  ),  where  S  is  comprised  of  all  other  bound- 
state  vectors.  It  is  now  possible  to  write  the  general  T-matrix  element  connecting 
the  initial  and  final  states: 


(11.73) 


Equation  (11.73)  contains  all  the  information  relevant  to  a  chemical 
reaction,  including  the  effect  of  the  radiation  field.  It  may  at  first  sight 
appear  analogous  to  the  DWBA  picture  of  a  chemical  reaction,  but  it  is  not. 

It  does  in  fact  represent  the  exact  solution  to  the  Lippmann-Schwinger  equa¬ 
tion,  with  the  complication  that  the  effective  (optical)  potential,  HpRaHRp, 
is  not  only  complex  but  also  nonlocal  in  the  coordinates  R  and  r.  The  problem 
hinges  on  the  representation  of  a.  This  is  given  in  terms  of  the  width  matrix 
r,  the  level -shift  matrix  F  and  the  matrix  V  whose  off-diagonal  elements  describe 
radiative  bound-bound  transitions.  The  off-diagonal  elements  of  r  consist  of 
the  (nuclear)  bound-continuum  interaction  (both  radiative  and  nonradiati ve) . 

In  other  words,  r  couples  the  open  and  closed  channels.  Under  certain  condi¬ 
tions  we  expect  these  matrices  to  be  representations  of  local  operators,  where 
these  conditions  might  hinge  on  the  validity  of  the  Markov  approximation  used 
in  time-dependent  many-body  theory.  Here  a  nonlocal  kernel  K(t.t')  in  the  rate 
equation  for  a  subsystem  is  replaced  by  a  local  function  of  time: 

K(t.t')  ;  K(t)6(t-t' ).  This  means  physically  that  the  motion  of  the  subsystem 
at  one  time  is  completely  uncorrelated  with  its  motion  at  some  other  time, 
i.e.,  there  is  no  memory  in  the  motion  (due  to  randomization  on  a  very  short 
time  scale).  Similar  situations  are  expected  to  pertain  for  long-lived  resonances 
in  three-body  collisions,  which  is  quite  different  from  what  would  be  expected 
in,  say,  a  direct  collision  or  stripping  reaction  where  the  initial  and  final 
states  are  much  more  correlated  and  the  nonlocality  of  the  potential  is  more 

important.*^ 

We  could  initially  assume  a  local  optical  potential.  Then  the  matrix 
representation  of  u  in  the  space  projected  by  R  will  consist  of  radiative  bound- 
bound  transitions  between  reactant  and  product  states;  it  will  also  consist  of 
nonradiative  predissociation  widths  which  couple  the  bound  and  scattering  states 
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for  either  the  product  or  the  reactant  channels  taken  separately.  Roughly 
speaking,  the  probability  of  reaction  may  be  regarded  as  the  probability  uf 
association  multiplied  by  the  probability  of  dissociation  into  a  particular 
product  channel.  [This  picture  should  not  be  carried  too  far,  however,  be¬ 
cause  it  is  essentially  perturbative  and  we  are  not  using  perturbation  theory.] 

The  rather  complicated  effect  of  nonlocality  will  eventually  have  to  be  addressed, 
especially  since  there  is  already  evidence  that  nonlocality  becomes  more  important 
at  low  energies and  this  is  just  the  scattering  regime  we  are  interested 
in.  A  simplifying  feature  is  that  nonlocality  appears  to  be  separable  at  low 
energies,  i.e.,  fE(R,R' )^gE(R)hE(R' ) 

Finally,  we  would  like  to  say  a  brief  word  about  transition-state  spectroscopy 

in  the  light  of  what  has  been  discussed  in  this  section.  There  has  recently  been 
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interest  in,  and  experimental  evidence  for,  such  spectroscopy.  ’  We  wish  to 
emphasize  that  the  transition  state  may  be  regarded  as  a  set  of  resonance  states 
formed  during  a  collision.  Spectroscopy  of  the  transition  state  then  involves 
radiative  transitions  between  such  states.  Spectroscopy  and  energy  transfer 
are  distinguishable  only  according  to  what  one  wishes  to  detect:  if  one  chooses 
to  detect  emitted  (or  absorbed)  photons,  then  it  is  spectroscopy,  and  if  particles 
are  detected  then  it  is  energy  transfer  (scattering).  The  spectroscopic  "probe" 
may  thus  be  regarded  merely  as  one  contributor  to  the  total  S-matrix  for  particle 
scattering.  Thus,  the  methodology  for  treating  spectroscopy  of  the  transition 
state  would  be  essentially  what  we  have  outlined  for  radiative  collisions.  We 
note  that  L2  acts  not  only  as  a  spectroscopic  probe  but  also  as  the  agent  of 
energy  transfer  within  the  manifold  of  transition  states.  Thus,  when  L2  has 
a  significant  effect  on  the  overall  scattering  event,  the  distinction  between 
these  two  functions  is  unclear.  They  are,  however,  to  be  thought  of  in  a  unified 
way,  as  we  have  mentioned. 


An  intriguing  potential  application  of  radiative  resonance  scattering 
is  to  isotope  separation.  Consider  the  radiative  velocity-changing  collision 
depicted  in  Fig.  10.  We  may  represent  the  event  depicted  in  the  top  part  as 

-fiu)-. 

A  +  B(Ei)  - U  AB(vta)  - AB(v'a')  - -+  A  +  B(Ef )  (11.75) 

where  +fiu)  denotes  absorbtion  and  -hto  denotes  emission.  The  vibration-rotation 

level  of  the  AB  complex  is  denoted  by  (v,Jt).  The  selectivity  of  such  a  process 

arises  from  the  different  vibration-rotation  spacings  for  systems  of  different 
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isotopic  composition.  Clearly,  only  certain  colliding  pairs  could  then 
undergo  the  above  three-photon  event.  For  such  systems  the  final  translational 
energy  has  been  increased  by  an  amount  equal  to  the  energy  of  a  single  photon 
from  L2  (i  .e. ,  iiu^) .  From  the  previous  considerations  we  might  expect  such  an 
event  to  be  extremely  efficient  at  low  collision  energies.  The  extra  transla¬ 
tional  energy  picked  up  would  be  very  much  greater  than  the  relative  translational 
energy  of  the  unexcited  molecules  in  the  beam,  and  comparable  to  the  forward 
translational  energy.  Thus,  the  excited  molecules  will  tend  to  be  scattered 
out  of  the  beam.  Experimentally  the  collision  partners  can  be  prepared  in  the 
same  beam  which  is  formed  by  supersonic  expansion.  At  some  point  downstream 
of  the  source,  where  the  kinetic  energy  of  relative  motion  is  extremely  low 
(  IK),  the  beam  may  be  crossed  by  the  two  lasers  LI  and  L2.  The  desired  products 
which  are  scattered  out  of  the  beam  may  then  be  skinrod  off  and  collected. 

We  have  used  the  formalism  of  the  previous  section  to  calculate  the  import¬ 
ance  of  this  effect  for  the  XeF  system.  Let  us  define  a  separation  p  for 
separation  of  the  heavy  isotope  ^Xe: 
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where  o^g  and  o^q  are  the  kinetic-energy-averaged  cross  sections  for  Xe 

1 30  19 

and  Xe  collisions  with  F,  respectively.  When  L2  is  tuned  to  the  appro- 
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priate  frequency  for  resonance  in  Xe  F,  p  is  found  to  be  as  high  as  0.996. 

While  certain  aspects  not  taken  into  account  in  our  calculations  might  reduce 

the  efficiency,  such  as  finite  laser  bandw'dths  and  possible  elastic  collisions 

experienced  by  the  deflected  molecules  before  leaving  the  beam,  it  is  clear  that 

these  velocity-changing  radiative  collisions  could  provide  an  efficient  means 

of  isotope  separation.  Furthermore,  the  high  cross  section  of  the  event  (about 
°2  73 

20  A  in  our  specific  example  )  is  achieved  with  the  relatively  lower  laser 

2 

intensity  of  5  kW/cm  . 

The  advantage  such  a  process  has  over  more  conventional  techniques  like 
photodissociation  or  photoionization  is,  in  the  main,  one  of  flexibility:  first, 
there  is  almost  unlimited  control  over  the  initial  scattering  energy,  and 
second,  depending  on  the  tunability  of  L2,  there  is  a  large  choice  of  bound 
intermediate  states.  This  flexibility  allows  for  optimization  of  conditions. 

The  other  techniques  are  limited  to  using  initial  states  which  belong  to  a 
limited  number  of  bound  states.  Under  these  circumstances,  the  available  dis¬ 
sociation  pathways  may  often  be  unfavorable,  requiring  intense  lasers.  We  would 
like  to  add  that  a  chemically  reacting  system  (such  as  was  outlined  in  the  pre¬ 
vious  section)  might  be  even  more  appropriate  for  isotope  separation.  In  this 
instance,  only  one  laser  is  required  to  drive  the  reaction.  Finally,  we  sum¬ 
marize  by  noting  that  all  of  the  ideas  we  have  discussed  in  Part  B  of  Section 
II  lead  to  our  belief  that  laser-induced  chemistry  is  eminently  viable  provided 
one  stresses  resonance  formation  and  low  collision  energies. 
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III.  Surface  Processes 

While  the  influence  of  laser  radiation  on  homogeneous  gas-phase  sys- 
terns  has  been  intensely  studied,  ’  much  less  has  been  done  on  heterogeneous, 
e.g.,  gas/solid,  systems  either  theoretically35’74"93  or  experimentally.94"101 
The  effects  of  laser  radiation  on  adspecies/surface  systems  depend  upon  the 
nature  of  the  surface  (metal,  insulator  or  semiconductor,  smooth  or  rough,  etc.),  the 
coherent  properties  of  the  laser  radiation  (frequency,  intensity,  pulse  duration, 
polarization,  etc.)  and  the  electronic  and  vibrational  structure  of  the  adspecies/ 
surface  system.  Depending  upon  the  laser  frequency,  the  electronic,  vibrational 
or  vibronic  states  of  the  adspecies  and/or  the  substrate  may  be  perturbed  with 
susequent  energy  flow  among  the  various  modes  of  the  system.  In  this  section, 
we  shall  focus  on  the  infrared  vibrational  excitation  of  the  adspecies/surface 
system  in  which  the  direct  electronic  excitation  of  the  adspecies,  which  could 
only  be  affected  by,  e.g.,  visible  radiation,  will  be  ignored.  However,  indirect 
types  of  excitations  such  as  substrate  heating,  migration  and  desorption  of  the 
adspecies  will  be  considered. 

We  shall  first  investigate  the  field-free  system  where  species  (atoms  or 
molecules)  are  adsorbed  (physisorbed  or  chemisorbed)  on  a  solid  surface.  The 
combined  interface  system,  adspecies/surface,  is  significantly  different  from  the 
homogeneous,  i.e.,  gas-phase  or  clean  surface,  systems  as  follows: 

(i)  The  adspecies  (consisting  of  n  atoms)  has  3n  so-called  "frustrated"  libra- 
tional  degrees  of  freedom  instead  of  3n-6(5)  as  in  the  gas  phase.  For  example, 
there  are  six  vibrational  modes  for  CO  adsorbed  on  a  solid  surface  whereas  the 
surface-free  CO  molecule  has  only  one  vibrational  mode. 

(ii)  Due  to  the  substrate,  the  excited  adspecies  usually  has  a  shorter 
lifetime  as  compared  with  that  in  the  gas  phase,  because  of  its  strong  coupling 
to  the  phonon  modes,  particularly  for  chemisorbed  species. 
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(iii)  The  energy  spectrum  of  the  species  is  changed  in  the  presence  of  the 

substrate.  The  frequency  shift  and  level  broadening  may  be  caused  by  many  surface- 
induced  phenomena  such  as:  (1)  energy  flow  between  the  adspecies  and  the 

substrate  and  among  the  internal  modes  of  the  adspecies;  (2)  fluctuation  of 

the  effective  dipole  moment  of  the  adspecies  caused  by*  e.g.,  the  charge 
transfer  between  the  adspecies  and  the  adsorbent  surface  and  the  dynamic 
coupling  between  the  adspecies  dipole  and  its  image;  (3)  fluctuation  of  the 
conformations  of  the  adspecies  caused  by  the  libration  of  the  adspecies, 
particularly  when  the  system  is  subjected  to  laser  radiation;  (4)  migration- 
induced  broadening  via  inelastic  scattering  among  the  adspecies  (particularly 
for  physisorbed  states)  and  their  collision  with  the  surface  atoms;  and  (5) 
other  phase-relaxation-induced  broadening,  i.e.,  broadening. 

( i v )  The  quasi-continuum  region,  available  in  the  gas  phase  only  for 
polyatomic  molecules,  may  be  usually  achieved  even  for  diatomic  molecules  when 
adsorbed  on  a  solid  surface.  The  frustrated  degrees  of  freedom,  as  defined  in 
(i),  have  very  high  densities  of  states,  and  the  overlap  of  these  individual 
states  will  lead  to  a  quasi -continuum  of  the  overall  vibrational  potential  of 
the  adspecies  and/or  the  surface  potential. 

When  the  adspecies/surface  system  is  subjected  to  IR  radiation,  item  (iii) 
above  will  most  likely  play  a  particularly  important  role,  due  to  processes  such 
as  laser-enhanced  migration  and  desorption.  Furthermore,  when  the  adspecies  is 
highly  excited  and  reaches  the  quasi-continuum,  further  excitation  resulting  in 
bond-breaking  (desorption)  may  be  achieved  by  quasi-resonant  laser  excitation 
in  this  region  and/or  thermal -phonon-assisted  excitation.  The  concept  of  the 
quasi-continuum  is  essential  for  understanding  the  mechanism  of  desorption; 
without  this  concept,  one  must  resort  to  unrealistically  high  laser  powers  to 
overcome  the  anharmonicity  of,  e.g.,  a  Morse  potential.  Detailed  discussion 
on  this  aspect  will  be  given  later. 
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The  influence  of  IR  radiation  on  the  adspecies/surface  system  is  more 
complicated  to  describe  than  for  gas-phase  unimolecular  dynamics  as  described 
in  Part  A  of  Section  II.  The  following  important  mechanisms  will  be  addressed: 

•  How  is  the  laser  photon  energy  selectively  transferred  to  the  active  mode 
of  the  adspecies  with  consequent  relaxation  to  the  bath  modes? 

•  What  is  the  energy  relaxation  dynamics  and  how  are  the  (single-phonon 
and  multiphonon)  relaxation  rates  governed  by  systems  parameters,  such  as 
frequency  spectrum,  interaction  potential  and  masses  of  the  adspecies/surface 
system? 

•  What  is  the  nature  (selective  and  nonselective)  of  the  photon  energy 
deposition  in  the  system?  How  is  the  selectivity  characterized  by  the  damping 
rate,  pumping  rate  and  the  rate  of  intramolecular  vibrational  relaxation? 

•  How  does  the  nature  of  energy  flow  from  the  excited  mode  to  the  other 
modes  depend  upon  mathematical  treatments,  e.g.,  Markoff  approximation, 
conti nuum-density-of-states  approximation  and  the  feedback  mechanisms  achieved 
by  memory  effects  and/or  the  substrate  heating  effects? 

•  How  are  the  isotope  effects  for  a  system  consisting  of  a  mixture  of 
isotopic  adspecies  characterized  by,  e.g.,  the  frequency  difference  of  the 
isotopes,  the  direct  dipole-dipole  interaction  and  the  phonon-mediated  coupling? 

•  How  does  the  laser  radiation  affect  the  rate  processes  via  field-induced 
dynamics  such  as  desorption,  migration,  decomposition,  recombination  and 
desoption?  In  other  words,  how  do  the  laser  and  surface  synergistical ly  affect 
heterogeneous  catalytic  reactions? 

None  of  the  above  addressed  aspects  of  laser-induced  surface  processes 
(LSSP)  has  been  completely  investigated,  either  experimentally  or  theoretically. 
However,  some  qualitative  and/or  semi-quantitative  results,  will  be  discussed 
in  this  section,  which  consists  of  two  parts:  (A)  energy  flow  in  adspecies- 
surface  systems  where  microscopic  treatments,  types  of  selectivity,  photon 


energy  population,  isotope  effects  and  dynamic  effects  (migration  and  desorption) 
will  be  investigated;  and  (B)  laser/surface-catalyzed  reactions  where  applications 
of  the  theoretical  aspects  of  LSSP  and  some  recent  experimental  work  on  LSSP 
will  be  reviewed. 
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A.  Energy  Flow  in  Adspecies-Surface  Systems 

Considering  a  group  of  species  (atoms  or  molecules)  adsorbed  on  a  solid 
surface  and  subjected  to  IR  radiation,  the  most  general  total  Hamiltonian 
describing  the  energy  flow  within  the  adspecies/surface  system  may  be  written  as 


H  «  Hn  + 


ifj 


Hij 


l  H 


1  1 


Hiji 


i/Lk 


ikj 


HiF 


(t) 


(III. 1) 


Hg  represents  the  unperturbed  Hamiltonian  of  the  system  (adspecies,  substrate 
phonons  and  laser  photons);  H. .  and  Hikj  represent  the  direct  interaction  and 
the  k-mode-mediated  interaction  between  the  i-th  and  j-th  modes,  where  i,  j  and  k 
may  be  the  adspecies,  adbond  or  phonon  modes;  and  H^p  is  the  interaction 
between  the  i-th  mode  of  the  system  and  the  field  (laser  photon).  In  the 
following  discussion,  we  shall  assume  only  one  active  mode  in  the  adspecies  to 
be  excited  by  the  IR  radiation,  except  for  when  we  discuss  isotope  effects, 
in  which  case  there  will  be  two  active  modes. 

The  microscopic  Hamiltonian  including  heat- bath  effects  will  be  used  to 
calculate  the  excitations  for  both  single-phonon  and  multi  phonon  processes.  A 
more  complicated  system  consisting  of  A,  B  and  C  subgroups  will  be  investigated, 
where  the  memory  function  of  the  interaction  between  A  and  B  and  their  common 
damping  to  C  is  included.  Energy  populations  in  a  multilevel  system  will  be 
studied,  and  the  mechanism  of  excitation  (selective  vs  nonselective)  will  be 
numerically  examined.  Numerical  results  will  also  be  shown  for  the  adbond 
excitations  in  the  H/W  system.  By  solving  the  rate  equation  in  energy  space 
in  the  quasi-continuum,  several  types  of  photon  energy  populations  and  the 
associated  desorption  probabilities  will  be  calculated.  Extension  of  the 
single  adspecies  to  a  system  consisting  of  two  isotopes  will  be  investigated, 
and  finally,  the  dynamic  effects  (migration  and  desorption)  will  be  studied. 

While  much  of  the  work  presented  in  this  section  has  already  been  published  or 
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is  in  press  (with  references  cited),  some  unpublished  concepts  and  results  are 
also  included,  e.g.,  the  discussions  in  Parts  A. 2. a,  A. 4  and  A. 5. 

1 .  Microscopic  Treatment 
a.  Single-Phonon  Processes 

Consider  a  heterogeneous  system  consisting  of  atoms  or  molecules  adsorbed 
on  a  uniform  solid  surface  and  subjected  to  infrared  laser  radiation.  The 
vibrational  degrees  of  freedom  of  the  adspecies-surface  system  can  be  divided 
into  two  groups,  namely,  the  pump  mode  (resonant  infrared-active  vibrational 
mode  of  the  selectively-driven  adspecies)  and  the  bath  modes  (all  other  modes 
including  adspecies  inactive  modes  and  the  surface  phonon  modes).  The  radia¬ 
tion  feeds  energy  into  the  pump  mode,  and  the  heat  bath  provides  a  relaxation 
mechanism.  The  microscopic  model  Hamiltonian  describing  this  relaxation 
dynamics  can  be  written  in  the  following  second-quantized  form,  corresponding 

82 

to  the  simplest  case  of  Eq.(III.l)  consisting  of  only  HQ,  •  and  H..p  with  i=A: 


H  =  HA  +  HB  +  HAB  +  H'(t>  • 

(III. 2a) 

1.  =-f!u,ata  +  I-tlB  (a++a)p  , 

*  *  p*3  p 

(III. 2b) 

HB  =  H“jbjbj 

J 

(III. 2c) 

HAB  ■  H^V^ahj)  • 

J 

(III. 2d) 

H'(t)  ='tTVcos(ut)  (a++a)  , 

(III. 2e) 

V  =  (21TmAuA)'V(0)E  cos ( e) 


(III .2f ) 
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H^  and  Hg  are  the  unperturbed  Hamiltonians  (vibrational  energy)  of  the  pump 
mode  and  the  bath  modes,  respectively;  H^g  is  the  interaction  Hamiltonian 
coupling  the  pump  mode  and  the  bath  modes;  and  H'(t)  is  the  adspecies-field 
effective  interaction  Hamiltonian.  The  operators  a+,a  and  b+,b  are  the  usual 
harmonic  vibrational  ladder  operators  (with  fundamental  frequencies  and 
of  the  pump  mode  and  the  bath  modes,  respectively;  6p  is  the  anharmonicity  of 
the  nonlinear  quantum  oscillator  (pump  mode);  u'(0)  is  the  derivative  of  the 
effective  dipole  moment  of  the  pump  mode  evaluated  at  the  equilibrium  point; 

E  is  the  electric  field  of  the  radiation  with  frequency  u>,  linearly  polarized 
at  an  angle  e  with  respect  to  the  effective  dipole  moment;  and  m^  is  the  reduced 
mass  associated  with  the  active  mode.  In  Eq.  (III. 2d),  K.  represents  the  coupling 

J 

between  the  active  mode  and  the  j-th  bath  mode.  We  note  that  the  laser  field 
defined  in  the  expression  of  Eq.  (III. le)  was  assumed  to  be  a  single-mode  classical 
field;  the  case  of  a  quantized  multimode  laser  excitation  will  be  investigated 
later.  [We  should  also  mention  that  8*0  in  Eq.  (II. 12e) ,  which  is  why  the  cosine 
term  is  absent  there.] 


The  equation  of  motion  for  an  operator  0(t)  in  the  Heisenberg  picture, 

0(t)  =  exp(iHt/iT)  0  exp(-iHt/tf)  ,  (III. 3) 

where  the  time- independent  operator  0  is  defined  in  the  Schrodinger  picture,  is 

=  4,  [0(t)  ,  H]  ,  CHI-4) 


It  Tff 

14 

employing  the  operator  algebra 

[a(t)  ,  H]  =  3H/3a+(t)  , 

[b(t)  ,  H]  =  3H/3b+(t) 

we  obtain  the  following  set  of  coupled  equations: 

a(t)  =  iweffa(t)  -  i  J  Kjbj(t)  -iVcos(wt), 

J 


(III. 5a) 
(III. 5b) 

(III. 6a) 


bj(t)  =  -iujbj(t)  -iKja(t) 


(III. 6b) 
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where  weff  is  the  effective  frequency  obtained  by  the  contact  transformation,82 
including  anharmonicity  up  to  fourth  order  in  Eq.  (III. 2b), 

“eff(t)  =  “A‘2e*a+{t)a(t)  *  (III.6c) 


e*  =  3063/^  -  6S4 


(111. 6d ) 


We  solve  for  the  phonon  operators  b.(t)  by  formally  integrating  Eq.  ( I I I . 6b) 

J 

to  obtain 


bj(t) 


bj(Q)  exp(-iwjt) 


rt 

dt1  a(t')  exp[-1iD.(t-t')]. 

J 

0 


(III. 7a) 


Substituting  Eq.  (III. 7a)  in  Eq.  (III. 6a),  we  get 


a(t)  =  iwe^(t)a(t)  +  Aj^(t)  +  A2(t)  -  iV  cos(ut)  , 

(III. 8a) 

where 

A1(t)  =  -i  l  Kjbj(0)exp(-ia)jt) 

J 

ft 

(III. 8b) 

A2(t)  =  -  j  dt’A(t')  M(t-t') 

(111.8c) 

u 

and  M(t)  is  the  memory  (kernel)  function  given  by 

M(t)  *  l  iKj|2exp(-i<Djt)  . 

(III- 8d ) 

To  solve  the  above  equation  and  hence  the  average  excitation,  the  unknown 
function  A(t)  and  the  many-body  effects  due  to  the  phonon  modes  in  A2(t)  must  be 
decoupled.  Depending  on  the  memory  function  M(t),  which  is  characterized  by  the 
time  evolution  of  the  active  mode  and  the  density  of  states  (or  frequency  spectrum) 
of  the  phonon  modes,  the  excitation  and  relaxation  may  behave  quite  differently. 

For  example,  the  relaxation  of  the  active  mode  may  be  irreversible  in  one  way  or 
reversible  in  another,  i.e.,  there  may  be  energy  feedback  to  the  active  mode  in 
the  other  way.  There  are  several  techniques  to  deal  with  the  aforementioned 
many-body  problem.  We  shall  first  discuss  those  for  irreversible  relaxation,  and 
the  feedback  mechanism  resulting  from  different  features  of  the  memory  function 
will  be  investigated  later. 


62 


Within  a  Markoff  process  where  the  correlation  time  of  the  phonon  opera¬ 
tors  is  much  faster  than  that  of  the  active  mode  and  the  excitation  rate,  the 
zeroth-order  active-mode  operator  a(t)=a(0)exp[-i  dt'c^U' )]  may  be  used  to 
evaluate  AgU).  There  are  several  alternatives  to  evaluate  the  integral  in 
Eq.  (II I. 8c):  (i)  extending  the  upper  time  limit  from  t  to  infinity  and  working 
out  the  time  integral  first;  (ii)  assuming  a  continuum  phonon  spectrum  and  fac¬ 
toring  out  the  coupling  strength  K  and  the  density  of  states  p;  (iii)  assuming 

14 

the  memory  function  is  a  delta  function  in  time.  Either  of  these  methods  gives 

-■ -  (y^/2  +  i5w)a(t)  ,  (III. 9) 

where  the  damping  factor,  Yj»  and  the  frequency  shift,  6co,  are  given  by 

YX  *  2ir|K(o>A)|2p(o)A)  , 

=p[  l  |  Kj  |2/(uA-o,j)]  , 

«J 

where  P  stands  for  "principal  part". 

Substituting  the  result  of  Eq.  (III. 9)  into  Eq.  (III. 8a),  we  decouple 
the  phonon  modes  and  obtain  a  single  equation 

a(t)  »  i(o,eff+5u)-iY1/2)a(t)  +  Aj(t)  -  iV  cos(wt)  .  (III. 11) 

The  effects  of  the  phonon  modes  have  thus  been  incorporated  into  a  damping 
factor,  y^»  and  a  frequency  shift,  6u>,  of  the  active  mode.  Another  technique 
to  treat  the  multimode  phonon  effects  is  the  Wigner-Weisskoff  single-pole 
approximation,  which  we  shall  discuss  later  when  dealing  with  multiphonon 
relaxation  processes. 


(III. 10a) 
(III. 10b) 
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The  above  nonlinear  equation  cannot  be  solved  analytically  due  to  the 
excitation-dependent  anharmonicity  in  u>eff(t)  [see  Eq.  (III. 6)].  However,  for 
low  excitation,  we  may  ignore  the  anharmonicity  and  solve  Eq.  (I II. 11)  to  ob¬ 
tain  the  ensemble-averaged  fan  the  active-  and  bath-mode  coordinates)  excitation 

‘t* 

<n(t)>  i  «a  U)a(t)>>  for  the  harmonic  case  (e*=0). 


<n(t)> 


(V/2)2 

..  a^+(Yi/2)2 


' 

-Y,t  -y, t/2 

1+e  -2e  cos(At) 

* 

_Ylt 
1-e  1 

+  n 

(III. 12) 


where  is  the  phonon-shifted  detuning  and  n  is  the  phonon  occupation 

number  given  by  the  Bose-Einstein  function,  evaluated  at  the  active-mode  frequency. 


n  » 


exp(ficuA/kT)-l 


(III. 13) 


Here  we  assumed  that  the  phonon  modes  provide  an  infinite  energy  "sink"  and  the 
initial  surface  temperature  (T)  remains  constant.  This  is  not  necessarily  always 
true  when  one  considers  a  more  complicated  system  consisting  of  several  sub¬ 
groups  of  bath  modes,  where  some  of  the  subgroups  may  have  finite  numbers  of 
modes  rather  than  an  infinite  sink.  Further  discussion  on  this  will  be  given 
later. 


So  far  we  have  discussed  only  the  situation  where  the  damping  or  level 
width  of  the  excitation  is  governed  by  the  so-called  T^(energy)  relaxation 
rate,  Yj*  To  include  the  effects  of  Tgtphase)  relaxation  on  the  excitation, 
we  investigate  the  ensemble-averaged  (over  the  phonon-bath  coordinates)  equa¬ 
tion  of  motion  for  the  active-mode  operator,  0(t)=a(t)  or  a+(t)a(t),  in  the 
Heisenberg-Markoff  picture  (HMP).^’®* 


d<Q(t)> 

dt 


[iVcos(ajt)/ti] 


30(t) 

3a+(t) 


+  a+  »3(t) 

aa+(t) 


+  Y1n 


d20(t) 

aa(t)sa+(t) 


(III. 14) 
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The  last  term  involving  the  dephasing  (T2  processes)  is  characterized  by  the 
dephasing- induced  broadening  factor  y2  as  follows: 


=  y2  ^  a+(t)a(t),0(t)ja+(t)a(t)  -  a+(t)a(t)  ^a+  ( t ) a ( t )  ,5(t)  ^  , 


(III. 15) 


which  is  mathematically  constructed  such  that,  for  0(t)=a(t)  and  a+(t)a(t). 


MiT> 


-  -  r2a(t). 


d[at(t)a(t)]\ 
at  /T, 


=  0  . 


(III. 16) 


This  assures  that  the  T2  dephasing  changes  only  the  phase  of  the  active  mode 
without  changing  its  vibrational  energy.  By  analogy  with  the  above  phenomena, 
in  colli  si onal  phenomena  the  T^  and  T2  relaxation  correspond  to  inelastic  and 
elastic  scattering,  respectively,  and  the  overall  collisional  broadening  is 
then  given  by  y^+y2. 

By  using  Eq.  (III. 14)  and  the  rotating-wave  approximation,  the  ensemble- 
averaged  equations  of  motion  in  the  HMP  are  found  to  be 

<a(t)>  =  -  [i (A-2e*<a+(t)a(t)>)  +  (Yi+y,)/2]  <'a(t)>-iV/2  ,  (III. 17a) 


<n(t)>  =  (iV/2)  <<a+(t)-a(t)»  -  y^(<n(t)>-n) 


(III. 17b) 


4 
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i 


where  we  have  used  the  rotating  frame  defined  by  a(t)=a(t)exp(iwt) .  For  long 
laser-pulse  excitation,  we  may  employ  the  adiabatic  limit,  | a | << ( y^+yg ) |a| ,  to 
obtain  the  energy  rate  equation  from  Eq.  (III. 17)  with  a=0, 

tt  ■  &  -  n  •  OH-18) 


where  a  and  I  are  the  absorption  cross  section  of  the  whole  system  (active  plus 
bath  modes)  and  laser  intensity  (with  photon  energy^).  aI/>fJu  therefore  denotes 
the  total  excitation  rate  of  the  system  given  by  81 


AI(yj+y2) 

U-2e*<n>]  “>(y  ^2^/4 


(III. 19) 


2  ? 
and  A=  ( V/E )  (4ir/c)  is  a  constant  resulting  from  I=E  /(8tt/c). 

The  steady-state  average  excitation,  X,  for  a  cold  surface,  i.e.,  neglecting 

n,  is  seen  from  Eq.  (II I. 18)  to  be  governed  by  a  cubic  equation, 


x  _  Al(r+/Yl)  ,  (III. 20) 

(A-2e*X)2+(r+/2)2 

where  r+=Y1+Y2  is  the  total  damping  factor  resulting  from  the  Tj  (energy)  and 
7  2  (phase)  relaxation.  The  optimal  detuning  then  occurs  at  the  maximum  (dX/dA  »  0) 
and  is  given  by  A*  =  2e*X*  (note  -  the  single  asterisk  which  was  already  attached 
to  e  does  not  signify  an  optimal  condition  as  it  does  for  A  and  X).  At  the  other 
extreme  where  dX/dA  -*■  ~,  we  obtain  a  quadratic  equation  for  the  detuning,  whose 
two  roots  correspond  to  a  "bistability"  in  X  as  a  function  of  A.  By  equating 
the  two  roots,  we  obtain  the  critical  pumping  rate  I*=Yir+/(8Ae*) ,  implying  that 
the  existence  of  the  bistability  is  a  consequence  of  the  condition  I>I*.  For 
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p 

a  fixed  laser  intensity,  which  is  proportional  to  V  (or  the  pumping  rate), 
the  bistability  criterion  may  also  be  stated  in  terms  of  the  anharmonicity  as 
e*  >  e**  =  ( y^T+)/( 4AI ) This  "bistability"  feature  of  the  steady-state 
excitation  is  shown  in  Fig.  11.  It  is  seen  that  when  the  anharmonicity  e*  is 
larger  than  the  critical  value,  e**,  the  excitation  profile  shows  the  bistable 
transition  from  P  to  Q  as  the  detuning  increases,  and  from  R  to  S  as  the  detuning 
decreases.  We  note  that  the  maximum  excitation  is  red- shifted  to  a*  >  0,  which 
is  a  general  property  of  any  nonlinear  oscillator  with  e*  >  0.  A  classical 
analogy  of  this  nonlinear  quantum  oscillator  will  be  shown  later. 

A  remark  on  the  relaxation  rate  of  the  single-phonon  processes  must  be 
made  before  we  discuss  multi  phonon  processes.  We  have  shown  that  the  many-body 
effects  of  the  phonon  modes  on  the  excitations  of  the  active  mode  may  be  reduced 
to  an  energy  relaxation  factor  (y^)  and  frequency  shift  (5u>)  as  shown  by  Eq. 

(III. 10).  In  the  Markoff  approximation  the  single-phonon  relaxation  rate  is 


given  by 


14,78 


oo 

Yj  =  Re  ^  |K(ojj)|2  ^dx  exp[-i(uj-uA)x] 


00 

I 


doJj  |  K(uj)  |  2p(u>j)TT<$(w..-a>A) 


(III. 21) 


which  vanishes  when  a  situation  in  which  the  active-mode  frequency  is 

usually  greater  than  that  of  the  bath  modes.  This  unrealistic  zero  relaxation 

rate  may  be  avoided  if  one  includes  the  finite  lifetime  of  the  phonon  modes,  y^, 

87 

which  turns  the  delta  function  in  Eq.  (III. 21)  into  a  Lorentzian.  For  an  Einstein 
spectrum  p(w.)  =  6(u>.-u>D),  we  obtain 


1% 

(“a) 

V2 

\ms  / 

W) 

^a'wd)2+(yb/2)2. 

(III. 22) 


A 
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which  is  significantly  different  from  the  result  for  zero  phonon  width  =  0. 

2  3 

For  a  Debye  model,  p(wj)  =  3u>j/tdg  and  yg=  0,  we  have 


(III. 23) 


where  m^  and  m^  are  the  reduced  masses  of  the  adspecies  and  the  surface  atoms, 

..  ,  78,93 

respectively. 

b.  Multi  phonon  Processes 

We  have  shown  that  a  single-phonon  relaxation  rate  may  be  achieved  by  the 
finite  lifetime  (or  level  width)  of  the  phonon  modes  if  the  frequency  difference 

is  not  significantly  greater  than  the  phonon  width  yg.  This  is  for  the 
situation  that  the  active-mode  frequency  is  close  to  the  phonon  spectrum.  How¬ 
ever,  for  an  IR-active  system  the  active-mode  frequency  is  usually  much  higher 
than  that  of  the  phonon  modes,  e.g.,  for  CO/Ni,  u>A~2irxl000cm~  *,  toj.-2irx300cm"1 
and  Yg~50cm"*,  which  results  in  a  very  small  single-phonon  relaxation.  Thereto/*, 
for  those  systems  with  a  big  energy  gap  between  the  active  and  bath  modes,  the 
relaxation  mechanism  is  dominated  by  multi  phonon  processes  and  neither  Eq.  (II 1. 22) 
nor  (III. 23)  is  an  appropriate  expression  for  yj.  To  fint.  the  power  law  for  the 
ratios  (mA/ms)  and  (u)A/wD),  we  shall  investigate  the  thermally-averaged  (on  both 
the  active-  and  bath-mode  coordinates)  golden  rule  for  the  multi  phonon  relaxation 
rate,87 


00 

rl  =^rjRe^dt  hab^°^  ’ 


(III. 24) 
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which  is  governed  by  the  autocorrelation  function  of  the  interaction  Hamiltonian 
between  the  active  mode  and  the  bath  modes. 

Consider  a  Morse  surface  potential 


V(Q.U)  =  0  e 


-2«(Q-U-r0)  (Q.0) 


=  2aOQ[e2“U  -  e^] 


(III. 25a) 


(III. 25b) 


where  D  and  a  are  the  appropriate  parameters  of  the  Morse  potential,  and  Q  and 
U  are  the  position  of  the  adatom  and  of  the  surface  atom  with  equilibrium 
position  rg.  Since  the  frequency  of  the  adatom  is  much  higher  than  that  of 
the  phonon  spectrum,  we  have  used  the  linear  expension  exp  (aQ)  »  1  +  aQ  and 
dropped  the  constant  terms  which  do  not  conserve  energy  in  Eq.  (III.  25b).  The 
interaction  Hamiltonian  is  then  given  by  the  fluctuation  of  the  adatom-surface 


potential,  i.e., 


H«d  =  V(Q»U)  -  <V(Q,U)> 


(III. 26) 


where  <V(Q,U)>  is  the  thermally-averaged  (over  the  bath-mode  coordinates) 
surface  potential.  Employing  the  identities 


JaU(t)eau(0)>  =  <eoU(t)><eaU(0)>e<aZU(t)U(0)> 


(III. 27a) 


<eaU>  =  e<a2U2/2> 


(III. 27b) 


we  obtain  the  autocorrelation  function 


<<HAB(t)HAB(0)>>  =  (2a0)2«Q(t)Q(0)»  {EQ(Et+l)  +  E*(eJ+1)-2E*/2(E2+1  )} , (III. 28a) 
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=  exp(a2<<U2>>) 


(III. 28b) 


Et  =  exp(«  «U(t)U(0)»)  .  (111.28c) 

Eq.  (III. 28)  together  with  Eq.  (III. 24)  enables  us  to  evaluate  the  multiphonon 

relaxation  rate  up  to  any  order.  Note  that  the  linear  term  of  the  potential 
2 

V(Q,U)  =  4a  DQU  gives  us  the  single-phonon  relaxation  rate,  Eq.  (III. 22),  when 
2  2 

2a  D=m^.  Furthermore,  the  relaxation  rates  may  behave  differently  for  different 

choices  of  the  interaction  potential  and/or  phonon  spectrum.  To  evaluate  the 

multiphonon  relaxation  rate  for  a  very  strong  coupling  case,  the  steepest 

87  89 

descent  method  is  usually  used. 

For  a  p-phonon  process  with  a>^*qwq+(p-q)ujp,  where  and  wp  are  the 

frequencies  of  the  active  mode  and  of  the  acoustic  and  optical  phonons,  respecti¬ 
vely,  the  interaction  Hamiltonian  includes  the  term  [a2^p+1^QUp/p!].  Transforming 

the  coordinates  Q  and  U  into  the  second  quantization  operators,  the  integral  of 

77 

Eq.  (I II. 24)  may  be  worked  out  for  an  Einstein  spectrum  to  give 


2  /  \q, 

>  f  UA  \  f _ ^ 

>*  (p-q)  n  4 \qwqy  ^(p-q) 


vvVv 


(III. 29) 


where  m.,  m  and  m  are  the  reduced  mass  of  the  active  mode  and  of  the  acoustic 
A  q  p 

and  optical  phonons.  We  note  that  the  multiphonon  relaxation  rate  is  strongly 

decreasing  with  increasing  the  order  parameters  p  and  q,  since  mA<<mn’mn*  There- 

h  P  q 

fore,  for  systems  with  a  small  mass  ratio,  ni./m  ,m  «l,  and/or  a  high  active- 

m  q  q 

mode  frequency,  u«>>u_,u)  ,  we  expect  a  long  lifetime  of  the  excited  adspecies. 
a  p  q 

The  above  results,  which  provide  us  with  a  quantitative  description  of  the  power 

laws  of  the  mass  and  frequency,  are  in  agreement  with  that  of  the  energy-gap 
37 

law.  The  relaxation  rate  for  the  N/Si  system  has  been  estimated  to  be 

Yl  a  lons-1,  where  for  H/Pt,  Yl  =  loV1.  The  latter  is  much  lower  due  to  the 

smaller  mass  ratio,  ar|d  =  2000  cm  wq  =  150  cm  and  wp  *  400  cm  . 
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To  obtain  the  average  excitation  of  the  active  mode  and  temperature 

dependence  of  the  multiphonon  relaxation  rate,  we  shall  consider  the  microscopic 
81 

Hamiltonian 

H  =  Hft  +  Hb  +  Hab  +  H'(t)  ,  (III. 30} 

which  is  the  same  as  that  of  the  single  phonon  case  except  that  the  interaction 
Hamiltonian  H^B  describing  multi  phonon  relaxation  is  now  given  by 

HAB  =  l  fr6Aa+  +«a)  ’  (III. 31a) 

v 

where  we  define  a  multi  phonon  operator 

J  v 

B  5  n  b.  b.  ---b.  ,  (III. 31b) 

V  jj.J2.--j,,  J1  J2  JN 

and  Gv  is  the  coupling  strength.  It  is  seen  that  H^g  reduces  to  Eq.  (III. 2d) 

and  describes  single-phonon  relaxation  when  N=1  and  G  =K..  The  Heisenberg  equations 

^  J 

of  motion  are 


a(t)  =  i[ueff(t)+u(t)]a(t)  -  i 

1  GyBv(t)-iV  cos(ojt) , 

V 

(111.32a) 

B  (t)  =  -in  B  (t)-iN  G*a(t) , 
vv  vv 

(III. 32b) 

where 

L 

%  -  CBv(t).HB]/Bv(t) 

=  I  “j  » 

(III. 33a) 

Nv  -  CBvtt).B+{t)3  . 

(III. 33b) 

and  5s(t)  is  a  stochastic  frequency  modulation  which  accounts  for  dephasing  effects. 

By  employing  the  Markoff  approximation  as  used  in  the  single-phonon  case, 

Eq.  (III. 32)  is  decoupled  and  results  in  the  ensemble-averaged  equations  of  motion 
in  the  HMP14,84 

<a(t)>  =  -Ci(A-2e*<a+(t)a(t)>)  +  (Y1+Y2)/2]4(t)>-iV/2  ,  (III. 34a) 

<n(t)>  =  -(iV/2)<<a+(t)-a(t)»  -  Yj(<n(t)>-N)  ,  (III. 34b) 

which  have  the  same  structure  as  for  the  single-phonon  case  [Eq.  (III. 17)], 

except  that  we  now  have  the  multiphonon  equilibrium  occupation  N  *  nn.,  and  the 

j  J 
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relaxation  factor  is  given  by 

YX  =  2ir  l  |Gv|2Nv6(u>A-nv)  ,  (III. 35a) 

v 

where  .  . 

JN  JN 

N  =  n  (n.+l)  -  n  n..  (III. 35b) 

In  deriving  Eq.  (III. 35),  the  phonon-induced  frequency  shift  was  neglected 
and  we  assumed  that  the  stochastic  frequency  obeys  the  simple  correlation 
<aj(t)aj(tl  )>  *  We  note  that  the  significant  difference  between  single 

phonon  and  multiphonon  processes  lies  in  the  nature  of  the  relaxation  factor, 

Yj,  which  is  temperature  independent  [Eq.  (III. 21)]  for  single-phonon  relaxation 
but  strongly  temperature  dependent  for  multi  phonon  relaxation.  For  example, 
for  an  Einstein  spectrum  with  p ( oi ^ )  =»6 ( to ^ -toD )  and  fiv=pa)D  (p-phonon  process),  we 
find  the  temperature  dependence  of  the  relaxation  rate  as 

Yl  «  [exp(1io)A/kt-l]  /  [exp(fia>D/kt)  -l]p  ,  (III. 36) 

which,  for  low  temperature  (kT<«fTwg)  becomes  almost  temperature  independent, 

Yj  -  1  +  pexp(  fiwD/kT)  ,  (III. 37a) 

whereas  for  high  temperature  (iTuiA>kT>yffwD)  ,  (III. 37b) 

YX  «  (kT/%!D)p[exp(1Tu>A/kT)-l]  , 
which  shows  a  strong  temperature  dependence. 

We  have  also  assumed  the  coupling  factor  to  be  independent  of  tempe¬ 
rature  and  the  spectrum  of  the  phonon  modes  to  be  governed  by  a  delta  function. 
Furthermore,  only  the  single  term  in  Eq.  (III. 35)  which  provides  the  dominant 
contribution  to  the  relaxation  rate  is  considered.  When  considering  many  atoms 
of  the  substrate,  the  order  parameter  p  should  be  split  into  several  parts,  i.e. 
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o)^=  J  qojq,  with  Iq=p,  and  the  associated  compound  density  of  states  should  be 

36 

expressed  in  terms  of  a  convolution  of  the  single-phonon  density  of  states. 

To  evaluate  the  multi  phonon  relaxation  rate  for  an  actual  system,  a  realistic 

phonon  spectrum  is  needed,  and  information  on  the  crystal  structure  of  the 

102 

substrate  is  required  for  a  first-principles  calculation. 

We  have  so  far  investigated  the  microscopic  Hamiltonian  based  on  a 
single-mode  classical  field,  invoking  the  assumption  that  the  number  of  photons 
in  the  laser  field  is  so  large  such  that  the  electric  field  (or  intensity)  of 
the  radiation  is  kept  constant  even  in  the  presence  of  the  absorbing  medium. 
Depending  on  the  initial  states  of  the  system,  the  active  mode  may  absorb  or  emit 
photons  corresponding  to  the  excitation  and  de-excitation  of  the  adspecies. 
Similarly,  the  phonon  modes  may  behave  as  a  thermal  field  when  the  adspecies  is 
initially  "cooler"  than  the  bath  modes.  For  instance,  if  the  substrate  is  heated 
by  a  laser  or  other  thermal  source  but  the  adspecies  is  not,  then  energy  could 
flow  from  the  phonons  to  the  adspecies.  To  study  these  cooling  or  heating 
processes  of  the  adspecies,  we  shall  consider  a  microscopic  Hamiltonian  similar 
to  that  of  Eq.  (III.l),  except  that  the  adspecies-laser  interaction,  H'(t),  is 
non-quantized  into  a  multimode  field  form, 

H'(t)  £  (Vkcja  +  Vjcka+)  ,  (III. 38) 

4-  ^ 

where  ck  and  ck  are  the  harmonic  ladder  operators  for  the  k-th  mode  of  the 
quantized  field,  and  V k  is  proportional  to  the  electric  field  due  to  radation 
and  may  be  referred  to  as  the  Rabi  frequency  of  the  excitation.  The  Heisenberg 
equations  of  motion,  previously  given  by  Eq.  (III. 32),  now  become 

a(t)  »  i«o  ff(t)a(t)  -  i  l  G  B  (t)  -  i  l  V{ck(t)  ,  (III. 39a) 

v  k 

Bv(t)  =  -  10vBv(t)  -  iNvG*a(t)  ,  (III. 39b) 

cR(t)  =  1wkck(t)  -  iVka(t)  . 


(III.39c) 
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The  above  coupled  equations  can  be  solved  by  using  the  Markoff  approximation 
to  include  the  many-body  effects  due  to  the  phonon  modes  and  the  laser  as 
developed  in  the  previous  sections.  Here  we  shall  present  an  alternative 
technique. 

Taking  the  Laplace  transform  of  Eq.  (2.35),  we  get 

a(s)  3  U(s)a(s)  +  l  Vy(s)Bv(s)  +  l  Wk(s)ck(s),  (III. 40) 
v  k 

where  a(s),  Bv(s)  and  ck(s)  are  the  Laplace  transforms  of  a(t),  Bv(t)  and  ck(t), 
respectively,  and  U(s),  Vy(s)  and  Wk(s)  are  the  appropriate  functions  obtained 
from  the  Laplace  transform  of  Eq.  (III. 39).  Employing  the  Wigner-Weisskopf 
single-pole  approximation,  i.e.,  saO  in  Eq.  (III. 40),  the  inverse  Laplace  trans¬ 
form  of  Eq.  (III. 40)  can  be  obtained  to  evaluate  the  average  excitation  of  the 
active  mode, 

<n(t)>  s  «a+(t)  a(t)»  3  <n(0)>  e~rt  + 

2  V 

l  ,N-  -  i+e-rt-2e-rt/2cos(A  t) 
v  42+(r/2)2 
.  v  JL 

2 

+  l  ~  fl+e"rt-2e"rt/2cos(Akt)lnk(0)  , 

k  A2+(r/2)2  L  k  J  k 

J  (III. 41) 

where  <n(0)>,  n\(0)  and  ^(0)  are  the  initial  occupation  numbers  of  the  active, 

phonon  and  photon  modes,  respectively,  defined  by  a  Bose  function  with  the  same 

temperature  but  at  a  different  frequency.  The  detunings  are  defined  by 

A  and  A.  =o), -u.,  where  n  is  the  multi  phonon  frequency  and  u.  is  the  k-th 

vAvKAK  v  ^ 

mode  laser  frequency.  Finally,  r=v1+Yg  is  the  total  damping  factor  describing 
the  effects  of  the  parallel  nonradiative  (phonon)  and  radiative  (photon) 
relaxation  of  the  quantum  oscillator  into  two  independent,  noninteracting 
multimode  baths.  We  note  that  the  multi  phonon  rate  Yl  has  a  temperature  dependence 


Hn.(O) 

j  J 
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given  by  Eq.  (III. 36),  while  the  radiative  factor  Yq  is  independent  of  tempera¬ 
ture  since  we  consider  dipole  transitions  of  the  active  mode,  i.e.,  the 
interaction  Hamiltonian,  Eq.  (III. 38),  contains  only  linear  coupling  terms  in 
the  resonance  excitation  with  frequency 

Assuming  an  Einstein  spectrum  for  both  the  phonon  and  photon  modes,  the 
average  excitation  at  resonance,  ^=4^=0,  becomes 

2 

1  -  e  "rt/2l  +  <n(0)>e"rt,  (III. 42) 


which  is  characterized  by  the  Incoherent  phonon  field  and  the  coherent  laser 
field  with  initial  occupations  and  ppiy,  respectively,  and  the  initial 
occupation  of  the  active  mode  is  decaying  with  a  total  rate  r  through  two  parallel 
channels,  the  phonon  and  photon  bath.  We  note  that  the  active-mode  excitation, 
with  initial  state  <n(0)>=0,  is  built  up  by  both  the  laser  field  and  the  phonon 
field,  particularly  when  the  initial  phonon  occupation  is  high  and/or  is 
"heated  up"  by  the  laser  field  either  directly  or  indirectly  via  mode-mode 
coup! ing. 

2.  Selective  Excitation  and  Thermal  Effects 

So  far  the  adspecies/surface  system  has  been  described  microscopically  by 
a  Hamiltonian  in  which  the  active  mode  is  singled  out, with  all  the  other  (inactive) 
modes  treated  as  a  heat  bath.  The  many-body  effects  of  the  bath  modes  were 
reduced  to  the  damping  factors  Yj  and  y2  based  on  a  Markoff  approximation  in 
which  the  kernel  function,  M(t)  in  Eq.  (111.8c),  was  assumed  to  be  a  delta 
function.  We  shall  now  investigate  the  memory  effects  of  the  kernel  function  on 
the  energy-flow  dynamics  (which  up  to  now  were  assumed  to  be  governed  by  an 
irreversible  process  with  an  infinite  heat  bath).  To  discuss  the  types  of 
selectivity,  the  energy  population  of  a  multilevel  system  will  be  studied,  and 
finally,  the  excitations  of  the  a'bond  in  the  H/W  system  will  be  examined 
numerically. 
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a.  Memory  Function  and  Feedback  Mechanisms 

As  illustrated  in  Fig.  12,  in  a  simple  heat  bath  model,  the  system  consists 
of  the  active  (A)  mode  and  the  bath  (B)  modes  in  which  the  B  modes  are  treated 
as  an  infinite  energy  sink  due  to  their  large  number  of  modes  and  the  continuum 
phonon  spectrum.  The  energy  flow  from  A  to  B  is  irreversible  (with  a  relaxation 
rate  y)  within  the  spirit  of  Markoff  processes.  In  a  modified  model,  the  system 
is  divided  into  three  groups:  the  active  (A)  mode,  the  quasi-continuum  (B) 
modes  and  the  true  continuum  (C)  modes.  For  CO/Ni  as  an  example,  the  v3  mode 
(with  frequency  'vlOOO  cm"*)  is  the  active  mode,  and  the  adspecies  internal  and 
adbond  modes  (^300  cm"*)  are  the  B  modes.  We  shall  demonstrate  that  in  the 
A-B-C  model,  although  the  energy  flow  from  A  to  C  and  B  to  C  are  irreversible, 
the  energy  flow  between  A  and  B  is  not  necessarily  a  one-way  process.  Depending 
upon  the  kernel  function  and  the  density  of  states  associated  with  the  quasi - 
continuum  B  modes,  the  interaction  between  A  and  B  may  be  a  two-way  process, 
i.e.,  there  can  be  a  feedback  from  A  to  B.  Furthermore,  we  shall  show  that 
there  is  an  interference  effect  which  causes  an  indirect  coupling  between  A  and  B 
via  the  C  modes. 

The  microscopic  Hamiltonian  describing  the  energy  flow  dynamics  of  the 
A-B-C  system  subjected  to  IR  radiation  may  be  best  expressed  as  [Eq .  ( 1 1 1 . 1 ) 
with  i=A,  j=B  and  k=C,  where  the  C  modes  are  inactive] 

H  =  H0  +  HAB  +  HAC  +  HBC  +  HAF  +  HBB  +  HAA  +  HBF,  (III. 43) 

where  hq=^A+^B+HC  1S  the  unperturbed  Hamiltonian  of  the  system,  H. j( i .j^A.B.C) 
represent  the  interactions  among  them,  and  H^p  and  Hgp  represent  the  A-mode-laser 
and  B-mode-laser  interactions.  In  the  following  discussion,  we  shall  focus  on 
the  energy  flow  between  A  and  B,  and  the  last  two  terms  in  Eq.  (III. 43)  will  be 
ignored  since  our  model  system  will  consist  of  only  one  adspecies  and  one  active 
mode.  For  simplicity,  we  further  assume  single-phonon  processes  (the  extension 
to  multiphonon  processes  may  be  done  by  techniques  used  earlier). 
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The  equations  of  motion  for  the  operators  a»  bj  and  ck  for  the  A,  B  and  C 

14  84 

modes  are  given  in  the  rotating  frame  by  ’ 

•  sv  .  h  p  iV  iAt  (III. 44a) 

t  9j  bj  EJ  du>k  KAk  pk  ck  Ek  ‘  2  e 

A  a  i.  *1  i  fH  „  .  C  zr  ,  r  „  01Ajj,1:  (III. 44b) 

bj  “i9JaEj  d“k  Kjk  pk  ck  Ek/Ej 

,  JrJ 

ik--iit*kiEk‘l-ii|{jkbJEj/Ek  •  <in-44c) 

J 

Ei  =  exp(iA.t),  etc.  (III. 44d ) 

J  4 


V  is  given  by  Eq.  (III. 2f ) ;  g^  Kjk  and  1/^ ^ ,  represent  the  appropriate 

coupling  factors  associated  with  HAB>  HAC>  HgC,  and  Hgg,  respectively;  the 
detunings  are  defined  by  a^-oj,  Aj=t»^-Wj  a^^-oj^,  Akj=aij"“k*  etc‘  ’  w^ere 

and  oj^  are  the  frequencies  of  the  A  mode,  the  j-th  B  mode  and  the  k-th  C  mode 
with  the  continuum  spectrum  and  density  of  states  pfc.  Employing  the  many-body 
techniques  previously  discussed,  the  effects  due  to  the  C  modes  may  be  replaced 
by  the  frequency  shift  and  the  damping  factor,  whereby  Eq.  (I II. 44)  becomes 


iV  jAt 


5  .  -v-  I  IGj  bj  Ej  -  f  e 

J 


ii„,t 


'  -,GjaEj  '  *  |Kjol  °0]e  JJ  • 

J'  J 


(III. 45a) 


(III. 45b) 


G-i=91*+(G1+iG?)  has  two  components--g.j  resulting  from  the  direct  coupling  and 
(Gj+i^)  resulting  from  the  indirect  coupling  between  A  and  B  via  their  coupling 
to  the  C  modes,  where  Gj  2  «  KjOKAOpO’  and  YA  is  the  c_mocle‘‘induceCl  damping 
factor.  We  note  that  in  Eq  .  (III. 45b),  the  B-B  coupling  is  governed  by  the 
direct  coupling  factor  V..,  (with  j^j'),  the  C-mode-induced  indirect  coupling 

J  J 

tfjKj ,o|2po(with  j'^j)  and  the  damping  factor  yb=tt1  Kj , Q ) 2pQ(wi th  j'=j).  For 
tractable  results,  we  shall  keep  only  the  j'=j  term  in  Eq.  (III. 45b),  so  that 


-iA.t 

bj  s  'YBbj  "  iGjae 


(III.45c) 
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Eq . ( I II . 45)  has  the  same  functional  structure  as  Eq. ( III .6) •  except  that  now 

the  damping  factor  y^g  has  been  introduced  in  the  operators  for  the  A  and  B 

modes  due  to  their  energy  relaxation  to  the  C  modes.  Substituting  the  formal 

solution  of  Eq.  (III. 45c)  into  Eq.  (III. 45a),  we  obtain,  for  b.(0)=0, 

J 

*  =  "V  '  Jtdt'a(t')M(t-t')-1|  elAt  ,  (III. 46) 

where  the  kernel  (memory)  function  is  given  by 


M(t)  •  l  |Gj|zexp[(lAj-YB)t] 
J 


(III. 47) 


To  solve  for  a(t)  with  the  boundary  condition  a(0)=0,  we  obtain  the  Laplace 
transform  of  a(t)  from  Eq . ( III.46)  as 


a(s)  = 


-i(V/2)/(s-iA) 


(III. 48) 


(S+yA)+M(s) 

which  is  characterized  by  the  Laplace  transform  of  the  kernel  function,  M(s), 
given  from  Eq.  (Ill- 47)  by 


(III. 49) 


N 

H(s)  *  I  — !-li - 

j(snB)-Uj 


In  our  earlier  discussions  involving  the  Markoff  or  continuum  spectrum  approxi¬ 
mation,  M(t)  *  <S(t),  the  dynamics  of  the  B  mode  were  simply  damped  due  to  the  A 
mode.  To  investigate  non-Markoff  processes  or  reversible  processes  between 
A  and  B,  we  shall  now  evaluate  M(s)  without  assuming  a  continuum  spectrum. 

From  the  concept  of  the  energy-gap  law  and  the  multiphonon  relaxation  rate 
discussed  earlier,  we  know  that  the  coupling  strengt.i  |G.|2  is  strongly  decreasing 

J 

with  the  increasing  of  the  multiphonon  order  parameter.  For  only  finite  numbers 

of  modes  in  B,  we  then  may  truncate  the  sum  in  Eq .(III. 49 )  and  further  assume 

the  mean  value  A.=A  and  J  | G .  | 2  *  N|G|2  to  obtain 
J  j  J 


which  yields 


,  %  -1(V/2)(s+y  -iA) 

a(s)  = - 2 - 

(s-x1)(s-x2)(s-x3) 


9 


where  A3=iA  and  Xj  2  are  the  roots  of 


(III. 51a) 


s2+(WiA)s  +  YA(yB‘iA)  +  NiGi2  =  0  *  (III. 51b) 

We  shall  show  that  the  nautre  of  the  energy  flow  between  A  and  B  is  characterized 
by  these  roots.  For  near  resonance,  a«ya  g,  when  the  damping  of  the  B  mode 
is  in  the  range 

YA-2vff|G|  <  yb  <  YA+2*tflG|  ,  (III. 52) 

Eq . ( I II . 52)  has  two  complex  roots  (which  otherwise  are  real).  These  complex  roots 
yield  the  oscillatory  behavior  of  a(t),  and  accordingly  the  active-mode  exci¬ 
tation  <n(t)>  =  «a+(t)a(t)»  is  an  oscillating  function.  Therefore,  Eq .  ( 1 1 1 . 52 ) 
provides  the  criterion,  in  terms  of  the  damping  factors  ya  b  and  the  coupling 
strength  Yffi|G|,  for  the  energy  flow  from  A  to  B  and  back  to  A.  In  the  limit  of 
yA  b<<^GI  ’  we  readily  obtain  X12*±1’^l^l»  which  provides  the  frequency  of 
this  back-and  forth-flow.  We  note  that  in  the  single-mode  limit,  s  »  o, 

M(0)=N|G|  /Yg»  for  A  =  0,  which  gives  the  same  results  ac  that  of  the  Markoff 
processes  in  which  energy  is  irreversibly  flowing  from  A  to  B  and  C  with  a 
total  rate  of  ya+M(0).  The  active-mode  excitation  for  arbitrary  values  of 
yr  Is  given  by 


<n(t)>  =  | a(t) | 


(III. 53a) 


where  are  the  roots  of  Eq.  (III. 52),  X3  =  iA  and 


-i«/2)(yiA+xn) 

rIVV«VV) 


(III. 53b) 


with  n,n' ,n"=l,2,3. 

In  addition  to  the  above,  there  is  a  second  type  of  energy  feedback  based 

on  the  periodicity  of  the  kernel  function  in  the  quasi-continuum  and  the  memory 

103 

effects  of  the  bound-continuum  coupling,  and  a  third  type  which  we  now 
consider  corresponds  to  the  feedback  effects  resulting  from  the  direct  heating 
of  the  substrate  or  indirect  heating  via  phonon  coupling.  The  substrate 
temperature  obeys  the  heat  diffusion  equation,  given  in  a  one-dimensional 
approximation  as 


(III. 54) 


D  is  the  diffusivity  and  the  heating  source  is  S(z,t)  =  (a/pCy)(l-R)I(z,t) , 
where  a,  p  and  Cy  are  the  absorption  coefficient,  mass  density  and  specific 
heat  of  the  substrate,  respectively,  and  R  is  the  reflectivity.  We  note 
that  the  substrate  temperature  is  governed  by  the  volume  source  of  the  laser 
radiation,  I(z,t),  with  the  boundary  condition  DpCv(3T/3z)  *  0  on  the 
surface  (z  =  0).  The  solution  of  Eq.  (III. 54)  is  in  general  only  available 
numerically.  Except  for  the  situation  where  the  substrate  is  heated  by  a 
high-power  short  pulse,  e.g.,  gigawatt  picosecond  pulse,  the  above  diffusion 
equation,  may  be  replaced  by^ 


It  =  az  (D  W 


(III. 55a) 


and  the  boundary  condition 


k  (~)  =  -d-R)  1(0, t)  , 


(III. 55b) 


where  K  =  DpCy  is  the  thermal  conductivity  of  the  substrate. 

In  an  indirect  heating  process  for  the  aforementioned  A-B-C  system 
(Fig.  13),  the  boundary  condition  is  given  by 
[instead  of  Eq.  (III. 55b] 


K  (f£)  =  -F(t), 

32  z=0 


(III. 56) 


where  F(t)  is  the  rate  of  energy  flux  (W/cm2)  from  the  (A+B)  modes  to  the 
C  modes  whose  temperature  increases  by  an  indirect  channel  through  the  laser 
excitation  of  (A+B).  We  note  that  the  energy  flux  rate  F(t)  in  general  is 
time  dependent  since  the  excitation  and  relaxation  of  (A+B)  are  time  dependent 
However,  when  the  excitation  of  (A+B)  saturates  to  steady  state,  the  rate 
of  energy  flow  from  the  radition  into  (A+B)  equals  the  rate  of  eakage  from 
(A+B)  to  C.  In  this  situation,  the  steady-state  energy  flow  F  (t)  = 
ol/tlu),  where  a  is  the  steady-state  total  absorption  cross  section  of  (A+B) 
and  may  be  written  as  [see  Eq.  (III. 19)] 


.  f  Air 

s's-  (A-2e*X)2+(r/2)2 


(III. 57) 


where  f  is  a  fraction  relating  the  microscopic  number  of  modes  to  the  macro¬ 
scopic  thermal  parameter,  e.g.,  thermal  conductivity  of  the  substrate,  and 
r  is  the  damping  factor  of  (A+B). 
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In  general,  the  damping  factor  is  time  dependent  and  given  by  the 
thermally-averaged  golden  rule,  similar  to  Eq.  (III. 24), 


r(t)  =  (fj)Re/d,  «HA6c(t)  Habc(0): 


(III. 58) 


where  the  interaction  Hamiltonia  H^g(,  is  related  to  the  power  series  of  the 
correlation  function  of  the  C-mode  coordinate,  «Qr(t)  Q  (0)>>,  which  generally 
depends  on  time  through  the  time-dependent  occupation  number  <n^(t)>.  In  view 
of  the  fact  that  r(t),  <n^  g(t)>  and  <n^>  are  all  varying  in  time,  the  func¬ 
tional  form  of  F(t)  should  be  quite  complicated.  To  find  the  temperature  of 
the  C  modes,  one  should  simultaneously  solve  the  rate  equation  of  (A+B)  with 
a  time-dependent  damping  factor  r(t)  and  the  heat  diffusion  equation.  However, 
when  the  excitation  of  (A+B)  Is  saturated  and  when  r(t)  is  assumed  to  be  a 
constant  (or  given  by  its  mean  value),  the  maximum  surface  temperature  of  the 
substrate,  subject  to  the  boundary  condition  of  Eq.  (III. 56)  with  a  steady- 
state  flux  F  ,  may  be  obtained  by  a  Green's  function  technique  to  be105,106 


1  -  To  +  2Fs.s.lV'D)]*/pCv 


(III. 59) 


From  this  result  and  the  expression  of  F  in  Eq.  (III. 57),  we  can  draw 

s  •  s  • 

the  following  conclusions  regarding  the  temperature  behavior  of  the  substrate 
due  to  indirect  laser  heating:  (i)  in  contrast  to  direct  heating,  the  indi¬ 
rectly-heated  surface  temperature  is  sensitive  to  the  laser  frequency  through 
the  Lorentzion  of  Fe  _  vs  detuning;  (ii)  Ti  is  linearly  proportional 

to  the  laser  intensity  at  low  excitation,  whereas  T*  <*  I13  with  1  /3<p<  1 ,  in  genera 
(i i i ) the  increase  of  the  substrate  surface  temperature  will  be  significant 

only  when  the  diffusivity  is  small;  otherwise  the  thermal  energy  will  diffuse 

into  the  bulk;  (iv)  the  increase  of  T*  results  in  a  higher  occupation  number 


of  the  C  modes, which  in  turn  may  provide  energy  feedback  via  their  thermal 
phonons.  Furthermore,  the  increase  of  <ng(t)>  leads  to  the  decrease  of  the 
rate  of  energy  flow  from  (A+B)  to  C  given  by  r(t)[<n^g(t)>  -  <n^(t)>]  . 
However,  we  also  note  that  the  overall  effects  of  <ng(t)>  on  the  energy- 
flow  dynamics  of  the  system  are  not  yet  clear  since  the  relaxation  rate, 
r(t),  is  also  time  dependent.  The  above  results  based  on  a  square  laser 
pulse  tend  to  overestimate  the  surface  temperature  compared  with  that  for 
an  actual  pulse  which  is  close  to  a  Gaussian.  For  better  results,  we  may 
approximate  the  Gaussian  by  an  isosceles  triangle,  which  gives  T£  =  4/(27)* 
(=  0.77  of  that  for  the  square  pulse).1*^ 

To  further  investigate  the  nonequilibrium  behavior,  we  consider  the 
rate  equation 


d<nAB>  .  »I 


-  r(t)  [<n..(t)>  -<nc(t)>] 


(III. 60) 


which  has  the  same  structure  as  that  of  Eq.  (I II. 18)  except  that  the  time 
dependence  of  r(t)  and  <nc(t)>  are  given  by107 


Jt  i“nT  + 

r(t)  -J  dxe  0  «[lKJ.ct(T),  lK.C.(0)]»  , 
n  j  J 


(III. 61a) 


t  iw_T 


<nc(t)>  =Jdre  0  «I|Kj  |2Cj(0)Cj(t): 

n  J 


(III. 61b) 


where  wg  is  the  mean  frequency  of  the  (A+B)  subsystem,  is  the  multiphonon 

operator  for  the  C  modes,  and  K.  is  the  coupling  constant  between  (A+B)  and  C. 

J 

We  note  that  Eq.  (II I. 60)  reduces  to  the  equilibrium  expression,  Eq.  (II I. 18), 
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when  the  C-mode  correlation  time  is  much  shorter  than  the  inverse  relaxation 
time  and  the  upper  time  limit  of  Eq. (III.61 )  is  extended  to  infinity.  We 
finally  note  that  the  multiphonon  relaxation  rate  given  by  Eq ♦ ( I I I . 58)  is  a 
more  general  form  than  that  of  Eq.  ( 1 1 1 . 61 ) . 
b.  Multilevel  System  --  Selective  vs  Nonselective 

So  far  we  have  assumed  the  active  mode  to  be  a  harmonic  "oscillator" 
without  energy-level  structure.  For  a  true  quantum  system,  the  active  mode 
should  be  treated  as  a  multilevel  system.  As  shown  in  Fig.  13,  we  consider  a 
multilevel  system  in  which  the  active  mode  is  coupled  to  the  B  mode,  which  in 
turn  couples  to  the  continuum  C  modes  with  the  damping  factor  y.  In  the 
Heisenberg-Markoff  picture  (HMP)  and  within  the  rotating-wave  approximation, 
the  equations  of  motion  for  the  Bose  operators  are  given  by®^ 

dAn(t)  1 r 

■1“  *  ■«!  U(t)  “"t-W*!  •  Uii.ez) 

where  Hmn  are  the  matrix  elements  of  the  interaction  Hamiltonian 
H'  =  HAB+HBC+HAF(t)  g1ven  by 


Hmn  *  <*„!"’ I V  =  +k  «HAB+  HBC+  HAF>l**.*j.*k'>  . 


(III. 63) 


ABF 

being  the  eigenstates  of  the  A,  B  and  photon  modes  with  the  corres- 

o 

ponding  probabilities  Pm(t)s|Am(t)|  .  For  example,  in  a  two-photon  process 
the  eigenstates  j002>,  |200>  and  |020>  represent  A  and  B  both  in  their  ground 
states,  only  A  excited  and  only  B  excited,  respectively.  Considering  the  A 
mode  as  a  three-level  system  with  excited  states,  coupled  to  the  m-th  and 
2m-th  level  of  B  which  is  in  turn  decaying  to  C,  the  equations  of  motion  for 


the  Base  operators  —  al  2  3  for  the  ^  mode  and  B1  2  ^or  the  B  mode  ln 
the  HMP  are  given  from  Eq.  (III. 62)  as 


a^t)  =  -  l'Vj  a2(t)  exp(-iAt)  ,  (III. 64a) 

a2(t)  =  -1  Vj  a^t)  exp(iA1t)  -  igj  B1  (t)  exp(-iAjt)  -  iV2a3(t)  exp(-iA2t), 

(III. 64b) 

a3(t)  =  -  iV2  a2( t )  exp(-iA2t)  -  ig2  B2(t)  exp(-iA2t),  (III. 64c) 

B^(t)  *  -ig|  a2(t)  exp(iAjt)  -  (yx/2)  B^t),  (III. 64d ) 

B2(t)  -  -  ig|  a3(t)  exp(iA^t)  -  (y2/2)  B2(t)  ,  (III .64e) 


where  Aj  =  o>A  -  u 


A2  =  A1  "  2e*’  A1  =  m“e  ~  “A*  A2  ”  ma)B  +  2e*  "  UA’  e*  is 
anharmonicity  of  the  A  mode,  and  V,  g,  Y  are  the  pumping  rates,  coupling  factors 

and  the  damping  rates  for  the  related  levels,  respectively.  The  above  coupled 

equations  are  numerically  solved  for  Vj  *  V2  =  V,  g^  =  g2  =  g  and  Yj  =  y2  =  Y 

to  obtain  the  level  populations 


Pj  *  |  at(t)  | 2 ,  PA  ■  |  a2(t)  | 2  +  1  a3(t)  | 2  and  Pg  =  iB^t))2  +  |B2(t)|2.  (III. 65) 

PA  and  Pg  describe  the  population  dynamics  of  the  photon  energy  deposited  in 
the  A  and  B  modes,  respectively,  while  P^  =  1  -  ( )  describes  the 
population  loss  of  the  (A+B)  modes  and  represents  thermal  heating,  i.e.,  the 
portion  of  the  photon  energy  randomized  in  the  phonon  bath  C  modes.  The 
energy  populations  are  shown  in  Fig.  14(A)  for  selective  excitation  of  the 
A  mode  with  (V,g,Y)  =  (4, 0.1, 0.4)  and  in  Fig.  14(B)  for  nonselective  heating 
of  the  C  modes  with  (V.g.y)  =  (4,1,1).  We  see  that  for  fixed  laser  pumping 
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rates  the  selective  excitation  of  the  active  mode  (A)  requires  a  small  multi - 
phonon  coupling  factor,  g,  and  small  energy  leakage  rates,  y,  out  of  the  8 
mode,  while  appreciable  nonselective  thermal  heating  of  the  bath  modes  (C) 
is  achieved  when  the  coupling  factor  and  the  damping  rate  are  comparable  to 


the  pumping  rate. 

c.  Adbond  Excitation  —  Hydrogen/Tungsten 

In  the  heat-bath  treatments,  we  have  focused  on  the  excitation  of  the 

active  mode  and  treated  all  the  other  modes  as  a  heat  bath,  which  leads  to  a 

damping  factor  and  the  frequency  shift  for  the  active  mode.  To  investigate 

the  energy  flow  within  the  adspecies  or  the  adbond,  the  vibrational  modes 

should  be  rigorously  treated  rather  than  singling  out  the  active  mode  from 

the  bath  modes.  The  intramolecular  vibrational  relaxation  (IVR)  rate  of 

SFg  molecule  in  the  gas  phase  has  been  addressed  in  Part  A. 3  of  Section  II, 

and  now  we  shall  extend  that  study  to  the  excitation  of  a  heterogeneous  system 

namely  hydrogen  atom  adsorbed  on  a  tungsten  surface.  As  shown  in  Fig.  15, 

there  are  total  of  six  (3n,  with  rrZ)  possible  librational  modes  in  the  C2v 

point  group  bridge  site  adsorption  of  H/W( 100) .  These  consist  of  two-fold 

degenerative  pairs  of  lateral  modes  (v2  and  v3)  and  vertical  modes 

in-phase  and  out-of-phase  modes).  In  view  of  the  frequency  spectrum 

Vj  =  1038.6,  \>2  =  645.3  and  v3  =  1290.6  cm-1,  we  shall  consider  the 

2 

fourth-order  anharmonic  coupling  potential,  proportional  to  w*llc*1 

dominates  the  vibrational  relaxation  rate  with  the  internal  detuning 
Aj  *  2tt(2v1-v2-v3)  ,  where  Qj  2  3  are  the  norma1  coordinates  of  the  Vj  2  3 
modes.  By  the  concept  of  the  energy-gap  law  and  the  multiphonon  relaxation 
rate  discussed  earlier,  we  presume  that  the  fourth-order  anharmonic  coupling 
is  the  lowest  order  we  can  use  since  the  detuning  Aj  is  not  far  off  resonance. 

The  Hamiltonian  describing  the  fourth-order  anharmonic  coupling  among 
the  intra-bond  modes  of  H/W  in  which  the  mode  is  subjected  to  IR  radiation 


is  given  by 


3 

H=I  rfu.aja.  +  X  Hanh  + 'HVcos(a)t)(aJ'  +ax)  ,  (III. 66) 

i=l 

where  V  is  given  by  Eq.  (III. 2f ) ,  and  the  anharmonic  coupling  is  Hflnh  » 

(Qj+Q2+^3'  *  ^he  diagonal  terms  of  this  coupling  give  the  anharmonic  components 

of  the  modes,  whereas  the  off-diagonal  terms  provide  the  mode-mode 

o 

coupling  with  the  fourth-order  coupling  term  x'QjQ2Q3,  where  X1  =  12\. 

Referring  to  Eq.(III.41),  we  can  write  the  equation  of  motion  for  the  ampli¬ 


tude  functions  B.  as 

J 

mi  =  X'cJ  BJ  B22  +  2X'CB*  B2B3E+  “(V10/2)2C1Ef  ,  (III. 67a) 

fhB2  =  X'C2  BJ  B22  +  X'CB^BJ  E.  ,  (III. 67b) 

WB3  =  X'C*  Bj  B32  +  X,CB12B|  E_  ,  (III.67c) 

E+  =  exp(±iA^t),  Ef  =  exp(iA^t)  ,  (III. 67d) 

a|  =  2tt(2v|  -vj  -\>3  ),  A^  *  2irvj  -a)  ,  (III. 67e) 


_  2 

where  C  =  C^C-j  and  other  parameters  are  defined  as  in  Eq.  (11.41).  The 
above  coupled  equations  have  been  numerically  solved  by  the  fourth-order 
Runge-Kutta  method  (using  a  DEC-10  computer  connected  to  the  NCAR  graphics 
software).  The  dynamical  feature  of  the  excitations  which  show  the  energy 
flow  among  v,,  v,  and  the  active  mode  (v,)  are  illustrated  in  Fig.  16. 


3.  Photon  Energy  Population 


In  the  heat-bath  treatments,  we  have  shown  that  the  laser  radiation  may 
selectivity  excite  the  active  mode  while  the  bath  modes  are  kept  "cold"  if 
the  pumping  rate  is  greater  than  the  relaxation  rate;  otherwise,  there  is 
nonselective  thermal  heating  of  the  whole  adspecies  and/or  the  substrate 
surface.  We  have  also  demonstrated  that  in  the  H/W  system,  the  IVR  rate  could 
be  very  fast  (^10  s  )  when  the  internal  resonance  condition  is  not  within 
the  low-order  anharmonic  coupling.  For  systems  with  a  very  fast  IVR  rate, 
the  photon  energy  initially  deposited  in  the  active  mode  may  be  randomized 
among  all  the  adspecies  and  adbond  modes.  In  this  situation,  the  master  (rate) 
equation,  which  treats  all  the  vibrational  modes  in  the  subsystem,  e.g.,  (A+B) 
in  Fig.  13,  on  an  equal  footing,  is  more  appropriate  than  that  of  the  heat-bath 
treatments.  We  should  note  that,  from  the  concept  of  energy-gap  law,  the  sub¬ 
system  involved  in  the  master  equation  is  defined  by  those  modes  which  are 
strongly  coupled,  e.g.,  (A+B)  modes  in  Fig.  13,  and  usually  excludes  those 
low-frequency  phonon  modes  (C)  which  are  separated  far  in  frequency  from  the  A 
and  B  modes.  Therefore,  the  energy  space  of  the  master  equation  is  defined 
by  the  number  of  photons  absorbed  by  the  whole  adspecies,  or  (A+B)  modes, 
rather  than  the  vibrational  quantum  number  of  the  active  mode. 

Before  presenting  the  master  equation  which  involves  the  pumping 
(relaxation)  rates,  we  shall  first  calculate  the  associated  absorption 
(stimulated  emission)  cross  section  for  the  adspecies.  This  may  be  obtained 
by  a  generalization  of  the  total  cross  section  of  an  anharmonic  oscillator 
in  Eq.  (III. 19), 

A(n+l)rn 

n  U-2e*X)2+(r/2)2 


(III. 68) 
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where  the  n  dependence  of  the  transition  dipole  matrix  element 
un  n+1  *  <n+l,  ' x 1 n>  *  /(n+l)u0  j  has  been  included.  an  represents  the  cross 
section  for  a  transition  between  the  m-th  and  (m+l)-th  levels  with  level 
widths  rn  and  rn+jj  rn  3  (rn+rn+^ )/2  is  the  mean  width;  and  X  is  the  steady- 
state  excitation  of  the  active  mode.  Noting  that  both  F  and  X  in  general 
are  n-dependent,  e.g.,  rn=(n+l)rg  for  the  transition  pair  (n,  n+1)  and 
X  «  (n+1)  at  low  active-mode  excitations,  the  cross  section  may  be  written 
in  the  general  form 

on  =  (n+l)aa*  ,  (III. 69) 

where  and  a  are  correlated  parameters  depending  on  the  relative  magnitudes 

2  2 

of  the  anharmonicity  and  the  bandwidth.  For  example:  o=l  and  a*=rQ/[A  +r0]. 
for  e*=0  (harmonic  oscillator)  and  F  »r  (constant  bandwidth);  a=-l  and 
a*®rg/(2e*)2  for  an  anharmonic  oscillator,  Trn=r0«2e*X,  X=n+1  and  A=0;  a=0 
and  a*=A/rQ  for  e*=A=0  and  ?n  =  (n+l)rQ. 

The  master  equation  describing  the  photon  energy  population  in  the 
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energy  (n)  space  can  be  written  as 

~i  -  -0«n«»)[V(3„/Vl>Pnnl  *  <IVl*T“>CPn-r<Vl/9n>'>J  .  <nI-7°> 

2 

Pn  is  the  population  (adspecies/cm  )  of  the  level  of  energy  ntTw,  i.e.,  absorbing 
n  laser  quanta,  and  gp  is  the  degeneracy  of  the  n-th  level,  related  to  the 
number  of  vibrational  modes  S  in  the  adspecies  [or  (A+B)  modes]  by 


S-l 

» 


g„  -  g0(n+D 


(III. 71) 
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and  on  n-1  are  the  quantum  cross  sections  for  transitions  between  the  level 

pairs  (n,  n±l)  given  by  Eq.  (III. 69).  We  note  that,  as  discussed  in  Part  A. 2. a, 

the  transition  rates  or  the  cross  sections  are  in  general  time  dependent,  and 

the  above  master  equation  with  constant  cross  sections  is  the  limiting  case 

84 

of  the  generalized  master  equation 


Hvn<t-f>  V(t') 

0  n1 


onn.(t-f)  pn(t')]  ’  (HI. 72) 


which  reduces  to  Eq.  (III. 70)  when  we  assume  just  the  dipole  transition,  i.e., 
°n'n  =  °n6n'  n+1’  and  a  Marko^  process,  i.e.,  crn(t-t')  =  6 ( t-t ' ) .  We  note 
that  within  the  harmonic  model  the  dipole  transition  is  a  selection-rule-allowed 
process,  but  if  the  anharmonic  potential  is  considered,  the  high-order  transi¬ 
tions  with  on,n  =  on<5n,  n±j  (j*2)  are  also  allowed  . 

The  exact  solution  of  the  quanta!  master  equation  for  general  forms  of 
a„  and  9„  is  not  available.  However,  we  shall  discuss  two  limiting  cases  which 
are  physically  interesting  and  can  be  analytically  solved. 

(1)  a=S-l  (single-mode  harmonic  oscillator).  Eq.  (III. 70)  becomes 


-3?  =  -  (Ic*/JM[(n+l)Pn+1  +  nPn_1  -  (2n+l)Pn]  , 


(III. 73) 


which  has  the  solution,  with  the  initial  condition  Pn(t=0)=N0<5(t) , 


Pn(t)  =  N0Wn(t)/[l+W(t)]r 


(III. 74a) 


W(t)  =  0*<P  /"tfco  , 


(III. 74b) 


<j>  ■  I  dt  I  =  laser  fluence  (J/cm^) 


(III.74c) 
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The  corresponding  average  excitation  (quanta/adspecies)  is 

1 

<n(t)>  =  jr-  l  nP  -  W  «  •  (III. 75 

Hq  n 

(2)  a=0  (constant  cross  section).  S=l.  For  this  case  Eq.  (III. 70)  becomes 

dP 

-j-  .  -{Ia*/«o)[2Pn-Pn+1-Pn_1]  .  (111.76: 

With  the  same  initial  condition  as  in  case  (1)  and  by  using  the  recurrence  rela¬ 
tion  of  the  modified  Bessel  function,  we  obtain  the  population  function 

Pn(t)  =  NNQexp(-2W)In(2W)  ,  (III. 77] 

where  N  is  a  normalization  constant,  Ip  is  the  modified  Bessel  function,  and  W 
is  again  given  by  Eq.  (III. 74b).  The  corresponding  average  excitation  for  this 
population  is 


<n(t)>  *  2(W/ir)ls  « 


(III. 73] 


We  note  that  the  average  excitation  is  proportional  to  the  square  root  of 
the  laser  fluence,  whereas  in  case  (1)  it  is  linearly  proportional. 

To  investigate  the  population  function  and  the  associated  average  exci¬ 
tation  for  the  general  forms  of  o„  and  g  ,  we  shall  now  assume  that  P  , 
ard  gn  are  smooth  functions  in  n-space.  Within  this  continuum  assumption, 

which  is  a  good  approximation  when  n  is  large,  the  quanta!  master  equation, 

31 

Eq.  (III. 70),  is  then  converted  into  the  classical  diffusion  equation 


91 

(III. 79) 


Substituting  Eqs.  (III. 69)  and  (III. 71)  inco  Eq.  (III. 79),  a  particular 
solution  of  Eq.  (III. 79)  consistent  with  the  initial  condition  Pp(t=0)  a 
NQ5(t)  is 

pn(t)  =  NN0gnexp[-nB/(e2W)]  ,  (III. 80) 

N  *  (e/g0)(e2w)"s/s/r(s/B)  ,  (iii.sob) 


B  =  2  -  a 


(III.80c) 


where  W  is  again  given  by  Eq.  (III. 74b)  and  N  is  the  normalization  factor. 

The  average  excitation  for  this  classical  population  function  can  then  be 
calculated  as 

oo 

<n(t)>  =  N"1  fdn  P  n  «  (B2W)1/Br[(S+l)/8]r(S/B)  ,  (111.81) 

0  Jo  n 


1/6 

which  is  proportional  to  a  '  (since  W  <*  $)  and  consistent  with  the  quantal 
results  of  Eqs.  (IH.’S)  and  (III. 78)  for  B  =  1  (a=l)  and  B  =  2  (a=0).  It  is 
worth  noting  that  for  a  =  -1  (anharmonic  oscillator),  <n(t)>  *  4  in  this 
classical  diffusion  model,  whereas  the  steady-state  excitation  of  the  active 
mode  <n.(t)>  *  I^2  in  the  quantal  Heisenberg-Markoff  model  Eq.  (III. 20). 

Combining  Eqs.  (III.71),  (I 11.80)  and  (III.81),we  can  express  the  population 
function  in  terms  of  the  average  excitation  as 


(III. 82a) 


F(S)  =  r[(S+l)/B]r(S/8) 


(III. 82b) 
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For  a  comparison  with  the  above  classical  diffusion  model,  we  now  consider 
a  Boltzmann  thermal  distribution  P*( t)  which  is  characterized  by  the  effective 
temperature  and  the  quanta!  degeneracy  g*  [gn  in  Eq.  (III. 71)  is  a 
classical  degeneracy]  as  follows: 

H(t)  ■  NN09;e*P[-"fr»meff]  _  (III. 83a) 


where 


N  = 


l  P*/N 
n=0  n  J 


-1 


=  [1  -exp(+Ta /k r 


9n  -  (n+S-l) !/[n! (S-l) !]  , 

and  the  effective  temperature  is  defined  as  the  average  excitation 
energy  per  vibrational  mode  and  is  governed  by  energy  conservation  as 


(III. 83b) 
(III.83c) 


-tTcun  =  *ttu)<n(t)>/S  ,  (III. 84a) 

n  =  [ex-\SWkTeff)  -  l]"1  .  (III. 84b) 

For  a  multiphoton  process  (<n(t)>  »1),  we  obtain  the  high  effective  tempe¬ 
rature  limit  (kTg^>> ifo) ,  and  for  n»S,  Eqs.  (III. 83)  and  (III. 84)  reduce  to 

7)  »  kT  ^eff  =  -ttcj<n(t)>/S  ,  (III. 85a) 

N  =  (kTeff/tfui)S  ,  (III. 85b) 

<  »  nS'VtS-l)!  =  gn  . 


(III.85c) 
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which  is  identical  to  the  result  of  the  classical  diffusion  equation, 

Eq.  (III. 82a),  for  the  case  a  =  1  (harmonic  oscillator). 

Another  important  population  function,  the  Poisson  distribution,  which 

••  OO 

follows  from  a  solution  of  the  Schrodinger  equation,  has  the  form 

Pn  “  (<n>n/n! )exp(-<n>)  ,  (III. 87a) 


where  <n>  is  the  average  excitation  of  the  adspecies  assumed  to  be  harmonic 
with  equal  energy  spacings.  We  note  that  <n>  in  general  is  time  dependent  and 
is  given  by14 

AI 

<n(t)>  =  — - -  [l+e'rt-2e"rt/2cos(At)]  ,  (III. 87b) 

/r+(r/2)z 

where  A  is  a  constant  proportional  to  the  transition  dipole  of  the  adspecies 
and  r  is  the  overall  level  width.  The  above  result  is  identical  to  Eq.  (II I. 12) 
with  n  =  0;  however,  <n>  here  represents  the  photon  quanta  absorbed  by  the 
adspecies  (as  a  whole)  rather  than  that  of  the  active  mode.  Moreover, 

<n(t)>  “I  in  the  Poisson  distribution  which  is  quite  different  from  the 

S  •  5  • 

fluence  dependence  <n(t)>  «  <j>^  in  the  other  types  of  distributions  aforemen¬ 
tioned.  For  a  comparison  of  the  Poisson  function  and  those  obtained  from 
the  diffusion  equation  and  the  Boltzmann-type  distribution,  we  plot  the 
results  for  various  values  of  S  and  a  in  Fig.  17.  The  associated  desorp¬ 
tion  probabilities  P^,  shown  in  Fig.  18,  are  defined  by 
1  00 

PH  =  tt-  £  P  (quanta!  population),  (III. 88a) 

a  "0  n=n*  n 


(classical  population), 


(III. 88b) 


where  n*  is  the  desorption  threshold.  The  adspecies  will  be  desorbed 
from  the  solid  surface  through  the  channel  of  a  direct  bond  breaking,  indirect 
migration-induced  desorption  or  Fano-type  autodesorption  via  mode-mode  coupling 
within  the  adspecies.  We  note  that  with  a  knowledge  of  the  energy  population 
together  with  the  measured  quantity  n*,  we  should  be  able  to  evaluate  the 
desorption  probability.  The  reverse  procedure,  namely  knowing  the  measured  n* 
and  the  desorption  yield  (probability  and/or  rate  of  desorptiGn)  will  provide 
the  information  about  P  ,  i.e.,  how  good  are  the  calculated  populations  for 
various  systems.  More  discussion  on  the  dynamics  of  desorption  will  be  given 
later. 

4.  Isotope  Effects 

Infrared  isotope  separation  of  species  in  the  gas  phase  has  been  success- 

p 

fully  carried  out  both  theoretically  and  experimentally.  However,  the  isotope 
effects  of  species  adsorbed  on  a  solid  surface  have  only  been  recently  studied 
either  experimentally100,109  or  theoretically.86  Unlike  gas-phase  isotopes 
whose  absorption  cross  sections  are  governed  mainly  by  their  frequency  difference, 
the  collective  excitations  of  isotopic  adspecies  are  characterized  by  many  other 
surface-induced  effects:  (i)  direct  dipole-dipole  interactions  among  the 
identical  and  isotopic  adspecies;  (ii)  indirect  adspecies-adspecies  inter¬ 
action  due  to  their  common  coupling  to  the  substrate  phonons,  i.e.,  the 
phonon- mediated  interaction;  (iii)  frequency  shift  and  level  broadening  of 
the  adspecies  induced  by  such  effects  as  substrate  heating,  multiphonon 
relaxation  and  dephasing;  and  (iv)  dynamic  effects  such  as  migration  and 
desorption.  The  above  effects  all  influence  the  absorption  cross  sections 
of  the  isotope  adspecies  and  usually  smear  out  the  "frequency  separation"  of 

different  isotopes  when  they  are  co-adsorbed  on  a  solid  surface.  For  systems 
12  16  12  18 

like  CO  and  CO  adsorbed  on  a  copper  surface,  the  frequency  difference 

of  the  C-0  stretching  is  about  47  cm"1  and  their  dynamic  dipole-dipole 

109 

coupling  is  rather  large. 


Therefore,  it  is  more  difficult  to  observe  the 
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isotope  effects  for  adsorbed  CO  than  for  adsorbed  and  Dg,  which  have 
weaker  intermolecular  coupling  and  a  larger  frequency  difference.  We  shall 
suggest  below  some  possible  mechanisms  for  isotope  separation  of  adspecies 

QC 

using  a  model  based  on  a  previously-developed  theory. 

The  ensemble-averaged  equations  of  motion  for  a  system  consisting  of 
a  mixture  of  two  isotopic  adspecies  may  be  obtained  by  extending  those  for 
the  single  adspecies,  Eq.  (III. 34),  as  follows: 

ia(t)  =  [u>gff(t)-^]a(t)  +  Va/2  +  D'b(t)  ,  (III. 89a) 

ib(t)  =  ^ff(t)-u)  b(t)  +  Vg/2  +  D'a(t),  (III. 89b) 

<nA(t)>  =  -VA«Im[a(t)]»  -  2D«Im[a(t)b+(t)]»  -  y^  <n^(t)>  -  nc/2j ,  (III. 89c) 
<nB(t)>  =  -Vg«Im[b(t)]» +2D«Im[a(t)b+(t)]»  -  Yi|<ng(t)>  -  n^/2  ,  (III. 89d ) 

where  VA  and  Vg  are  the  pumping  rates  for  the  A  and  B  modes,  nc  is  the  Bose- 

Einstein  distribution  for  the  C  modes,  Im  denotes  the  imaginary  part,  and 
A  B 

weff  anc*  “eff  are  t*ie  e^ect^ve  frequencies  of  the  isotopic  adspecies  A  and  B 
whose  ensemble  averages  are  given  by 

<<a)eff>>  =  UA,B  "  2EA,B<nA,B^^  '  irA,B/2  ’  (III. 90a) 

rA,B  =  Yl’B  +  y2,B  ‘  (III. 90b) 

In  the  above  equations,  a  new  adspecies-adspecies  coupling  strength,  D1,  is 
introduced, 

O'  *  D-i^+iDg), 


(III. 91) 


which  consists  of  two  components,  the  direct  dipole-dipole  coupling  strength  D 
and  the  phonon-mediated  coupling  (Dj  +  iD2),  where  Dj  and  D2  are  given  by  the 
real  and  imaginary  parts,  respectively,  of  the  integral110 


L 

5^rJdt‘  b(t,)  I  "csjGj  )] 

n  J 


(III. 9 


For  the  Markoff  process  or  for  a  continuum  C-mode  spectrum,  the  real  part  of 
the  integral  reduces  to  the  simple  expression 


°1  =  7rnCp0°0°0 


(III. 9 


which  is  proportional  to  the  multi  phonon  occupation  number,  nr  »  n(n .+l)-im^  , 

J  J 

the  density  of  states  of  the  C-modes, Pq,  and  product  of  the  coupling  strengths 

A  B 

between  the  isotopic  adspecies  and  the  C-modes,  GqGq.  Here  we  have  approximated 
pjGjGj  ^  P0G0G0  * 

The  coupled  equations  [Eq.  (III. 89)],  which  are  highly  nonlinear  due  to 
the  anharmonic  corrections  2e*<nA  B(t)>  and  the  isotope  coupling  strength  D', 
can  only  be  solved  numerically.  However,  one  can  obtain  the  steady-state 
excitatations  for  the  weak-coupl ing  case,  D'=0,  as 


which  show  that  one  of  the  adspecies  may  be  selectively  excitated  without  sig¬ 
nificant  excitation  of  the  other  when  the  laser  frequency,  u,  is  tuned  to  one 
of  the  optimal  values,  i.e.,  =  «Ug^>>  -u  =  0,  for  either  adspecies  A  or  B. 

To  demonstrate  the  effect  of  the  coupling  strength,  D',  on  the  dynamics 
and  the  steady-state  excitations,  we  plot  the  numerical  solutions  of  Eq . ( 1 1 1 . 89) 
for  the  harmonic  case  (e^  B=0)  in  Fig.  19.  It  is  seen  that  <n^(t)>  is  higher  than 
<nn(t)>  for  aa  <  Ag,  where  A^  b  *  g-u,  with  D'*0  [Fig.  19(A)].  As  D' 
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increases,  both  excitations  decrease  [Fig.  19(B)].  Increasing  the  coupling 

D  A 

strength  to  the  transition  value,  i.e.,  D'  =  D*  3 (A  +a  )/2»  causes  the  steady- 
state  excitations  to  become  identical  [Fig.  9(C)].  For  large  coupling  strength, 
D‘>D*,  both  excitations  are  low  and  <ng(t)>  is  higher  than  <n^(t)>  (Fig.  19(D)]. 
These  numerical  results  for  the  steady-state  excitations  are  seen  to  be  in 
accord  with  analytical  results.  Such  results  can  be  obtained  by  defining  the 
"difference  excitation"  N_=  X  -  Y,  where  X  and  Y  are  the  steady-state  excitations 
of  the  adspecies  A  and  B,  respectively.  We  obtain,  from  Eq.  (III. 89),  for 

eA,B=0,  VAsVB=V’  Y1=Y1=Y1  and  rA=rB=r’ 


n  *  \rrn 


JV2D,>/[4yi(Z1+Z2  ’ 


(III. 95a) 


Zj  =  aaab  -  D'^  -  (r/2)‘ 


(III. 95b) 


Z2  =  rn+/2  , 


(III.95c) 


n±  =  ab  *  aa 


(III. 95d) 


*A,B  =  “A.B  "  “ 


(III. 95e) 


The  above  expression  for  the  steady-state  "difference  excitation"  N_  displays 
the  following  important  features:  (i)  isotopic  selectivity  increases  with 
decreasing  coupling  strength;  ( i i )  when  the  coupling  strength  reaches  the 
transition  value  D'  =  D*  =  n+/2,  there  is  zero  selectivity,  i.e.,  N_  =  0  as 
shown  in  Fig.  19(C). 

For  further  investigation  of  the  selectivity,  we  consider  the  reduced  selec¬ 


tivity  S  defined  by 


r 

J  8  J  dt<nft( t ) >/<ng ( t ) >  , 


(III. 96) 
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where  tp  Is  the  laser  pulse  duration,  and  <n^(t)>  and  <ng(t)>  are  the  excita¬ 
tions  of  the  isotopes  to  be  solved  numerically  from  Eq.  (III. 89).  The  numerical 
results  for  non-interacting  (D'=0)  and  interacting  ( D ' /0)  isotopes  are  shown 
in  Fig.  20,  which  reveals  the  following  features:  (i)  for  the  harmonic  case, 
with  g=0,  S=1  for  AA=|Ag|  [curve  (E)],  (ii)  for  the  anharmonic  cases, 

B^0,  a  higher  selectivity  is  expected  [curve  (B)]  than  that  of  curve  (E) 

due  to  the  nonlinear  feature  of  the  excitation;  (iii)  the  optimal  condition 
for  high  selectivity  is  that  the  laser  frequency  has  positive  detuning  for 
both  the  A  and  B  adspecies  [compare  curve  (C)  with  (E)  and  curve  (A)  with  (B)]; 
(iv)  curve  (E‘)  shows  the  effects  of  D'  on  the  selectivity  for  the  case  of 
eA  B=0  at  D’ =D*=(A^+Ag)/2,  where  Ssl  at  steady-state  as  predicted  by  the 
analytical  expression,  Eq.  (III. 95);  (v)  for  ej  g?<0,  D*  is  "blue  shifted" 
toward  "S=l  [curve  (F')];  (vi)  an  increase  in  D1  shows  a  decrease  in  S’  [curves 
(A)  to  (G1)],  and  S<1  when  D'<0*,  corresponding  to  the  situation  of  <nft>  <  <ng> 
[curve  (D)  of  Fig.  19  and  curves  (F1)  and  (G1)  in  Fig.  20].  We  finally  note 
that  for  a  non-interacting  system,  D'*0,the  greater  the  frequency  difference, 
AA-Ag,  the  higher  the  selectivity;  however,  for  an  interacting  system,  an 
increase  of  the  coupling  strength,  D1 ,  does  not  necessarily  result  in  a  decrease 
of  the  excitations.  Numerical  solutions  of  Eqs.  (III. 89  and  (III. 96)  show  that 
there  is  an  optimal  set  of  values  (D**,  e^  b)  which  yield  maximum  excitations; 
e.g.,  for  e^=1,  Va=10,  =  y*=2  and  A^=8-3»  D**  =  2.9  yields  the  maximum 

value  of  <nA>. 

For  further  demonstration  of  the  effects  of  the  coupling  strength,  we 
look  at  the  total  steady-state  excitation  N+*X+Y.  Some  results  are  shown  in 
Fig.  21. 


99 


5.  Dynamic  Effects 

We  shall  end  this  section  by  discussing  some  dynamic  effects  resulting 
from  the  excitation  of  the  adspecies  and/or  the  substrate.  These  effects 
are  associated  with  phenomena  such  as  migration,  recombination,  decomposition 
and  desorption.  Here  we  shall  focus  on  theoretical  aspects  of  migration  and 
desorption,  while  other  processes  will  be  looked  at  in  Section  IV. 
a.  Migration 

To  describe  the  individual  behavior  of  the  adspecies  occupied  at  different 
lattice  sites,  we  introduce  another  parameter,  namely  the  site  operator,  in 
addition  to  the  ladder  operators  for  the  normal  modes  in  the  dynamic  Hamiltonian 
given  by 


H(t)  *  HgtQj.Qj.-Qj)  ♦  I  .VtQj.Qj.—  QjJeJft.  +  HAF(t>  • 

K  j  K 


(III. 97) 


Hg  is  the  unperturbed  Hamiltonian  of  the  system  (with  normal  coordinates  Qj), 

Vkk ' (Qj )  1s  the  lattice"s1te"dePendent  interaction  potential  of  the  system,  and 
cfk  and  c^,  are  the  site  operators  of  the  Bloch  states  |k>  and  <k'|,  respectively, 
which  can  be  expressed  in  terms  of  Wannier  functions  in  the  site  representation 


t  1  1  * 

c t  =  |k>  3  4  I  e  | n> 


(III. 98a) 


ck.  =  <k' !  »  --  I  e 
K  JR  m 


-*■ ,  ->■ 

•  i  k  •  R_ 


(III. 98b) 


Taylor  expansion  of  the  interaction  potential  , 


i  QiQ?'+"''  (iu-99) 

U  l  '  J  A 


gives  us  the  general  forms  for  the  intramolecular  couplings,  p  being  the  order  of 


the  multi  phonon  processes.  Using  the  second-quantization  expression  of  Eq. 
(II I. 99)  and  the  Wannier  site  representation  in  Eq.  (I 11.98),  we  obtain,  from 
Eq.  (II 1. 97),  the  microscopic  model  Hamiltonian  as  follows:^ 


H(t)  =  +  H°  +  H°  +  HJ.  +  Hj  ♦  H2  +  H3  +  H4  +  HAF(t)  , 


ha  ■•*>*»••  “c  ■  E0  I  'k  • 

j  n 


(III. 100) 


(III. 101a) 


HC  *  I  *nJ<c»<bW  +  l  VTncn(>T*a)  +  l  Vncn(aT+a)<b?bj>'  (UMOlb) 
n,j  n  nj  v 


/-t 


i,  =  y  J  c!cl, 

1  .  nh  m  "  m  * 


nliSn  c" 


nCm(a++a)  * 

n  m 


Ho  *  y  y  c!tc  (bt+b.)  , 
3  j  X  m  n  m  J  J 


4  =  l  ^  Wmncncm^a++a^bj+bj^  +  [higher-order  terms] 


HAF(t)  “  _L  V(t)  Knm  cncm^a++a)  * 


(III. 101c) 


(Ill.lOld) 


(Ill.lOle) 


(Ill.lOlf) 


(Ill.lOlg) 


X,  Y,  Z,  J,  K  and  G  are  the  appropriate  coupling  parameters  proportional  to  the 

first  derivative  of  the  interaction  potential  (V^. ), while  W  is  proportional 

to  the  second  derivative  of  V^,.  The  above  dynamic  Hamiltonian  provides  us 

with  the  following  features:  (i)  the  ground  state  site  energy  of  the  adspecies 

En  is  perturbed  by  H',  which  includes  changes  in  site  n  due  to  direct  interac- 
u  C 

tions  with  the  lattice  and  due  to  active-mode  excitation,  as  well  as  an  indirect 
interaction  with  the  phonons  via  the  active  mode;  (ii)  the  terms  H^  and  H^ 
represent  parallel  motion  of  the  adspecies,  |J  |  in  H^  being  related  to  the 
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intersite  transition  probability  from  site  n  to  site  m  due  to  coherent  motion;  (iii) 
the  terms  and  H3  represent  the  perpendicular  vibrational-motion-induced  intersite 
transitions  due  to  the  active  mode  A  and  the  bath  modes  B  of  the  adspecies/ 
surface  system,  respectively;  (iv)  H^  is  the  A  and  B  mode  coupling-induced 
intersite  transition;  (v)  finally,  the  excitation  of  the  active  mode,  governed 
by  H^t),  provides  the  dominant  driving  force  for  intersite  migration  of  the 
adspecies. 

For  the  case  of  chemisorption  on  a  lattice  site,  the  equilibrium  position 
of  the  adspecies  is  shifted  due  to  the  distortion  of  the  lattice.  In  this  case 
the  adspecies-phonon  interaction  is  very  strong  and  perturbation  theory  cannot 
be  used.  It  is  posssible,  however,  to  use  a  canonical  transformation  to  go  to 
lattice-dressed  operators  A^,  A  ,  Bt,  B.  and  cj"c  which  take  into  account  the 
shifted  equilibrium  position  of  the  adspecies  and  the  lattice  distortion.  The 
transformed  total  Hamiltonian  is  quite  complicated,  and  for  tractable  results 
we  shall  investigate  the  single-phonon  case,  where  the  dressed  Hamiltonian  is 


H(t)  *  Hg  +  HAg  +  Hc  +  HAp(t)  , 

(III. 102a) 

"0  '*0  l  CnCn  +  «v/A  +  5  ’ 

n  J 

(III. 102b) 

Hfto  =  I  'tiWL  CnCn  [A^.+ABt]  , 

AB  ntf.j  nm  n  n  l  ^ 

(III.102c) 

H:  =  l  cfc  [a+B -+ABt| 

C  nj  n  n  j  jj 

n  9  j 

( III . 102d) 

SAF(t)  =  ^  CnCn  (A%A)  ' 

(III.  102e) 

0)^^^  =  -  2e*  A  A 

(III.  102f ) 

(III. 102f ) 
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Eq  is  the  distortion  energy,  W^m  and  Znj.  are  the  transformed  coupling  constants 
of  and  Znj,  respectively,  and  V(t)  is  the  transformed  laser-adspecies  coupling 
constant  of  V(t)Knm.  Employing  the  many-body  technique  described  in  the  previous 
sections,  we  obtain  the  equations  of  motion  for  the  ensemble  average  «<•••»> 
(over  the  A  and  B  modes  and  the  lattice  site  coordinates)  of  the  active  (A)  mode 
and  the  lattice  site  transition  probability: 

d<A> 


dt 


=  ■i<“eff(t)>  <A>  -  iV(t)  - 


['fl+Y2^H]/2 


<A> 


(III. 103a) 


(III. 103b) 


dt 


2y 


m 


+1]nB  +  (v1)  <^A» 


4$  ■  4^  -  (t 

(III".  103c) 


Here  y2  is  a  dephasing  factor, and  y^  is  the  phonon-coupling-induced  damping 
factor  given  by 

i  -  1 2 


T1  =  ^  tb/[4j  +  (v2)2]  • 

J 

4j  ■  •  “j  1 


(III. 104a) 

(III. 104b) 


whera  yB  is  the  decay  factor  of  the  phonon  (B)  modes  due  to  anharmonic  coupling. 
In  Eq.  (III. 104)  and  in  what  follows,  we  assume  only  nearest-neighbor  contribu¬ 
tions,  namely,  Zj*m  =  zjj  n+^  =  ZJ,  etc.  The  migration-induced  damping  factor 
ym  is  given  by 


,  2 


I  |wJ| 


fv2l 

COS 

M  - 

A^sin 

L 

+ 

[vB/2l 

I2 

i 


1.  ■  1,»  ■-  .  a 


(III. 105) 
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where  A..  and  d  are  the  wavelength  of  the  B-mode  vibration  and  the  lattice 


spacing  of  the  substrate,  respectively.  We  note  that  a  simplified  model 
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which  assumes  an  infinite  phonon  lifetime,  i.e.,  yb=0,  results  in  y1*ym=0  for 
single-phonon  processes  and  for  >  wj.  For  systems  with  great  energy  gap, 
(u’eff"UJj)>>V  damping  factors  for  the  single-phonon  processes  are  very 
small  and  multiphonon  processes  are  dominant.  Employing  techniques  similar 
to  those  used  for  obtaining  Eq.  (III. 35),  the  expressions  of  and  for 
multi  phonon  processes  may  be  obtained  by  setting  A .  equal  to  <w  ff>-£  w .  and 

J  CTT  .  J 

-  i  ^  \  7  J 

multiplying  the  coupling  parameters  |ZJ|  and  jWJ|  by  the  multiphonon  occu¬ 
pation  number  of  the  heat-bath  modes  in  Eqs.  (I II. 104)  and  (II 1. 105). 

From  the  coupled  equations,  Eq.  (III. 103),  we  can  calculate  the  lattice 
site  occupation  probabi1ity^Pn(tj^  ,  which  in  turn  gives  us  the  mean-square 
displacement  of  the  adspecies 


<r2<4  •  “2 1 


(111.106) 


and  the  diffusion  (migration)  coefficient 


0  =  lim 

t-KO 


l  2t 


(III. 107) 


The  site  probability  function  ^PnU)^  is  in  general  not  analytically  available 
due  to  the  time-dependent  excitation^n^(t^  which  is  nonlinearly  coupled 
[Eq.  (III. 103)].  For  tractable  results,  we  investigate  the  large  damping  case, 
Y1  2>>yM’  suc^  adspecies  reaches  its  steady-state  excitation 

X  s  which  is  governed  by  a  cubic  equation,  Eq.  (III. 20).  Using 

this  steady-state  excitation,  we  may  solve  Eq.  (III.103C)  to  obtain  the 
quasi -steady-state  site  probability. 
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^Pn(t)^  =  In(4Wt)e"4Wt  ,  (III. 108) 

where  I  is  a  modified  Bessel  function  and 
n 

W  =  [2nBX  +  X  +  nB|vM  ,  (III. 109) 

Thus,  from  Eq.  (III. 106)  and  using  the  recursion  relation  for  the  modified 
Bessel  functions,  we  obtain  the  mean-square  displacement,  which  in  turn  yie  ^ 
the  migration  coefficient 

0  =  4Wd2  .  UII.110) 

This  is  related  to  the  laser  intensity  by  a  power  law,  Ip,  l<p<3,  since  W  «  Ip; 
p*l  for  low  excitations  for  the  harmonic  case  (e*=0),  and  P=l/3  for  high 
excitations.  We  note  that  the  above  laser-enhanced  migration  constant,  D  *  W, 
is  governed  by  an  Arrhenius  form  for  the  high-temperature  limit,  kT>>fjoo^ , 

D  =  Dq  exp|-EA/kT]  ,  (III. 111a) 

EA«||wj|  +  |Zj|  jtfa)A  (III. 111b) 

is  the  "activation  energy"  for  migration.  We  note  that  the  above  Arrhenius  form 

for  the  migration  coefficient  resulting  from  single-phonon  processes  may  be 

replaced  by  a  non-Arrhenius  function  for  multiphonon  processes. 

The  master  equation  given  in  Eq.  (III. 103c)  involves  the  assumptions  of 

nearest-neighbor  contribution  and  constant  migration  rate.  Removing  these  two 
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assumptions,  we  have  the  generalized  master  equation 


( t-t ' )P(m,n,j ,t 1 ) 


WJJ  (t-f  )P(m'  ,n‘ , j  *  ,t ' )  -  WJ  3 
mm'.nn1  m'nun'n 

(III. 112) 

where  P(m,n,j,t)  is  the  site  probability  of  finding  the  adspecies  at  the  lattice 

site  (m,n)  and  at  the  j-th  vibrational  level  at  time  t,  wf£,  nn(t)  is  a  time- 

dependent  transition  rate  from  state  (m'.n'.j1)  to  (m,n,j),  and  W^  ._(t)  is 

m  m,n  n 

the  associated  reverse  rate.  We  note  that  in  the  above  general  form  the  transition 
rates  are  not  only  time  dependent  but  also  governed  by  both  the  site  coordinate 
(m,n)  and  the  vibrational  states  of  the  adspecies  (j)  which  are  activated  by 
the  laser  radiation.  Depending  on  the  migration  barrier  and  the  vibrational 
energy  of  the  adspecies,  there  are  at  least  three  types  of  dynamical  transitions. 

Type  (I):  horizontal  migration.  The  laser  induces  the  transition  from 
(m,n,j)  to(m',n',j)  with  a  change  in  lattice  site  but  leaving  the  vibrational 
state  unchanged.  This  is  the  situation  that  the  adspecies  is  vibrational ly 
excited  to  a  mobile  site  followed  by  a  horizontal  migration,  but  returning  to 
its  initial  vibrational  state  by  energy  relaxation  to  phonons.  Mathematically, 
this  process  may  be  described  by  and  in  Eq.  (III. 101). 

Type  (II):  vertical  transition.  This  involves  the  transition  from 
(m,n,j)  to  (m,n,j‘),  where  the  initial  j-th  vibrational  state  is  promoted  to 
the  higher  state  j 1 > j  when  the  adspecies  absorbs  photons,  but  the  energy 
level  is  still  not  high  enough  to  cause  intersite  transitions.  For  this  case  the 
migration  barrier  is  large,  and  the  description  involves  H£  in  Eq.  (III. 101b). 

Type  (III):  oblique  transition.  This  involves  simultaneous  intersite  and 
vibrational  state  transitions  and  may  be  described  by  and  H^p  in  Eq.  (III. 101). 

We  note  that  a  type  (I)  transition  usually  involves  a  physisorbed  state 
for  highly  mobile  adspecies  and  may  be  thermally  induced  by  either  direct  laser 
heating  of  indirect  phonon  excitation.  A  type  (II)  transition  may  cause  direct 
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desorption  when  enough  phonons  are  absorbed  by  the  adspecies.  The  last  type 
of  transition  may  cause  indirect  desorption  via  surface  rate  processes  such 
as  recombination  or  decomposition  or  by  the  assistance  of  phonons. 

We  finally  note  that  the  laser-induced  selective  surface  migration  of 

84 

adspecies  has  been  theoretically  proposed,  although  no  experimental  work 
has  been  reported.  In  the  model  system  of  Lin  and  George,84  the  moble 
adspecies  (usually  in  a  physisorbed  state)  could  migrate  in  a  preferential 
direction  which  has  a  lower  migration  energy  barrier  and/or  higher  effective 
temperature  [see  Eq.  (III. 111)].  This  laser- enhanced  selective  migration, 
either  direction-selective  or  adspecies-selective  (if  isotope  and/or  mixture 
adspecies  are  involved),  plays  an  essential  role  in  the  surface  rate  processes 
which  are  diffusion-limited, 
b.  Desorption 

Desorption  mechanisms  have  been  recently  studied,  both  theoreti¬ 
cally78,80"82,89"92  and  experimentally.94”101  First  of  all,  we  should  note  that 
none  of  the  theories  developed  as  far  nor  measured  data  are  able  to  completely 
describe  the  desorption  mechanisms  which  may  be  achieved  through  several 
channels.  Depending  upon  the  states  of  the  adspecies  (such  as  physisorbed  or 
chemisorbed,  adatom  or  admolecule,  and  the  potential  of  the  adbond)  and  the 
frequency  spectrum  of  the  system  and  the  laser  field,  as  shown  in  Fig.  22  the 
desorption  may  result  from  direct  laser  excitation  or  from  indirect  energy 
transfer  processes.  The  desorption  channels  from  type  (I)  to  type  (V)  will 
be  discussed  separately. 

Type  (I):  direct  desorption  via  active-mode  excitation.  This  type  of 

desorption  channel  has  been  investigated  by  Slutsky  and  George  in  a  harmonic 

78  80  82 

model,  Lin  and  George  in  an  anharmonic  model  and  Metiu  et  al  in  a 

on 

Horse  potential  model  [see  Fig.  22(1)].  In  these  models,  the  active  mode 
is  excited  by  multiphoton  processes,  while  energy  is  leaking  to  the  bath  modes 
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which  lead  to  a  damping  factor  or  level  width  of  the  active  mode.  When  enough 
photon  energy  is  absorbed  by  the  active  mode  [or  the  adbond  in  the  one-dimensional 
(ID)  Morse  potential  model],  the  adspecies  may  be  desorbed  from  the  surface  by 
overcoming  the  desorption  energy.  In  a  harmonic  model,  this  threshold  photon 
energy  (or  laser  intensity)  in  order  for  desorption  to  occur  has  been  under¬ 
estimated  by  the  harmonic  model,  where  the  absorption  cross  section  is  governed 
by  a  constant  detuning  A=w^-w,  and  has  been  overestimated  by  the  ID  Morse 
potential  model.  While  the  anharmonic  model  of  Lin  and  George  (AMLG)  should 
sometimes  provide  a  good  estimate  of  the  desorption  threshold  energy,  none  of 
the  above  treatments  will  be  appropriate  when  the  IVR  rate  of  the  system  is 
faster  than  or  comparable  to  the  laser  excitation  rate,  since  the  photon  energy 
is  more  likely  randomized  within  the  adspecies  and/or  substrate  phonons. 

Overcoming  the  anharmonicity  in  the  AMLG  model,  A,=A-2e*<n>,  or  in  the  Morse 
potential  model,  and  thus  reaching  the  continuum.is  quite  unlikely  to  occur  for 
this  fast  IVR  situation.  We  also  note  that  in  the  Morse  potential  model,  the 
important  concept  of  the  "quasi -continuum"  is  excluded.  Therefore, one  might 
look  for  another  type  of  desorption  channel,  such  as  the  following  suggestion. 

Type  (II):  direct  desorption  via  the  quasi-continuum.  The  concept  of  the 
quasi -continuum  and  the  mechanisms  for  level  broadening  have  been  proposed  by 
Lin  and  George.  As  shown  in  Fig.  22(11),  the  active  (A)  mode  is  pumped  by 
near-resonant  radiation  and  the  photon  energy  is  rapidly  randomized  within 
the  admolecule(M)  with  a  damping  rate  y  to  the  C  modes.  There  we  have  intro¬ 
duced  the  quasi-continuum  (QC)  in  the  surface  potential  M-C  bond.  We  note 
that  for  a  small  adspecies  (e.g.,  single  atom  or  diatomic  molecule),  the  QC 
may  not  be  achieved  when  one  assumes  an  ID  potential.  As  discussed  in  Part 
A. 3  of  this  section,  the  master  (rate)  equation  treatment  for  the  incoherent 
excitation  in  the  regime  of  the  QC  is  able  to  provide  a  quantitative  descrip- 
tion  of  the  desorption  mechanism.  We  should  also  note  that  the  level 
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broadening  of  the  adspecies  can  be  achieved  by  several  factors  such  as  the 
anharmonic  coupling  in  the  adbond,  the  surface-phonon-induced  Tj  and  l^  relaxa¬ 
tion  and  adspecies-adspecies  interactions  (static  dipole  coupling  or  dynamic 
scattering- induced  pressure  broadening).  All  of  these  broadening  effects  may 
achieve  the  QC  even  for  small  species  where  this  is  not  possible  in  the  gas 
phase. 

Type  (III):  indirect  desorption  via  tunneling.  Desorption  may  also  occur 

indirectly  when  the  adspecies  (particularly  for  a  physisorbed  state)  is  excited 

to  a  high  vibrational  level  whose  energy  is  degenerate  with  some  continuum  of 

the  adbond.  As  shown  in  Fig.  22(111),  the  active  mode  of  the  admolecule  is 

excited  with  its  excited  state  coupled  to  the  continuum  of  the  adbond  (B-C) 

to  cause  the  desorption  of  A-B.  This  type  of  desorption  channel  essentially 

involves  a  Fano-type  bound-continuum  coupling  and  has  been  employed  by  Lucas 

and  Ewing  in  their  predissociation  model  for  an  adspecies  and  later  by  Kreuzer 
92  91 

et  ar  and  Casassa  et  al .  Although  a  Morse  potential  has  been  used  for  the 

adbond  (B-C),  the  above  authors  have  assumed  either  a  single-photon  process  or 

a  multiphoton  process  but  with  a  harmonic  potential  for  the  active  bond  (B-C). 

This  harmonic  assumption  again  might  overestimate  the  transition  rate  (or  the 

92 

absorption  cross  section)  in  the  master  equation.  Furthermore,  a  ID  Morse 
potential  of  the  adbond  (B-C)  without  the  concept  of  the  QC  may  cause  errors  in 
the  estimation  of  the  desorption  rate. 

Type  (IV):  indirect  desorption  via  substrate  heating.  Shown  in  Fig.  22( I V) 
is  a  system  in  which  the  adspecies  is  not  active  to  the  radiation  but  the 
substrate  can  absorb  photon  energy  either  by  phonon  or  electron  excitation. 

The  adbond  may  be  broken  by  absorbing  thermal  energy  from  the  heated  substrate 
surface.  This  desorption  channel  may  occur  when  the  substrate  with  small 
diffusivity  is  highly  heated  and  the  adspecies  is  weakly  bounded,  e.g.,  physi¬ 
sorbed.  However,  this  type  of  desorption  is  much  less  sensitive  to  the  laser  fre¬ 
quency  than  types  (I)-(III)  where  the  dipole  transition  is  governed  by  the 


109 


detuning. 

Type  (V):  phonon-assisted  desorption  in  the  QC  regime.  We  propose  another 
type  of  desorption  channel,  as  shown  in  Fig.  22(V),  which  combines  those  mecha¬ 
nisms  described  for  types  (I)-(IV).  The  effective  adbond  potential  (M-C) 
consisting  of  regions  (i)  and  (ii)  is  coupled  to  the  substrate  phonons  which 
may  help  desorption  occur,  particularly  when  the  adbond  is  excited  to  region  (ii) 
and  near  the  true  continuum.  In  the  coherent  excitations  of  region  (i)  with 
strong  phonon  coupling,  the  energy  flow  from  the  adspecies  to  the  substrate  is 
more  likely  irreversible.  However,  when  the  adspecies  reaches  region  (ii)> 
i.e.,  the  QC,desorption  may  be  assisted  by  the  feedback  energy  from  the  thermal 
phonons.  When  the  adspecies  is  highly  excited  into  the  QC  region  and  in  thermal 
equilibrium  with  the  substrate,  the  desorption  rate  (kg)  may  be  approximated  by 

kQ  *  kQexp[-(E*-Wn>)/kBT*]  ,  (III. Ill) 

where  kg  is  the  pre-exponential  factor,  kg  is  Boltzmann's  constant,  and  <n>  is 
the  mean  number  of  photon  absorbed  by  the  adspecies  [given  by  e.g.,  Eq . ( I II . 75) 
or  (III. 78.)]  with  the  threshhold  desorption  energy  E*  and  maximum  surface 
temperature  T|  [given  by  Eq.  (III. 59)].  Another  way  of  describing  the  phonon- 
assisted  desorption  rate  is  by  the  expression 

kQ  =  k0exp[-E*/kB(Teff+T|)]  ,  (III. 114) 

where  Tg^  is  an  effective  temperature  of  the  adspecies  [given  by  Eq.  (III. 84)]. 

We  note  that  in  either  of  the  above  expressions,  Eq.  (III. 113)  or  (III. 114), 
the  desorption  rate  may  be  enhanced  by  the  laser  excitation  through  <n>  and 
Teff  and  further  assisted  by  the  thermal  phonons  via  T*  ,  where  T*  «  TQ 
(initial  surface  temperature)  when  the  substrate  surface  is  not  heated. 
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A  more  rigorous  expression  for  the  desorption  rate  given  now  as  time  depen¬ 
dent  is 

oo 

kD(t)  =  fdel  VI  (t)  Pn(t)  ,  (III. 115) 

0  n 

where  Pn ( t )  is  the  probability  that  the  adspecies  is  in  the  bound  vibrational 
state  |n>,  given  by,  e.g.,  a  Poisson  function  [Eq. (III. 87a)].  W  (t)  is  a  thermal 
transition  rate  from  the  bound  state  |n>  to  the  continuum  state  |e>  given  by 

Wne(t)  S  kQ  exp[-(e-En)/kBTs(t)]  ,  (III  116) 

where  En  is  the  energy  of  the  n-th  vibrational  state,  and  Ts(t)  is  the  time- 
dependent  surface  temperature  which  can  be  obtained  from  the  diffusion  equation, 

Eq.  (III. 55a),  with  the  boundary  condition  Eq.  (III. 56). 

Type  (VI):  indirect  desorption  by  dynamics,  The  final  desorption  channel 
to  be  discussed  is  also  an  indirect  type,  which  may  be  achieved  by  dynamic 
processes  such  as  migration,  scattering  and  reactions.  An  example  is  shown 
in  Fig.  22(VI),  where  the  Initially  chemisorbed  adspecies  may  make  a  transition 
to  the  mobile  physisorbed  P-state  by  absorbing  enough  photon  energy  to  over¬ 
come  the  transition  barrier,  but  still  not  have  enough  energy  to  overcome  the 
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desorption  barrier.  The  mobile  P-state  may  easily  migrate  on  the  surface 
and  eventually  desorb  from  the  substrate  by  further  absorbing  photon  energy 
via  the  direct  desorption  channel  [types  (I)  and  II)]  or  the  indirect  channel 
[types  (III)  and  (IV)]  or  the  phonon-assisted  channel  [type  (V)]. 

B.  Laser  Applications  to  Surface  Chemistry 

We  have  discussed  in  Part  A. 5  two  of  the  most  important  dynamical  processes 
associated  with  the  phenomenon  of  laser-stimulated  surface  processes  (LSSP) 
--migration  and  desorption.  In  addition  to  these  processes,  LSSP  in  general 
should  include  rate  processes  such  as  diffusion,  recombination,  decompo¬ 
sition,  dissociation  and  adsorption  (or  deposition).  Before  investigating 


the  experimental  aspects  of  LSSP,  we  display  in  Fig.  23  a  flow  chart  of  the 
methodology  of  LSSP.  Both  classical  and  quantum  treatments  have  been  developed 
for  the  absorption  cross  section  which,  with  the  master  equation  in  energy 
space, gives  the  energy  population  and  the  average  excitation.  Furthermore, 
by  solving  the  master  equation  in  the  lattice-site  space  we  may  investigate 
laser-induced  surface  migration  (which  in  the  field-free  case  has  been  studied 
by  random  walk  techniques. ^  We  have  shown  that  for  large  values  of  n,  the 
discrete  master  equation  reduces  to  a  classical  diffusion  equation  where  the 
energy  population  can  be  analytically  obtained.  The  energy  flow  and  population 
of  the  adspecies/surface  system  provide  the  fundamental  mechanisms  of  LSSP  as 
discussed  in  Part  A.  In  this  second  part  we  shall  study  the  applications  of 
LSSP  to  surface  chemistry,  which  includes  heterogeneous  catalysis,  chemical 
vapor  deposition  and  laser  annealing. 

Some  of  the  applications  of  lasers  to  surface  chemistry  and/or  physics  are 

(1)  Enhancement  of  surface  diffusion-limited  reactions. 

(2)  Enhancement  of  the  mobilities  of  selective  species  in  a  multi- 
component  environment. 

(3)  Control  of  the  concentration  of  reagents  by  selective  desorption 
or  excitation-induced  migration  of  the  species. 

(4)  Study  of  the  decomposition  and  recombination  rate  processes 
on  solid  surfaces. 

(5)  Study  of  the  catalytic  properties  and  heterogeneity  features 
of  the  adsorbents. 

(6)  Isotope  separation  and  mass  separation  of  adspecies  via  selective 
desorption  (laser  chromatography). 

(7)  Study  of  the  composition  and  location  of  the  active  sites  and 
the  conformation  structure  of  the  adspecies. 

(8)  Fabrication  of  microelectronics  via  laser-induced  chemical  vapor 
deposition;  microetching  and  laser  annealing. 
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We  first  discuss  laser  applications  to  heterogeneous  catalysis,  which 
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involve  the  following  processes:  (1)  adsorption  (physical  or  chemical) 

and  desorption  of  the  species  on  the  catalytic  surface;  (2)  migration  of 
adsorbed  species  and  subsequent  collisions;  (3)  interactions  (via  dipole- 
dipole,  charge  transfer,  etc.)  between  the  adspecies,  either  directly  or 
surface-mediated;  (4)  scattering  (reactive  or  nonreactive)  of  gas-phase 
species  by  the  clean  surface  or  adsorbed  species. 

A  laser  beam  might  influence  the  above  processes  in  a  number  of  ways 
depending  upon  the  physical  and  chemical  states  of  the  excited  species.  Some 
examples  are  (where  K  denotes  the  catalytic  substrate): 

(i)  Excitation  of  a  reactant  in  the  gas  phase 


A  +  K  -»•  A/K 


< 


— ►  AK  +  K 
-*  B  +  C  +  K 


(a) 

(b) 


(ii)  Excitation  of  an  intermediate  adsorbed  on  the  surface 


(A  +  B)/K  —  1— ►  (AB)/K 


- ►  AB  +  K 

- ►  C  +  D  +  K 


(c) 

(d) 


(iii)  Excitation  of  a  reaction  product  in  the  gas  phase 


A+B  +  K - ►  (AB)/K - ►  C  +  D  . 


(e) 


In  (i),  laser  radiation  may  enhance  the  adspecies-substrate  reaction 
(a)  and  the  decomposition  processes  catalyzed  by  the  substrate  (b).  Laser/ 
surface-catalyzed  decomposition  processes  have  been  investigated  for  the 

in  07 

chemisorption  of  CH^  on  Rh  and  HCOOH  on  Pt,  where  the  vibrational 
excitation  of  the  reactants  has  been  chosen  to  change  the  rate  of  adsorption 
and  the  amount  of  products  formed.  The  infrared  laser- induced  etching  of  a 
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semiconductor  has  been  studied  for  $Fg/Si,  where  the  gas-phase  reactant  SFg 

is  excited  and  then  chemsorbed  on  the  Si  surface  to  react  and  form  the  products 

SFg(g)  and  SiF^(g)  desorbed  from  the  substrate.100 

In  (if),  laser  radation  may  influence  the  overall  reaction  rate  in 

several  ways:  (1)  by  increasing  the  mobility  of  the  reactant  atoms  (A  or  B) 

through  photon  excitation  of  the  A-K  or  B-K  bond  with  subsequent  enhancement 

of  the  reaction  rate  kj  [see  type  (VI)  in  Fig.  22];  (2)  by  removal  of  the 

excess  energy  from  the  unstable  complex  (AB)  on  the  substrate  surface  via 

laser-stimulated  emission  accompanied  by  sv  face-phonon-mediated  relaxation, 

thereby  increasing  the  reaction  rate  kg;  and  (3)  by  breaking  the  AB-K  bond 

either  directly  through  laser  excitation  of  the  adspecies  or  indirectly  through 

thermal  desorption  by  laser  heating  of  the  surface  [see  types  (I)  and  (IV)  in 

Fig.  22].  We  note  that  the  direct  desorption  of  chemisorbed  species  from  a 

solid  surface  usually  requires  multiphoton  absorption,  necessitating  the  use 

of  high-power  radiation.  However,  much  lower  powers  may  be  sufficient  for 

the  desorption  of  a  diatomic  molecule  adsorbed  on  a  solid  surface  (A-B-K) 

if  the  photon  energy  absorbed  by  the  A-B  molecule  can  be  easily  transferred 

to  the  surface  to  break  the  B-K  bond  via  anharmonic  coupling  [see  type  (III) 

in  Fig.  22].  Laser  desorption  of  OH  radicals  from  a  silica  surface  has  been 

reported  in  which  a  low-power  (=10  W/cm  )  CW  COg  laser  was  used  to  excite  the 
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Si-OH  stretching  mode.  High-intensity  (=MW/cm  )  pulsed  COg  laser  desorption 

99 

processes  have  also  been  performed  for  CHgF/NaCl ,  CgHgN/Ag  and  the  isotopes 

CgHgN  and  CgDgN  coadsorbed  on  KC1. 100,101 

In  (iii),  we  illustrate  laser  excitation  of  a  reaction  product  after  it 

is  desorbed  from  the  substrate.  The  detection  of  the  OH  radical  desorbing 

from  a  Pt  and  Rh/Pt  catalyst  surface  by  laser-induced  fluorescence  has  been 
97 

recently  reported. 
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Another  important  aspect  of  laser  applications  is  in  laser-assisted 
chemical  vapor  deposition  (LCVD),  which  could  provide  a  new  method  in  the 
fabrication  of  microelectronics.  A  typical  sequence  of  events  occurring  in 
CVD  are: 

1.  Diffusional  transfer  of  the  gaseous  reactants  to  the  surface. 

2.  Adsorption  of  reactants  on  the  surface. 

3.  Events  on  the  surface,  e.g.,  reactions,  migration,  nucleation,  etc. 

4.  Desorption  of  products  from  the  surface. 

Laser  radiation  can  be  used  to  decompose  the  gaseous  molecules  or  heat  the 
substrate  in  order  to  enhance  adsorption.  Examples  of  systems  studied  so  far 
are:  the  decomposition  of  metal  alkyls  and  metal  carbonyls  to  deposit  on 
metals  and  the  decomposition  of  $iH4  to  deposit  on  silicon;****  reactions  of 
SiCl 4  with  H2  to  deposit  Si,  and  SiH4  with  N»3  to  deposit  Si3N4-116  We  note 
that  in  the  above  mentioned  examples,  laser  radiation  may  be  used  to  induce 
CVD  with  many  features  which  are  not  available  from  the  conventional  procedures, 
e.g.,  single-step  processes,  small-scale  deposition  with  high  homogeneity, 
remote  procedure,  and  selective  and  local  deposition. 

We  finally  mention  laser  annealing.  Both  CW  laser  solid-phase  recrystalli¬ 
zation  and  pulsed  laser  liquid-phase  epitoxy  have  been  reported.*  ^  Furthermore, 
the  controversy  regarding  the  physical  mechanisms  responsible  for  annealing, 

namely  the  thermal  melting  model  versus  nonthermal  plasma  model,  has  been 

118 

investigated  by  picosecond  laser  processes. 
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Fig.  1.  Power  law  for  the  average  number  of  photons  absorbed  per  molecule: 

(A)  fluence  dependence  <N>«^  and  (B)  intensity  dependence  <N>*I^.  Non¬ 
linear  effects  with  q<l  result  from  the  anharmonicity  of  the  absorption 
cross  section  [Eq.  (11.16)].  The  dots  are  the  experimental  data  based  on 
the  measurement  of  laser  f1uence*^a^  and  laser  intensity. 

Fig.  2.  Schematic  energy  diagram  for  infrared  multi  photon  excitation  of 
a  polyatomic  molecule.  The  discrete  states  of  the  active  mode  are  initially 
excited  and  coupled  to  the  quasi -continuum  which  may  be  further  pumped  to 
reach  the  true  continuum  by  absorbing  enough  photons.  Note  that  in  the 
quasi -continuum  the  energy  populations  are  usually  determined  by  a  set  of 
incoherent  rate  (master)  equations. 

Fig.  3.  Popluation  of  the  photon  energy  (Pn)  vs.  number  of  photons  absorbed 

per  molecule  (n)  in  the  quasi-continuum  region  characterized  by:  (A)  thermal 

distribution  [Eq.  (11.25)],  Poisson  distribution  in  a  harmonic  model  and  (C) 
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with  the  anharmonic  corrections. 

Fig.  4.  Red-shift  diagrams  of  the  optimal  detunings  for  infrared  excitation 
of  the  SFg  molecule:  (A)  maximum  active  (v^J-mode  excitation  (n*)  vs. 
external  detuning  (a£)  and  (B)  maximum  vg-mode  excitation  (n|)  vs.  internal 
detuning  (a{)  with  a  fixed  A^/Zir  *  2.7  cm”*. 

Fig.  5.  The  time-dependent  average  excitations  of  the  active  mode  (n^), 

the  \>2  mode  (r^)  and  the  Vg  mode  (n^)  for  the  detunings  (a^/2it,  aJ/2tt)  = 

(A)  (0,  0)  exact  resonance, (B)  (0,20)  far-off  internal  resonance  and  (C) 

(2.71,  1.00).  a1/2tt  =  2.71  is  the  optimal  external  detuning  for  the  SF, 

12  t  ° 

molecule. 

Fig.  6.  Time  evolution  of  the  average  excitation  of  SFg  subject  to  one  and 
to  two  infrared  lasers  at  various  scaled  [in  unit  of  (yj-h^)]  frequencies 
n..  =(A.jI.j  )1//2/(y1  +  y2)  and  scaled  detuning  Di  =  a./(y1+y2):  (A)  simultaneous 
excitation  by  two  lasers  with  (n,,Q?)  =  (0.5,  1)  and  single  laser  excitations 


with  (B)  (Q,  0)  =  (1,  8),  (C)  (fi,  D)  =  (0.5,  1)  and  (D)  (n,  D)  =  (1.5,  2). 

Here  the  time  scale  is  in  the  unit  of  laser  pulse  duration  (x)  with  intensity 

profiles  I.(t)  =  I.  sin^t/x),  i  *  1,  2.  47 
i  'o 

Fig.  7.  Two  radiative  mechanisms  for  curve-switching  ( j  =  +  j  -  j-)  in  an 

excimer  system.  Process  I:  during  the  collision  a  laser  LI  excites  a  bound 
state  of  the  system.  Subsequently,  a  stimulated  emission  brings  the  system 
into  a  new  channel.  Process  II:  while  in  the  initial  bound  state,  the  system 
absorbs  a  photon  from  a  second  laser,  L2,  and  enters  a  second  bound  state.  As 
in  I,  the  system  is  brought  into  the  final  channel  through  an  emission  stimu¬ 
lated  by  LI . 

Fig.  8.  Illustration  of  possible  collisional  conversions  involving  two  lasers. 
The  open  channels  |a>  and  je>  are  radiatively  coupled  to  the  bound  states  |s>, 

|t>,  and  |u>  by  the  first  laser.  These  bound  states  are  themselves  mutually 

coupled  by  the  second  laser,  which  induces  real  as  well  as  virtual  transitions 
(solid  and  dashed  arrows,  respectively).  Two  non-radiati ve  interactions  are 
also  included.  The  spin-orbit  interaction  couples  the  open  channels  |a>  and 
|8>  and  is  denoted  by  the  curved  arrow  in  (1);  the  second,  unspecified,  inter¬ 
action  couples  |t>  to  a  third  channel  |y>  and  is  depicted  in  (3). 


Fig.  9.  Schematic  representation  of  the  Hel2  molecule  in  the  X  and  B  elec¬ 
tronic  states.  The  motion  of  He  is  in  a  potential  defined  with  respect  to  the 
distance,  r,  from  the  He  atom  to  the  center  of  mass  of  I2>  The  other  coordinate, 
R,  is  the  separation  of  the  I  atoms.  In  the  absence  of  radiation  there  are 
certain  collisional  energies,  ,  for  which  the  system  undergoes  resonance 
scattering.  Here  the  motion  in  the  I2  bond  is  strongly  coupled  to  the  motion 
in  the  Hel  bond.  In  the  excited  excimer  state,  the  I-I"  bond  is  isoelectronic 
with  the  Xel  ground  state  (which  is  weakly  bound),  wherein  the  roles  of  strong 
bond  and  weak  bond  are  switched.  After  formation  of  the  vibrationally  predis- 
sociative  resonance  in  the  X  state,  a  photon  may  be  absorbed  exciting  the  system 
into  a  predissociati ve  excimer  state.  Subseguently ,  predissociation  occurs: 
vibrational  energy  from  the  I‘-He+  bond  is  transferred  to  the  weak  I-I"  bond 
resulting  in  bond  dissociation. 
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Fig.  10.  Radiative  scattering  pathways  with  initial  scattering  energy  E- . 

Top:  LI  (long  arrow)  effects  a  continuum-to-bound  transition.  The  level 
spacing  in  the  bound  state  manifold  is  such  that  a  subsequent  transition 
by  L2  (short  arrow)  can  be  made  to  a  higher  level),  with  subsequent  decay 
to  the  initial  scattering  state  but  at  an  increased  translational  energy 
Ef.  Bottom:  the  same  electronic  system  but  with  different  isotopic  compo¬ 
sition.  The  intermediate  level  spacing  does  not  permit  a  significant  transi¬ 
tion  by  L2,  therefore  only  elastic  scattering  is  permitted. 

Fig.  11.  The  normalized  steady-state  excitation  as  a  function  of  detuning  for 
the  harmonic  case  c*=  0  (solid  curve)  and  anharmonic  cases  with  e*=  critical 
value  (e**)  =  5/3  (dashed  curve)  and  e*=6  (dashed-dotted  curve).  The  bistable 
transition  points  are  shown  by  P,  Q,  R  and  S.  The  parameters  used  are 

(yi»Yo>a)  s  (1 ,10,2.5) ,  with  the  units  of  y,  in  10^s"\  y9,e*  and  a  in  cm-1, 

'  c  ?  4  1  _i  10-1 

I  in  W/cm  and  A  in  10  cm/W/s.  Note  that  1  cm  corresponds  to  3x10  s 

in  frequency. ^ 

Fig.  12.  Schematic  diagrams  for  the  heat-bath  models:  (I)  simple  case  with 
irreversible  energy  flow  from  A  to  B  and  (II)  modified  case  for  energy  flow  between 
A  and  B  while  being  damped  by  the  common  bath  C. 

Fig.  13.  Schematic  energy  level  diagrams  for  the  A,  B  and  C  modes,  where  V^ 

are  the  pumping  rates  between  the  i-th  and  the  (i+l)-th  vibrational  levels  of  the 

active  mode,  coupled  to  the  B  mode  via  multi  phonon  coupling  through  the 

factor  g^ ;  Ym  denotes  the  energy  relaxation  of  the  m-th  level  of  the  B  mode  due 

to  its  coupling  to  the  C  modes  which  are  condensed  modes  with  density  of  states 
79 

p. 


Fig.  14.  Energy  populations  of  the  A,  B  and  C  modes  for  two-photon  multiphonon 

processes,  given  by  Pft,  Pg  and  Pc,  respectively,  for  (pumping  rate,  coupling 

factor,  damping  rate)  =  (V,g,y)  =  (A)  (4,0. 1,0. 4)  selective  excitation  and 
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(B)  (4,1,1)  nonselective  thermal  heating. 

Fig.  15.  (A)  Hydrogenic  lattice  modes  of  H  chemisorbed  on  W(100)  at  saturation 

coverage  with  H  occupying  a  Cgv  point  group  symmetric  bridge  site.  Motions  of 
the  H  atoms  vertical  and  lateral  to  the  surface  are  indicated  by  the  arrows  t  and 
-►(*),  respectively.  (B)  Schematic  diagram  indicating  the  frequency  spectrum 


of  the  H/W(100)  system.  Note  that  the  IR  active  mode  (1048.6  cm  )  is 
coupled  to  the  v.  and  v,  modes  by  the  fourth-order  anharmonic  coupling 
x*QiQ2Q3. 

Fig,  16.  Dynamical  features  of  the  average  excitations  of  the  active  mode  (rfj), 
the  \>2  mode  (n^)  and  the  vg  mode  (rTg)  for  the  heterogeneous  system  H/W  with  detun¬ 
ings  (Ac/2it,  a}/2tt)  =  (A)  (0,0)  exact  resonance  and  (B)  (10,80).  Here  we 

*“  *  o  o 

have  used  the  following  parameters:  laser  intensity  I  =  10  W/cm  ,  anharmonic 
coupling  constant  A  =  5  cm"1  and  initial  temperature  TQ  =  300  K. 


Fig.  17.  The  distribution  functions  of  four-photon  excitations  (<n(t)>  =  4),  for 
a  Poisson  population  (--),  diffusion  model  population  (for  a  =  1)  with  S=1  ( — ) 

and  S=6  (••••)»  Boltzmann  population  with  S=1  ( - )  and  S=6  ( — ),  and  quantal 

population  with  S  *  a  =  1  (-•••-). 

Fig.  18.  The  desorption  probabilities  of  n*=5  [defined  by  Eq.(IV.88)]  associated 
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with  the  distribution  functions  shown  in  Fig.  17. 

Fig.  19.  Time -dependent  excitations  <n^  g(t)>  of  the  active  modes  for  the  harmonic 
case,  i.e.,  e*=0  with  (V,y,aa,ab)=  (10,1,4,8)  and  D' =  (A)  0,  (B)  2,  (C)  D*  and 
(D)  10.  D*=  (A^+Ag)/2=6  is  the  transition  value  where  N_=0.^ 

Fig.  20.  Time  evolution  of  the  reduced  selectivity  [Eq . ( I V . 96 ) J  for  the  non¬ 
interacting  [curves  (A)-  (B)]  and  interacting  isotopic  system  [curves  (B1)-  (G')]. 

The  parameters  used  are:  V  =  1 0 ,  y^,B=Y2’B=4,  e^=  eg=  e*  and  (O' ,e*,aa,ab) = 

(A) (0,e**,5,10) ,  (B)  (0,e**,5,-5) ,  (C)  (0,0,5,10),  (D)  (0,e**,5,0),  (E)  (0,0, 5, -5); 
and  (e*,AA,AB)=  (e**,5,10)  with  O'*  (B’)2,  (C')  4,  (D')  6,  (F1)  D*  and  (G' )  10. 
e**  =  1.28  is  the  critical  anharmonici ty  as  defined  in  Fig.  11,  and  0*  =  <W/2 15 
the  coupling  strength  at  the  transition  point  defined  by  Eq.  (Ill  .95 )  for  N_=0. 


Fig.  21 .  Total  steady-state  excitation  profiles  in  (N+,wA  B,D')-space  fc 
Yp«0.  V*  10  and  (A)  (15,5,5),  (B)  (15,7,3),’  (C)  (5,3,7) 

(D)  (5,3,7),  where  N  *  X  +  Y  is  the  total  steady-state  excitation,  X=<nfl> 

*  az  A  S 


for  e*=  0 
7)  and 


and  Y  =  <n_> 


B  s.s. 


given  by  the  numerical  solutions  of  Eq.(IV.89). 


Fig.  22.  Schematic  diagrams  of  adspecies-surface  systems  and  the  associated 
energy  leve’s,  where  A,  B  and  M  represent  the  adspecies  (adatom  or  admolecule), 

C  represents  the  substrate  (or  bath  modes),  and  the  laser  radiation  is  indicated 
by  the  wiggly  lines.  Several  types  of  desorption  channels  are  illustrated: 

(I)  direct  desorption  via  active-mode  excitation,  (II)  direct  desorption  via 
the  quasi-continuum,  (III)  indirect  desorption  via  tunneling,  (IV)  indirect 
desorption  via  substrate  heating,  (V)  phonon-assisted  desorption  and  (VI) 
indirect  desorption  via  dynamics. 

Fig.  23.  Flow  chart  of^the  methdd^ogy  of  LSSP^  indicating  both  (^[assjcal  and 
quantum  -treatments  of  tfte^bsorption  £r<*ss  section.^  By  solving  th^master 
equation  in  en£r§^space  and^fa^lattice-sit^  space,  we 'are  able  to  inVe^g^e 
the  ^ft&is^opulati&iv^desorption  and  migration,  which  are  essential  comp^rienls 
IK.  heterogeneous  rate  processes. 
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